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In  this  dissertation  techniques  for  dimension-driven  solid  modeling  have  been 
studied  with  the  goal  of  developing  reliable  and  efficient  approaches  and  algorithms  for 
design  intent  representation,  constraint  imposition  and  geometry  updating  in  dimension- 
driven  solid  modeling.  A  new  design  intent  representation  using  extended  CSG  tree  is 
presented.  Since  information  on  primitive  creation,  position  constraints  and  attach 
operation  is  integrated  into  extended  CSG  tree,  it  is  a  complete  representation  of  the 
design  intent.  The  major  contributions  of  this  dissertation  are  a  sequential  constraint 
imposition  method  and  a  non-regenerative  geometry  updating  algorithm  proposed  by  the 
author.  The  sequential  constraint  imposition  method  is  used  to  position  one  rigid  object 
with  respect  to  another  by  specifying  spatial  relations  or  position  constraints.  The 
constraints  are  specified  between  entities  of  these  two  objects  and  these  constraints  are 
imposed  sequentially.  Many  types  of  position  constraint  and  their  combinations  were 
studied  and  reduced  to  a  limited  set  of  constraint  types  and  combinations.  Based  on  this, 
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techniques  were  developed  to  impose  each  constraint  such  that  it  does  not  violate 
previously  imposed  constraints.  The  main  advantage  of  this  method  is  that  it  does  not 
require  solving  all  the  constraint  equations  used  in  the  model  simuUaneously.  By  using 
this  method,  constraints  can  be  imposed  faster  both  during  construction  and  updating  of 
the  solid  model  because  analytical  solutions  exist  for  most  constraint  combinations  when 
they  are  imposed  sequentially.  A  software  package  CONSPACK  was  implemented  based 
on  this  method  and  can  be  used  in  dimension-driven  solid  modeling  and  other 
applications.  The  non-regenerative  updating  algorithm  presented  in  this  dissertation  is  a 
novel  updating  algorithm  that  does  not  need  entity  recreation.  By  reusing  geometric 
entities  in  Boolean  attach  operation  and  recording  information  generated  in  Boolean 
operation,  the  algorithm  can  directly  modify  the  B-rep  of  a  solid  model  in  response  to  the 
dimension  change.  The  algorithm  is  robust  enough  to  handle  dimension  changes  that 
cause  a  change  in  the  topology  of  the  solid.  The  main  advantage  of  this  algorithm  is  the 
high  efficiency  achieved  by  avoiding  entity  recreation  and  persistent  naming.  The 
technique  for  localizing  the  geometry  modification  further  improves  the  efficiency  of 
geometry  updating.  A  feature-based  dimension-driven  solid  modeling  system  (FBM)  is 
implemented  to  test  these  methods  and  algorithms. 
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CHAPTER  1 
INTRODUCTION 

Overview 

Solid  modeling  is  the  fundamental  technology  used  in  all  modern  CAD  software 
systems  that  enables  designers  to  create  three-dimensional  models  for  visualization  of 
their  design  ideas.  A  solid  model  can  be  displayed  on  a  computer  screen  and  looks  very 
much  like  the  real  object,  and  can  be  manipulated  for  design  modification  by  editing  the 
mathematical  description  of  shape  and  its  dimension.  Based  on  this  visualization, 
designers  can  modify  their  design  many  times  until  a  satisfactory  design  is  obtained 
without  creating  many  physical  models. 

A  solid  model  is  a  complete,  valid  and  unambiguous  representation  of  a  three- 
dimensional  object  having  a  closed  volume.  It  contains  both  geometry  and  topology 
information  and  enables  any  point  in  three-dimensional  space  to  be  classified  as  inside, 
outside  or  on  the  object.  Based  on  the  data  structure  or  format  of  the  geometry  and 
topology  information  used,  we  can  classify  solid  models  into  many  representation  types. 
The  two  most  popular  are  Boundary  representation  (B-rep)  and  Constructive  Solid 
Geometry  (CSG)  representation.  B-rep  represents  a  solid  using  a  data  structure  containing 
boundary  entities  including  vertices,  edges  and  faces.  The  connectivity  information  or 
topology  of  these  entities  is  also  stored  in  the  data  structure.  CSG  representation 
represents  a  solid  using  a  tree-like  data  structure,  which  is  called  CSG  tree,  of  Boolean 
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operations  applied  on  primitives  used  to  construct  the  solid.  A  primitive  is  a  solid  with 
simple  shape  and  is  pre-defined  by  a  set  of  parameters. 

Since  a  B-rep  contains  a  large  amount  of  geometry  and  topology  information, 
direct  inputting  of  this  information  to  create  a  complex  solid  model  is  very  difficult  and 
therefore  this  approach  for  creating  B-rep  is  not  used  in  any  commercial  CAD  system. 
The  CSG  representation  can  represent  a  complex  solid  clearly  using  a  CSG  tree,  but  it  is 
usually  difficult  to  display  the  model  since  there  are  no  boundary  entities  in  CSG 
representation.  The  usual  way  to  create  B-rep  of  a  complex  solid  is  to  combine  many  pre- 
defined solid  models  together  using  Boolean  operations.  This  approach  is  called  dual 
representation  since  it  uses  both  CSG  tree  and  B-rep  data  structure.  The  Boolean 
operation  in  dual  representation  is  called  boundary  evaluation.  Some  details  about 
Boolean  operation  can  be  found  in  the  next  chapter. 

Dimension-driven  Editing  for  Solid  Model 

As  mentioned  before,  design  process  is  an  iterative  process.  New  design  is  often  a 
modified  version  of  existing  design.  The  final  design,  which  meets  new  design 
requirement,  can  usually  be  obtained  through  several  modifications.  When  using  a 
computer  aided  drafting  system,  designers  have  to  manually  modify  the  two-dimensional 
drawing  on  a  computer  screen  by  adding,  deleting  or  modifying  lines  and  curves,  and 
then  add  dimensional  annotations  to  match  the  changes.  In  traditional  solid  modeling 
systems,  the  designer  had  to  delete  and  recreate  the  model  to  make  changes.  Most  modem 
commercial  CAD  systems  use  dimension-driven  solid  modeling  technique  to  facilitate  the 
model  modification  and  reduce  product  design  period.  A  solid  model  can  be  considered 
to  be  dimension-driven  if  it  can  automatically  translate  changes  in  dimensions  to 
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corresponding  changes  in  geometry  (Gossard  et  al.,  1988).  This  translation  process  is 
usually  called  updating  or  regeneration  and  is  the  key  process  in  dimension-driven  editing 
of  solid  model.  For  performing  the  updating  process  to  obtain  a  modified  design,  the 
designer's  design  intent  must  be  recorded.  The  term  design  intent  is  used  here  to  refer  to 
specifications  regarding  how  the  designer  intends  the  design  to  be  updated  when  a 
parameter  or  dimension  is  modified. 

The  earliest  approach  to  implement  dimension-driven  editing  for  2D  drawings  is 
variational  geometry  (Lin  et  al,  1981;  Light  and  Gossard,  1982).  In  this  approach,  the 
design  intent  is  recorded  as  a  set  of  constraints  between  the  geometric  entities  of  the 
model.  The  geometry  of  the  model  is  updated  in  response  to  dimension  changes  such  that 
these  constraints  are  satisfied.  A  system  of  non-linear  simultaneous  equations  is  derived 
for  these  constraints  using  variables  that  describe  the  geometry  such  as  coordinates  of 
points,  radii  of  curves,  etc.  When  a  dimension  is  modified,  the  constraint  equations  are  re- 
solved to  obtain  new  values  of  the  geometric  variables  that  satisfy  both  the  specified 
dimensions  and  constraints.  The  geometry  of  this  model  is  then  updated  based  on  these 
new  values.  Variational  geometry  involves  solving  a  large  number  of  non-linear 
simultaneous  equations.  The  problems  in  this  approach  include  long  solution  time, 
muhiple  solutions,  over-constrained  and  under-constrained  equations  etc.  Many 
approaches  have  been  proposed  to  handle  these  problems  ranging  from  numerical 
methods  to  symbolic  manipulation  and  artificial  intelligence.  However,  direct  extensions 
of  these  methods  to  three-dimensional  solid  models  have  not  been  successful  because  the 
number  of  constraint  equations  required  to  fully  constrain  all  the  topological  entities  that 
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make  up  a  complex  solid  model  can  be  very  large.  This  aggravates  the  previously 
mentioned  problems  and  makes  equation  solving  slow  and  unreliable. 

A  practical  and  useful  approach  for  dimension-driven  solid  modeling  is  history- 
based  parametric  modeling  (Parametric  Technology,  1989).  Almost  all  modem 
commercial  CAD  systems  use  this  approach.  In  history-based  parametric  modeling,  the 
design  intent  is  represented  as  a  sequential  list  of  commands  or  modeling  operations 
including  primitive  creation  and  Boolean  operations.  When  the  dimensions  are  modified, 
the  history-based  parametric  modeling  system  will  find  the  first  operation  in  the  stored 
list  that  is  affected  by  the  dimension  change.  Then  the  system  will  retrieve  the  status  of 
the  modeling  process  before  that  operation  and  redo  that  operation  and  all  other  following 
operations  to  update  the  geometry  of  the  solid  model.  It  is  well  known  that  such  an 
approach  is  used  in  commercial  CAD  systems,  however,  there  are  few  papers  in  the 
literature  that  describe  this  algorithm  in  detail.  We  will  refer  to  this  updating  procedure  as 
"rollback  algorithm"  in  this  dissertation  because  the  procedure  is  very  much  like  rolling 
back  all  model  construction  commands  and  operations  affected  by  dimension  change  and 
redoing  them. 

It  is  very  popular  for  modem  CAD  systems  to  use  feature-based  implementation 
for  history-based  parametric  modeling.  The  term  "feature"  represents  a  high  level 
modeling  entity  and  has  a  different  meaning  in  different  systems.  In  feature-based 
systems,  the  operations  recorded  in  history  are  "feature  operations"  where  the  geometry 
of  the  feature  is  created  by  extmding,  revolving  or  sweeping  two-dimensional  profiles. 
The  feature  is  applied  by  attaching  its  geometry  to  prior  geometry  by  attach  operations 
such  as  Boolean  operations  or  merge  operations. 
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The  rollback  algorithm  has  some  obvious  disadvantages.  Firstly,  the  number  of 
operations  being  re-applied  depends  on  the  operation  that  is  affected  by  the  dimension 
change.  Therefore,  small  local  change  could  require  reconstructing  almost  the  entire  solid 
model.  Secondly,  when  the  operations  for  creating  primitive  solids  are  re-applied,  all  the 
geometric  and  topological  entities  of  these  solids  have  to  be  recreated  and  memory  for 
these  entities  has  to  be  relocated.  Thirdly,  when  the  re-applied  operation  is  a  Boolean 
operation,  redoing  it  is  very  time  consuming.  Fourthly,  since  the  result  of  updating  is  a 
new  B-rep,  the  correspondence  between  the  entities  before  and  after  modification  must  be 
created.  This  is  called  persistent  or  generic  naming  process  and  is  also  time  consuming. 
For  an  example,  one  operation  is  to  position  a  B-rep  A  wdth  respect  to  B-rep  B  by 
enforcing  one  planar  face  of  A, /a,  to  be  coincident  to  one  planar  face  of  B,  /b.  After 
updating  A  and  B,  they  are  now  new  B-reps  and  the  system  must  be  able  to  identify  the 
two  faces  corresponding  to /a  and  /b  to  reapply  the  positioning  operation.  The  purpose  of 
persistent  naming  process  is  to  create  and  find  this  kind  of  correspondence.  Finally,  in 
feature-based  implementation,  while  recreating  the  geometry  of  the  feature,  which  is 
generated  by  sweeping  a  profile,  variational  geometry  is  used  to  update  the  two- 
dimensional  profile  first.  Positioning  one  feature  with  respect  to  another  also  needs  to  use 
variational  geometry.  If  the  profile  is  very  complicated,  the  drawbacks  of  variational 
geometry  will  affect  the  efficiency  and  reliability  of  the  history-based  updating  process. 

Goals  and  Motivation 

Dimension-driven  solid  modeling  technique  takes  a  very  important  role  in  modern 
CAD  system  and  thus  has  attracted  many  researchers  working  on  this  area  in  recent 
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years.  However,  as  mentioned  in  last  section,  there  are  still  many  problems  in  the  existing 
methods.  Variational  geometry  method  has  many  fundamental  difficulties  when  applied 
to  3D  dimension-driven  solid  modeling  due  to  solving  large  number  of  simultaneous  non- 
linear equations.  These  difficulties  can  be  alleviated,  but  are  impossible  to  be  completely 
overcome.  History-based  parametric  modeling  is  very  valuable  in  practice  for  3D 
modeling.  Therefore,  the  objective  of  this  research  is  to  understand  how  dimension- 
driven  solid  modeling  works  and  try  to  solve  some  of  the  problems  listed  above  to 
improve  the  efficiency  and  reliability  of  history-based  parametric  modeling.  The  research 
is  mainly  focused  on  the  following  aspects: 

1 .  Develop  a  fast  updating  algorithm  for  dimension-driven  solid  modeling  that  does  not 
reconstruct  the  solid.  Instead,  it  modifies  the  geometry  and  the  topology  directly  to 
improve  the  updating  efficiency. 

2.  Develop  a  new  method  to  handle  position  constraint  problem  in  dimension-driven 
solid  modeling  to  avoid  using  variational  geometry  method. 

3.  Develop  a  new  method  to  represent  design  intent  in  dimension-driven  solid  modeling 
that  allows  us  to  achieve  goals  1  and  2. 

4.  Implement  a  feature-based  software  system  to  test  developed  methods  and 
algorithms. 

Scope  of  the  work 

In  this  dissertation,  a  design  intent  representation  using  extended  CSG  tree  for 
dimension-driven  solid  modeling  is  presented.  A  sequential  constraint  imposition 
algorithm  is  developed  for  sequentially  positioning  one  three-dimensional  object  with 
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respect  to  another.  This  algorithm  can  be  used  for  geometry  positioning  in  dimension- 
driven  soUd  modeling  to  avoid  using  variatonal  geometry  and  can  also  be  used  for  other 
applications  such  as  assembly  modeling.  A  new  updating  algorithm  is  proposed,  which 
does  not  require  recreation  of  the  geometry  and  topology  entities  in  updating  process  and 
thus  can  improve  the  updating  efficiency  and  avoid  generic  naming  problem.  To  test  the 
algorithms  and  methods  proposed  in  this  dissertation,  an  experimental  dimension-driven 
solid  modeling  system  is  implemented.  In  the  implementation,  the  sequential  constraint 
imposition  algorithm  is  implemented  as  a  stand-alone  package.  This  package  can  be  used 
for  any  other  applications  related  to  positioning  a  three-dimensional  rigid  body  with 
respect  to  another  by  specifying  the  geometric  relations  between  entities  of  the  two  rigid 
bodies  involved. 

In  chapter  2,  previous  work  related  to  dimension  driven  solid  modeling  is 
reviewed,  which  includes  the  history  and  background  of  solid  modeling,  variational 
geometry,  history-based  parametric  modeling  and  rigid  body  positioning. 

In  chapter  3,  we  introduce  a  new  design  intent  representation  using  the  extended 
CSG  tree.  The  definition  of  extended  CSG  tree  is  given  in  this  chapter  along  with  a 
detailed  explanation  of  the  new  information  added  into  traditional  CSG  tree  to  create  an 
extended  CSG  tree. 

In  chapter  4,  a  sequential  constraint  imposition  method  is  presented  for 
positioning  one  rigid  object  with  respect  to  another.  We  firstly  define  some  concepts  used 
in  this  method  and  then  the  methods  for  sequential  imposition  of  angular  constraints  and 
distance  constraints  are  introduced.  An  example  of  the  application  of  this  method  in  solid 
modeling  and  the  discussion  of  this  method  are  presented  in  the  last  part  of  this  chapter. 
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In  chapter  5,  we  introduce  a  non-regenerative  updating  algorithm  for  dimension- 
driven  solid  modeling.  We  firstly  discuss  geometric  entity  reuse  in  Boolean  operations 
followed  by  a  discussion  of  primitive  modification  and  positioning.  Then  a  detailed  non- 
regenerative  updating  algorithm  is  presented  for  cases  of  no-topology-change  and 
topology  change.  The  algorithm  is  further  improved  by  localizing  the  geometry 
modification  in  the  last  part  of  this  chapter. 

In  chapter  6,  implementation  details  of  the  proposed  approaches  and  algorithms 
are  introduced,  which  includes  definitions  of  major  classes,  data  structure  and  flowchart 
of  some  important  functions. 

Some  results  generated  in  the  implemented  system  are  presented  in  chapter  7  to 
show  the  user  interface  and  the  results  created  by  using  the  proposed  algorithms. 

Conclusions  about  this  research  are  made  in  chapter  8.  The  possible  extensions  of 
this  research  are  also  discussed  in  this  chapter. 


CHAPTER  2 
REVIEW  OF  PREVIOUS  WORK 

Overview 

In  this  chapter,  previous  work  related  to  this  research  has  been  reviewed.  Since 
soHd  modeling  is  the  basis  of  this  research,  the  history,  representations  and  Boolean 
operations  in  solid  modeling  are  reviewed  in  first  part  of  this  chapter.  The  techniques  for 
doing  dimension-driven  solid  modeling  can  be  classified  into  variational  geometry  and 
history-based  parametric  modeling.  The  approaches  for  variational  geometry  are 
reviewed  in  the  second  part  of  this  chapter,  which  include  the  concept  of  variational 
geometry  and  various  techniques  for  solving  geometric  constraint  system.  History-based 
parametric  modeling  is  introduced  in  the  third  part  of  this  chapter.  Since  it  is  very  popular 
to  implement  parametric  modeling  by  using  feature-based  modeling  approach  and  almost 
all  feature-based  modeling  systems  are  history-based,  feature-based  modeling  is  also 
introduced  in  this  part.  In  this  research,  sequentially  imposing  position  constraints  among 
features  is  a  key  step  in  solid  construction  and  regeneration  process.  Therefore,  previous 
work  related  to  positioning  a  rigid  body  by  spatial  relations  is  reviewed  in  the  final  part  of 
this  chapter. 

Solid  Modeling 

Visual  Models  in  CAD  System 

In  modem  engineering  design,  many  commercial  CAD  (Computer  Aided  Design) 

software  systems  are  used  to  help  engineers  to  visualize  their  design  ideas.  The 
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advantages  of  using  CAD  software  over  manual  drawing  in  design  are  obvious.  Firstly, 
design  visualization  using  CAD  software  is  much  faster  than  manual  drawing.  Secondly, 
modification  of  the  design  visualization  created  by  CAD  software  is  much  easier  than 
modifying  the  manual  drawing  since  the  visual  model  is  stored  as  electronic  information. 
Therefore,  using  CAD  software  system  greatly  facilitates  the  engineering  design  process. 

Early  CAD  software  systems  only  have  fiinctionality  of  two-dimension  drafting 
and  thus  are  also  called  Computer  Aided  Drafting  system.  Computer  Aided  Drafting  is  a 
direct  replacement  of  manual  drawing  and  enables  engineers  to  efficiently  create  and 
modify  standard  engineering  drawing  to  express  their  design  ideas  by  using  computers. 
Since  two-dimension  engineering  drawing  is  indispensable  in  industry,  Computer  Aided 
Drafting  is  still  an  important  part  in  modem  commercial  CAD  software  systems. 

In  the  design  process,  engineers  usually  create  physical  models  (prototypes)  of 
their  design  to  verify  and  analyze  their  design  idea  expressed  on  two-dimension  drafting. 
If  the  shape  or  dimension  of  the  physical  model  can  not  meet  the  design  requirement,  the 
two-dimension  drawing  is  modified  and  a  new  prototype  is  created.  This  process  is 
repeated  until  the  final  successfiil  design  is  obtained.  The  number  of  iterations  in  this 
design  process  is  determined  by  the  complexity  of  the  design  requirements  and  the 
experience  of  the  designers. 

To  avoid  this  costly  and  time  consuming  process,  the  geometric  modeling 
technique  was  developed  to  create  a  three-dimensional  visual  model  of  a  design  using  a 
mathematical  description.  The  visual  model  looks  identical  to  the  physical  model  on 
electronic  display  devices  but  is  much  easier  and  cheaper  to  create  and  manipulate.  Based 
on  the  mathematical  description,  designers  can  control  the  dimensions  of  the  visual 
model  and  easily  change  its  shape  by  cutting  and  adding  operations.  Therefore,  using  the 
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geometric  modeling  technique  enables  designers  to  directly  create  a  three-dimensional 
model  with  control  over  the  shape  and  dimensions  of  the  design,  greatly  reducing  the  cost 
and  number  of  iterations  in  the  traditional  design  process.  Due  to  these  advantages,  the 
solid  modeling  system  has  become  an  important  part  of  modern  commercial  CAD 
software  systems.  In  some  of  these  CAD  systems,  two-dimensional  engineering  drawing 
can  be  automatically  generated  by  projecting  the  visual  model.  With  this  functionality,  a 
three-dimensional  models  can  be  created  and  modified  first  and  then  an  engineering 
drawing  of  the  final  design  can  be  automatically  produced.  This  design  process  is  more 
natural  for  humans  because  a  three-dimensional  model  is  a  direct  visualization  of  design 
ideas  while  two-dimension  drawings  are  orthogonal  projects,  which  only  well  trained 
people  can  understand. 

Geometric  modeling  systems  can  be  classified  into  three  categories  based  on  the 
underlying  representation:  wire-frame  modeling  system,  surface  modeling  system  and 
solid  modeling  system.  A  wire-frame  modeling  system  represents  a  solid  using  lines, 
curves  and  points  to  describe  its  shape.  Wire-frame  model  is  a  natural  extension  of  the 
two-dimensional  drawing  and  is  easy  to  create.  The  CPU  time  and  memory  needed  for 
creating  wire-frame  model  are  relatively  small.  However,  the  information  included  in  a 
wire-frame  model  is  not  enough  to  explicitly  represent  a  real  three-dimensional  object. 
Sometimes  it  produces  ambiguous  representation  and  strongly  relies  on  human's 
interpretation.  A  surface  modeling  system  represents  an  object  as  a  collection  of  surfaces 
model,  which  is  an  extension  of  a  wire  frame  model.  The  additional  information  on 
surfaces  connecting  the  object  edges  is  provided  in  a  surface  model.  Since  more 
information  is  included,  a  surface  model  is  less  ambiguous  than  a  wire  frame  model. 
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Solid  modeling  techniques  (Requicha  and  Voelcker,  1982)  have  been  developed 
since  the  late  1960s  and  early  1970s.  A  solid  model  represents  a  three  dimensional  shape 
having  a  closed  volume.  In  wire-frame  and  surface  models,  only  geometry  information  is 
included  and  the  models  are  defined  by  a  list  of  curve  equations,  surface  equations  and 
the  coordinates  of  points.  A  solid  model  is  a  complete,  valid  and  unambiguous 
representation  of  a  three-dimensional  object,  which  contains  both  geometry  and  topology 
information  and  enables  any  point  in  3 -dimensional  space  to  be  classified  as  inside, 
outside  or  on  the  object. 
Representations  in  Solid  Modeling 

In  solid  modeling,  a  mathematical  description  of  an  object  can  be  represented  in  a 

variety  of  formats.  Various  such  formats  or  representation  types  of  solid  model  have  been 
developed  since  solid  modeling  technique  emerged.  The  two  most  popular  are  Boundary 
representation  (B-rep)  and  Constructive  Solid  Geometry  (CSG).  Other  representations 
include  Sweep  representation.  Analytic  Solid  Modeling  (ASM),  Primitive  Instancing, 
Cell  decomposition.  Spatial  Occupancy  Enumeration  and  Octree  Encoding  (Zeid  1991, 
Lee  1999). 

B-rep  represents  a  solid  as  a  data  structure  composed  of  boundary  entities: 
vertices,  edges  and  faces.  The  information  about  how  these  entities  are  interconnected  is 
also  stored  in  B-rep  data  structure.  Each  face  consists  of  a  linked  list  of  edges  bounding 
this  face  and  the  information  of  the  geometry  of  the  face  (surface  equation).  Each  edge 
consists  of  two  end  vertices  and  the  information  of  the  geometry  of  the  edge  (curve 
equation).  Based  on  how  to  represent  the  direction  of  edges  and  orientation  of  faces,  the 
data  structure  of  B-rep  can  be  classified  into  two  types:  half-edge  data  structure  and 
winged-edge  data  structure.  In  our  research,  half-edge  data  structure  is  used  and  will  be 
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briefly  introduced.  Winged-edge  data  structure  can  be  found  in  the  references  (Lee,  1999; 
Shah  andMantyla,  1995). 

The  convention  for  determining  the  face's  orientation  is  that  if  one  walks  along 
the  edge's  direction,  the  interior  of  the  face  is  always  on  the  left  side.  A  half-edge  data 
structure  assumes  the  model  is  a  manifold  model.  Therefore,  each  edge  has  two  adjacent 
faces  and  thus  has  two  opposite  directions.  To  represent  these  two  directions,  one  edge  is 
split  into  two  half-edges  which  represent  two  opposite  directions.  In  some  geometric 
modeler  such  as  ACTS  (ACTS  1996),  the  half-edge  is  called  COEDGE.  In  addition  to 
half-edge,  an  auxiliary  data  structure  called  loop  is  also  used  in  half-edge  data  structure. 

Figure  2-1  shows  a  block  with  a  square  hole.  The  B-rep  of  this  solid  has  a  total  of 
10  faces.  For  clarity,  only  visible  topology  entities  are  shown  in  the  figure.  The  half-edge 
data  structure  of  this  B-rep  can  be  represented  as  face  table,  loop  table,  half-edge  table 
and  edge  table  as  follows  (faces  F4, Fioare  omitted) : 


V4 


Figure  2-1  Half-edge  data  structure  in  B-rep 


Table  2-1 :  Face  table  of  the  B-rep  example 


Face 

Fi 

F2 

F3 

F4,  Fio 

Loops 

Li-o 

L2.0 

L3-0 ,  L3.1 
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Table  2-2:  Loop  table  of  the  B-rep  example  (The  loops  on  invisible  face  are  omitted  in 
this  table):   


Loop 

Li-o 

L2-O 

L3-O 

L3-I 

Half-edges 

C3.0—>  C443— > 
Cl-0 

C4-l^C7-l^ 
Cio-o^  Cs-O^ 
C4-I 

C3-l->C6-0-> 
Cg-o^  C7-0-> 
C3-I 

C13.0— >Ci2-0~^ 
Cii-o— >  Ci4^— > 
C13-O 

A  half-edge  table  is  a  mapping  between  half-edges  and  edges: 


Ci-o  and  Ci-i  — >  Ci    i=l,...,  n 
where  Ci-o  and  Ci-i  are  two  half-edges  corresponding  to  edge  ei ,  n  is  total  number  of 
edges.  The  direction  of  each  half-edge  is  shown  in  Figure  2-1 . 


Table  2-3 :  Edge  table  of  the  B-re 

p  example(l 

rhe  invisible  edges  are  omitted  in  t 

lis  table) 

Edge 

ei 

62 

63 

64 

66 

67 

68 

Vertices 

V4^  Vi 

Vi  V2 

V2^  V3 

V3->  V4 

V2->  V5 

V3^  V6 

V4->  V8 

Edge 

610 

611 

612 

613 

614 

Vertices 

V5-^  V6 

V6->  V7 

V9  Vio 

ViO^  Vii 

Vii^  V12 

Vi2^  Vi3 

CSG  representation  represents  a  solid  as  a  history  of  Boolean  operations  applied 


on  primitives  used  to  construct  the  solid.  A  primitive  is  a  solid  with  simple  shape  such  as 
block,  cylinder  and  sphere,  and  is  pre-defined  by  a  set  of  parameters.  The  history  of 
Boolean  operations  is  usually  stored  as  a  binary  CSG  tree.  Figure  2-2  shows  an  example 


(a)  (b) 
Figure  2-2  CSG  representation  of  solid  (a)CSG  tree;  (b)  represented  solid. 
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The  nodes  in  CSG  tree  represent  regularized  Boolean  operations  and  rigid  body 
transformations.  In  Figure  2-2(a),  and  -*  stand  for  regularized  Boolean  union 
operation  and  Boolean  difference  operation  respectively.  CSG  representation  only 
contains  the  history  of  Boolean  operations  and  thus  it  is  concise  and  requires  small 
storage  to  store  the  solid  definition.  However,  since  CSG  representation  does  not  contain 
any  boundary  entity,  it  must  be  converted  internally  into  B-rep  to  display  the  model  in 
applications.  Therefore,  a  solid  modeler  having  only  CSG  representation  is  rare.  Actually, 
most  modem  CAD  systems  use  dual  representation.  These  systems  only  store  the  CSG 
tree  and  can  always  reevaluate  the  corresponding  B-rep. 

Boolean  Operations 

Solid  model  is  a  complete  and  unambiguous  representation  of  the  geometry  of 

three-dimensional  object.  However,  when  B-rep  is  used,  it  is  diflficuU  to  directly  enter  all 
information,  such  as  vertices,  edges  and  faces,  of  the  mathematical  description  of  the 
solid  model,  especially  for  a  complex  solid.  Therefore,  dual  representation  is  often  used, 
in  which  simple  primitives  with  defined  boundary  are  pre-defined  and  various  Boolean 
operations  (union,  difference  and  intersection)  are  used  to  combine  them  into  a  complex 
solid. 

No  matter  what  representations  (B-rep  or  CSG)  are  used,  the  two  solids  involved 
in  Boolean  operation  must  contain  their  boundary  information.  For  B-rep,  the  boundary 
information  is  explicitly  stored  as  a  collection  of  vertices,  edges  and  faces,  whereas  for 
CSG,  the  boundary  information  is  not  stored  and  must  be  computed.  Based  on  the 
boundary  information  of  the  two  solids.  Boolean  operation  can  evaluate  the  boundary  of 
resulting  solid  and  thus  the  resuh  of  Boolean  operation  is  always  a  B-rep.  Boolean 
operation  is  called  boundary  merging  in  B-rep  and  boundary  evaluation  in  CSG  and  dual 
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representation  (Requicha  and  Voelcker,  1985).  Boundary  evaluation  can  be  classified  as 
non-incremental  and  incremental  boundary  evaluation.  In  the  non-incremental  boundary 
evaluation,  only  the  boundary  of  final  resulting  solid  is  evaluated.  In  the  incremental 
boundary  evaluation,  the  boundaries  of  intermediate  sub-solids  are  also  evaluated. 
Incremental  boundary  evaluation  is  much  more  popular  in  modem  solid  modeling 
systems  than  non-incremental  boundary  evaluation  and  is  actually  the  same  as  the 
procedure  of  boundary  merging  in  B-rep. 

Since  topology  entities  (vertices,  edges,  faces)  in  solid  model  can  be  thought  of  as 
point  sets.  Boolean  operation  in  solid  modeling  is  very  similar  to  that  in  set  theory.  The 
difference  is  that  Boolean  operations  in  solid  modeling  use  regularized  version  of  set 
operation  to  ensure  solidity  and  avoid  so  called  "dangling"  edges  or  faces.  Therefore, 
Boolean  operation  in  solid  modeling  is  usually  called  Regularized  Boolean  operation  and 
symbol  "*"  is  used  to  specify  the  regularity.  Figure  2-3  shows  the  difference  between 
standard  set  intersection  and  regularized  intersection  of  two  solids  (only  cross-sections 
are  shown  in  figure). 


Figure  2-3  Standard  Boolean  intersection  versus  Regularized  Boolean  intersection 

The  algorithms  implementing  Boolean  operations  in  solid  modeling  are  very 
complicated  and  are  different  fi-om  one  system  to  another.  Many  different  algorithms  for 
implementing  Boolean  operations  have  been  proposed  in  the  literature  (Mortenson,  1985; 
Mantyla,  1988;  Hoffmann,  1989;  Lee  1999).  However,  all  these  algorithms  conceptually 
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include  three  steps:  (1)  intersection  computation;  (2)  set  membership  classification  and 
subdivision  computation;  (3)  construction  of  connectivity  for  B-rep  entities.  In  stepl,  all 
intersection  vertices  and  edges  are  computed.  In  step  2,  some  edges  and  faces  are  divided 
into  edge  segments  and  face  segments  by  using  intersection  vertices  and  edges. 
Afterwards,  set  membership  classification  technique  (Tilove  1980)  is  used  to  decide 
which  segments  belong  to  resuhing  solid  and  which  B-rep  entities  of  two  original  solids 
belong  to  the  resuhing  solid  and  collect  them  in  an  entity  set.  In  step  3,  connectivity 
among  entities  in  the  created  entity  set  is  constructed  to  form  B-rep  of  resulting  solid. 


Variational  Geometry 

The  approach  using  variational  geometry  to  edit  geometry  was  first  proposed  by 
Lin  et  al.  (1981)  and  Light  and  Gossard  (1982).  In  this  approach,  geometry  of  an  object  is 
defined  by  a  set  of  characteristic  points,  which  can  be  represented  as  a  geometry  vector, 
V,  containing  the  coordinates  of  all  points: 

V  =  {x„y„z„-,x„,y„,zj''  (2.1) 

where  n  is  the  number  of  characteristics  points.  Geometric  constraints,  which  are  defined 
by  either  explicit  dimension  such  as  distance  or  implicit  dimension  such  as  parallelism, 
can  be  described  by  a  system  of  simuhaneous  equations: 

F,{V,d)  =  0  i=\,2,...,m  (2.2) 

where  d  is  the  dimension  vector  and  m  is  the  number  of  constraints.  The  equations  in 
(2.2)  are  mostly  non-linear.  Figure  2. 1  shows  a  2D  example  of  geometric  constraints. 
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Figure  2-4:  Example  of  geometric  constraints 


Only  two  geometric  constraints  are  shown  in  Figure  2-4,  which  are  distance  between 
points  6  and  5  and  line  1-2  being  horizontal.  These  two  constraints  can  generate  two 
equations  as  following: 


In  the  method  of  Lin  et  al.  (1981)  and  Light  and  Gossard  (1982),  variables  in  Fcan  only 
be  coordinates  of  points.  Late  works  (Owen,  1991;  Buchanan  1993)  added  new  types  of 
geometry  and  other  variables  to  equation  (2.2). 

Equation  (2.2)  is  the  basis  of  variational  geometry  method.  When  dimensions  in  d 
are  changed,  the  non-linear  simultaneous  equations  in  (2.2)  are  solved  to  find  the  new 
geometry  vector  V,  which  defines  the  new  geometry.  Light  and  Gossard(1982)  used 
numerical  method  (Newlon-Raphson  method)  to  solve  the  non-linear  simultaneous 
equations.  In  each  iteration  of  this  numerical  process,  change  of  vector  V  is  found  by 
using  Jacobian  of  constraint  equations  in  (2.2).  In  most  cases,  the  Jacobian  is  a  sparse 


^1  -yi  =0 


(2.3) 
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matrix  and  thus  some  of  the  constraints  can  be  decoupled  to  significantly  reduce  the 
computation. 

The  main  disadvantages  of  Newton-Raphson  method  are  unreliability  and  too 
much  computation  time.  When  m  ^  n  or  m  =  nhut  there  are  some  redundant  constraints, 
the  Jacobian  could  be  singular  and  numerical  difficulty  rises.  Light  and  Gossard  (1982) 
suggested  to  use  modified  Doolittle  method  to  handle  this  problem  and  other  methods 
(Light  and  Gossard,  1983,  Serrano,  1987)  were  also  proposed  to  detect  and  handle  the 
over-constrained,  under-constrained  and  redundantly  constrained  problem.  Another 
unreliability  of  Newton-Raphson  method  is  that  the  solution  of  the  equations  is  heavily 
dependent  on  initial  values  of  variables  in  geometric  vector.  Therefore,  the  solution  is 
usually  good  only  in  case  where  small  changes  are  made  to  the  initial  dimensions,  A 
proper  dimension  scheme  is  very  helpful  to  alleviate  the  unreliability  problem  of  the 
numerical  method.  Serrano  (1991)  proposed  a  graph  theoretical  approach  to  the  automatic 
generation  of  geometrical  constraints  to  do  automatic  dimensioning  in  design. 

Afterwards,  most  of  the  research  efforts  on  variational  geometry  have  focused  on 
how  to  solve  constraint  equations  so  that  the  drawbacks  of  Newton-Raphson  method  can 
be  overcome.  Owen  (1991)  developed  a  purely  algebraic  algorithm  to  solve  constraint 
equations.  He  firstly  presented  a  graphical  representation  for  the  equations  and  then  broke 
the  graph  into  pieces  each  of  which  can  be  solved  quadrafically  and  then  combined  the 
solutions  together  to  get  the  final  solution  of  the  constraint  system.  He  also  used  Galois 
theory  of  equation  to  classify  all  geometric  configurations,  which  can  be  solved  by  this 
method.  An  important  aspect  of  this  method  is  that  it  is  only  applicable  to  geometric 
configurations  that  can  be  constructed  by  ruler  and  compass. 
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A  symbolic  approach  was  proposed  by  Buchanan  and  De  Pennington  (1993)  to 
solve  geometric  constraint  system.  In  this  approach,  geometric  constraints  are  translated 
into  a  system  of  polynomial  equations.  These  equations  are  solved  simultaneously  by  a 
computer  algebra  technique  called  Grobner  basis  method.  Grobner  basis  of  a  polynomial 
equation  system  is  independent  of  redundant  equations  in  the  system.  Therefore,  the 
biggest  advantage  of  this  approach  is  that  it  can  handle  over-constrained  problem. 
Moreover,  in  over-constrained  case,  if  the  system  is  inconsistent,  Grobner  basis  can 
easily  detect  this  inconsistency  and  thus  avoid  a  misleading  result,  which  is  possible 
when  using  Newton-Raphson  method.  If  the  system  is  over-constrained  and  consistent 
without  redundancy,  this  approach  can  reduce  the  number  of  multiple  solutions  and  get 
the  right  answer.  Compared  to  numerical  method,  this  approach  is  reliable  and  complete. 
However,  it  is  too  slow  and  memory-consuming  to  solve  practical  problem.  This 
drawback  is  characteristic  of  computer  algebra  methods.  The  proposed  approach  has  been 
implemented  as  a  Constraint  Definition  System  (CDS)  by  the  authors. 

After  Owen's  method,  many  other  graph  analysis  approaches  were  proposed. 
Bouma  et  al.  (1995)  developed  a  graph-based  geometric  constraint  solver.  Similar  to 
Owen's  method,  the  constraints  are  also  translated  to  graph  and  then  solved  by  graph 
reduction.  The  differences  are  in  the  techniques  for  doing  graph  reduction.  If  the 
geometry  is  constructed  by  ruler-and-compass  constructions,  the  rest  of  the  solving 
process  is  same  as  that  of  Owen's  method.  However,  unlike  Owen's  method,  some 
construction  steps  are  permitted,  which  are  not  ruler-and-compass  based.  In  this  case,  the 
roots  are  found  numerically.  For  handling  multiple-solution  problems  in  a  well- 
constrained  geometric  problem,  this  solver  can  determine  all  possible  solutions  and  then 
find  the  solution  that  is  the  intended  one  with  high  probability.  In  case  the  found  solution 
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is  not  the  intended  solution,  user  interaction  is  necessary  and  the  solver  can  be  redirected 
to  the  intended  one.  Fudos  and  Hoffmann  (1997)  extended  this  approach  such  that  it  can 
more  efficiently  handle  over-,  under-  and  well-constraint  problem.  They  also  derived  the 
worst-case  time  complexity  for  the  algorithms  of  constraint  graph  analysis. 

Gao  and  Chou  (1998b)  presented  a  symbolic  approach  for  geometric  constraint 
solving.  Instead  of  Grobner  basis  method  in  (Buchanan  and  De  Pennington,  1993),  they 
used  Wu-Ritt's  decomposition  algorithm  to  solve  geometric  constraint  systems.  They 
discussed  how  to  compute  robust  numerical  solutions  for  a  set  of  equations  and  how  to 
decide  whether  the  constraints  are  independent  and  whether  a  constraint  system  is  over- 
constrained. 

Although  graph-based  and  symbolic  approaches  have  been  popular  in  recent 
years,  research  on  numerical  method  for  solving  geometric  constraints  is  still  continuing. 
Ge  et  al.  (1999)  discussed  using  optimization  method  to  solve  constraint  equations.  In 
their  method,  non-linear  equations  in  (2.2)  can  be  converted  to  unconstrained  non-linear 
optimization  problem  as: 

m 

rmn  a(V,d)  =  Y^F,(V,df  (2.4) 

1=1 

The  objective  function  of  this  optimization  problem  is  the  sum  of  squares  of  all  constraint 
functions.  After  the  optimization  model  is  formed,  well-known  optimization  methods 
such  as  BFGS  method  can  be  used  to  find  the  global  minimum  value  of  ofV.d),  which  is 
zero  in  case  the  solution  for  equations  (2.2)  exists.  It  is  obvious  that  there  is  no 
requirement  for  m  =  «  in  optimization  problem.  Therefore,  over-  and  under-constrained 
problem  can  be  handled  naturally. 


22 

Besides  numerical,  graph-based  and  symbolic  approaches  for  solving  geometric 
constraints,  another  type  of  method  is  called  propagation  approach  (Aldefeld,  1988, 
Shimizu  and  Numao,  1997;  Gao  and  Chou,  1998a).  In  this  type  of  approaches,  unknown 
variables  are  derived  from  known  ones  using  a  set  of  pre-defined  rules.  This  kind  of 
method  is  usually  implemented  using  techniques  of  expert  systems  or  logic  programming. 

Most  of  above  mentioned  approaches  are  focused  on  solving  2D-constraint 
system.  Although  some  of  them  (Shimizu  and  Numao,  1997)  claimed  their  method  can 
solve  3D  constraint  problem,  their  method  is  actually  combined  with  feature-based 
modeling  technique.  Direct  application  of  variational  geometry  method  to  3D  dimension- 
driven  modeling  is  not  successful.  This  is  because  the  number  of  the  non-linear 
simuhaneous  equations  is  very  large  and  the  solution  time  increases  dramatically. 
Moreover,  the  problems,  such  as  over-  or  under-constrained,  inconsistency  and  multi- 
solution  problem  become  more  serious  when  the  number  of  equations  is  large. 

In  the  reviewed  literature,  a  variety  of  geometric  constraints  have  been  used  for 
the  purpose  of  implementing  or  illustrating  the  constraint  solving  methods.  There  is  no 
standard  for  geometric  constraint  types  in  CAD  community.  To  enable  data  exchange  of 
different  constraint  solver,  Bettig  and  Shah  (1999)  derived  a  consistent  and 
comprehensive  set  of  geometric  constraints  for  shape  definition  in  Computer-Aided 
Design. 

History  Based  Parametric  Modeling 

History-based  parametric  modeling  is  a  more  practically  usefiil  method  than 
variational  geometry  method  in  modern  CAD  systems.  In  this  method,  a  sequential  list  of 
commands  or  modeling  operations  including  primitive  creation  and  Boolean  operations  is 
recorded.  When  the  dimensions  are  modified,  the  history-based  parametric  modeling 
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system  will  find  the  first  operation  affected  by  the  dimension  change  on  the  stored  list. 
Then  the  system  will  retrieve  the  status  of  the  modeling  process  before  that  operation  and 
redo  that  operation  and  all  other  following  operations  to  update  the  geometry  of  the  solid 
model. 

As  mentioned  before,  in  implementation,  history-based  parametric  modeling  is 
closely  related  to  feature-based  solid  modeling.  The  design  intent  recorded  in  model 
creation  history  is  usually  a  series  of  feature  operations  and  geometry  of  each  primitive  is 
parametrized  and  stored  in  feature  corresponding  to  the  primitive.  Therefore,  we  will 
firstly  review  previous  works  in  feature-based  modeling,  which  is  related  to  our  research, 
in  this  section  and  then  we  will  discuss  the  current  methods  and  problems  in  geometric 
updating  (or  regeneration)  in  history-based  parametric  modeling. 

Feature-based  Solid  Modeling 

In  solid  modeling,  B-rep  represents  a  solid  in  terms  of  low-level  entities  such  as 
faces,  edges  and  vertices,  which  are  difficult  for  engineering  designer  to  directly  create. 
CSG  representation  represents  a  solid  at  a  higher  level  of  abstraction  than  B-rep,  but  the 
primitives  and  CSG  operations  do  not  have  engineering  meaning  for  industry  designers. 
As  a  solution  to  this  problem,  feature-based  modeling  is  an  approach  where  higher-level 
entity  called  "feature"  is  used  to  provide  high-level  abstraction  of  the  part  that  directly 
relates  to  certain  design  functionalities  or  manufacturing  characteristics  (Shah,  J.  and 
Mantyla,  1995).  Feature-based  modeling  covers  many  topics  and  attracts  a  lot  of 
researchers  in  solid  modeling  area.  It  is  not  appropriate  to  discuss  all  previous  work  in 
feature-based  modeling  in  this  section  and  we  will  only  review  research  in  feature 
representation,  feature  operation  and  architecture  of  feature-based  modeling  system. 
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which  are  related  our  research.  Other  topics  such  as,  feature  recognition,  can  be  found  in 
the  Hterature  (Laakko  and  Mantyla,  1993;  Martino  et  al,  1994,  Han  and  Requicha,  1997). 

There  is  no  universally  accepted  definition  for  the  term  "feature".  One  definition, 
for  instance,  is  that  a  feature  represents  the  engineering  meaning  or  significance  of  the 
geometry  ofa  part  or  assembly  (Shah  and  Mantyla,  1995  ).  However,  it  is  safe  to  say 
that  a  feature  is  a  high-level  entity,  which  includes  its  geometry  and  other  design  or 
manufacturing  information,  and  the  relations  between  this  feature  and  other  features.  This 
kind  of  high-level  entity  can  be  seen  as  a  bridge  between  the  engineering  designer  and  the 
internal  representation  of  solid  models,  which  is  a  useful  tool  for  capturing  the  design 
intent. 

Chung  et  al.  (1990)  from  SDRC  developed  two  prototype  feature-based  solid 
modeling  systems  for  mechanical  design.  In  their  research,  they  defined  feature  as  object, 
which  may  contain  private  data  such  as  geometry,  constraints,  inherited  properties  and 
instances  of  other  features,  etc.  The  object  may  have  methods  such  as  geometry  creation 
and  modification,  manufacturing,  analysis  and  documentation.  Based  on  this  definition, 
they  developed  one  prototype  system  for  investment  casting  and  another  for  creating 
sheet  metal  product  model.  Their  results  show  that  geometry  representation  based  on 
features  allows  easy  integration  with  other  applications  such  as  analysis  and  cost 
estimation. 

Hoffmann  and  Juan  (1993)  developed  an  editable,  high-level,  generative 
representation  for  feature-based  solid  modeling,  which  is  called  Erep.  The  developed 
Erep  can  archive  geometric  designs  in  a  form  that  is  both  editable  and  translatable  to  any 
solid  modeling  system.  This  means  the  representation  is  a  high  level  abstraction  of 
design,  which  is  independent  of  any  solid  modeler.  The  geometry  corresponding  to  the 


25 

Erep  is  just  an  instance  of  the  generic  design  and  is  stored  as  a  particular  boundary 
representation  (B-rep),  which  is  created  by  the  solid  modeler  such  as  ACIS.  Therefore,  a 
bridge  between  Erep  and  solid  modeler  is  necessary  for  implementing  feature-based 
modeling  system  based  on  Erep.  For  this  purpose,  Chen  and  Hoffmann  (1995a) 
developed  a  design  compiler  to  translate  Erep  model  to  solid  modeler  ACIS  to  create  B- 
rep  of  the  corresponding  design.  The  compiler  can  perform  tasks  including  constructing 
and  matching  persistent  names,  constraint  solving  by  using  constraint  solver  and  feature 
attaching. 

Duan  et  al  (1993)  reported  that  they  developed  a  feature  based  solid  modeling  tool 
for  feature-based  design  and  manufacture  (FSMT).  The  feature  definition  in  the  system 
includes  geometry  description,  the  attributes  of  a  feature  and  the  relationship  between 
features  and  the  mapping  from  the  definition  to  applications.  The  2D  shape  in  the 
geometry  description  can  be  modified  and  updated  using  variational  geometry  and  thus  is 
dimension-driven.  Therefore,  a  3D  solid  generated  by  sweeping  the  2D  shape  is  also 
dimension-driven.  In  their  system,  mechanical  parts  are  classified  into  five  categories, 
which  are  shaft,  disk,  curved  surface,  box  and  bracket.  A  designer  can  set  up  one  or  more 
feature  libraries  for  each  category.  With  the  help  of  feature  libraries,  designer  can 
efficiently  construct  needed  parts. 

Similar  to  the  work  of  Duan  et  al  (1993),  Peters  et  al  (1994)  presented  a  feature- 
based  modeling  system  using  object-oriented  environment  called  ADAM.  In  this 
research,  a  feature  class  hierarchy  is  defined  using  object-oriented  technique.  They 
suggested  data  base  support  for  feature-based  modeling  and  provided  environment  for 
end-user  (  mechanical  designer)  to  define  and  add  new  feature  type  into  feature  class 
library.  The  authors  claimed  that  the  feature  editing  can  be  realized  by  propagating  the 
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changes  from  one  feature  to  another  and  "propagation"  is  an  action  or  a  collection  of 
actions  that  are  automatically  carried  out  by  the  database  system.  However,  they  only 
used  a  special  example  to  show  how  to  propagate  feature  changes.  The  propagating 
algorithm  can  not  be  found  in  this  paper. 

In  feature-based  modeling,  a  complex  solid  is  usually  constructed  by  combining 
several  features  and  the  operation  that  combine  the  features  is  usually  called  feature 
attachment  (Chen  and  Hoffmann,  1995b).  Since  feature  can  be  seemed  as  an  extension  of 
primitive  in  CSG  operations,  feature  attachment  is  conceptually  very  similar  to  Boolean 
operation.  The  difference  is  that  feature  construction  as  carried  out  in  commercial  CAD 
systems  is  strictly  sequential.  Each  feature  attachment  operation  adds  feature  whose 
shape  and  placement  inseparably  depend  on  the  existing  geometry.  Chen  and  Hoffmann 
(1995b)  thoroughly  analyzed  a  commercial  CAD  system  (Pro/Engineer)  and  proposed  a 
procedural  semantics  for  attaching  protrusion  and  cut  feature  using  regularized  Boolean 
operations.  They  implemented  the  feature  attachment  semantics  based  on  Erep 
(Hoffmann  and  Juan,  1993)  and  the  design  compiler  (Chen  and  Hoffmann,  1995a) 
mentioned  before. 

A  feature-based  modeling  system  usually  needs  standard  feature  library,  which 
includes  features  such  as  holes,  ribs  and  slots.  Although  a  large  number  of  standard 
features  can  be  stored  in  a  system,  it  is  not  possible  to  cover  every  possible  application 
need.  Peters  et  a!  (1994)  suggested  that  end-user  should  be  able  to  define  new  feature 
types  and  add  them  to  a  provided  database,  but  they  didn't  described  how  to  do  this. 
Hoffmann  and  Joan-Arinyo  (1998)  proposed  a  procedural  mechanism  for  generating  and 
deploying  user-defined  features  (UDFs)  in  a  feature-based  modeling  system  based  on 
Erep.  In  their  proposed  mechanism,  a  UDF  can  be  defined  using  standard  features  that  the 
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CAD  system  already  provides  or  using  previously  defined  UDFs.  The  end-user  must  also 
be  able  to  place  and  orient  the  defined  UDF  in  the  part  being  designed.  They  suggested 
that  the  existing  geometric  entities  might  be  used  as  reference  when  placing  and  orienting 
the  UDF,  They  proposed  using  constraints  such  as  mating  constraints,  however,  details 
about  this  process  are  not  included  in  this  paper  except  they  mentioned  using  mating 
constraints. 

Geometry  Updating  in  History-based  Parametric  Modeling 

History-based  parametric  modeling  technique  is  being  used  in  all  modem 
commercial  CAD  systems,  which  have  the  capability  of  dimension-driven  editing.  The 
geometry  updating  (regeneration)  algorithm  is  very  important  in  this  technique.  Since 
history-based  parametric  modeling  technique  was  firstly  developed  in  industry,  details  of 
the  updating  algorithm  in  commercial  systems  are  not  published.  However,  it  is  well 
known  that  geometry  is  updated  by  "rolling  back"  to  a  previous  state  and  reconstructing 
the  geometry  in  these  systems. 

One  of  the  problems  related  to  geometry  updating  in  history-based  parametric 
modeling  is  persistent  naming  (or  generic  naming).  Persistent  naming  is  a  technique  used 
to  establish  correspondence  between  topological  entities  of  B-reps  before  and  after  the 
model  has  been  updated.  For  an  example,  one  operation  is  to  position  a  B-rep  A  with 
respect  to  B-rep  B  by  enforcing  one  planar  face  of  A,  /a,  to  be  coincident  to  one  planar 
face  o{B,fB.  After  updating  A  and  B,  they  are  now  new  B-reps  and  the  system  must  be 
able  to  identify  the  two  faces  corresponding  to  /a  and  fs  to  reapply  the  positioning 
operation.  The  first  approach  directly  addressing  this  problem  was  proposed  by  Kripac 
(1995,  1997).  In  his  approach,  topological  entities  can  be  identified  by  using  editing 
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history  information  called  FaceldGraph,  topological  adjacency  information  and 
geometric  orientation  information.  An  object-oriented  software  architecture  with  a 
narrow  application  interface  was  also  developed  in  this  research.  Capoyleas  et  al  (1996) 
and  Chen  and  Hoffmann  (1995c)  also  studied  this  problem.  Their  method  starts  with  the 
names  of  the  primitive  faces  in  the  parametric  definition  and  constructs  the  generic  names 
for  topological  entities  using  extended  topological  adjacency  information,  called  context 
graph,  and  geometric  orientation  information.  In  Chen  and  Hoffmann  (1995c),  feature 
editing  operations,  which  includes  "roll  back"  updating  procedure,  are  also  introduced. 
Both  of  above  naming  algorithms  will  not  always  properly  handle  collapsing  maps, 
which  is  the  case  that  when  a  n-dimensional  boundary  entity  maps  to  a  m-dimensional 
one  (m<n),  such  as  when  a  face  maps  to  an  edge.  The  failure  of  the  algorithms  will  lead 
to  unexpected  resuh  in  the  updating  process.  Raghothama  and  Shapiro  (1998)  proposed  a 
formal  definition  of  Boundary  Representation  Deformation  (BR-deformation)  to  avoid 
unexpected  updating  results  based  on  the  assumption  of  continuity:  small  changes  in  solid 
parameter  values  should  result  in  small  changes  in  solid's  B-rep.  They  claimed  that  the 
necessary  conditions  that  must  be  satisfied  by  any  BR-deforming  mapping  can  be  used  to 
identify  invalid  update.  Therefore  the  model  change  being  a  BR-deformation  is  the 
necessary  condition  for  the  above  naming  algorithms  to  provide  a  correct  solution. 

Raghothama  and  Shapiro  (1999)  proposed  another  idea  for  geometry  updating  in 
history-based  parametric  modeling.  This  idea  is  to  modify  parametric  definition  and  B- 
rep  simultaneously  through  an  equivalent  semantics  for  every  edit  in  both  representations. 
Since  this  is  a  direct  modification  to  B-rep,  no  new  B-rep  is  created.  Therefore  persistent 
naming  is  not  necessary  as  long  as  we  know  that  the  modified  parametric  definition  is 
consistent  with  modified  B-rep.  The  authors  emphasized  the  consistency  test  and 
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presented  two  algorithms  to  determine  if  a  modified  parametric  definition  is  consistent 
with  the  corresponding  modified  B-rep.  However,  they  did  not  discuss  how  to  perform 
this  type  of  updating. 

Rigid  Body  Positioning  by  Spatial  Relations 

In  our  research,  a  soHd  is  constructed  by  assembling  primitives  together  and  then 
using  Boolean  attaching  operations  to  create  the  resultant  solid.  In  this  process,  we  need 
to  position  one  primitive  or  intermediate  geometry  with  respect  to  another  by  specifying 
spatial  relations,  such  as  "plane  to  plane  angle"  and  "vertex  to  plane  distance",  between 
geometric  entities  of  two  objects  involved  in  the  attaching  operations.  This  problem  can 
be  stated  as  Rigid  Body  Positioning  problem.  After  extensively  searching  relevant 
literature,  research  on  this  topic  was  found  to  be  surprisingly  few  and  only  four  papers 
were  found  closely  related  to  this  problem. 

Ambler  and  Popplestone  (1975)  studied  the  problem  in  the  context  of  robot 
manipulation  of  3D  assemblies.  The  allowed  types  of  primitives  to  be  assembled  are 
block,  cylinder  and  sphere  and  the  entities  allowed  are  face,  edge  and  vertex.  The 
constraints  defined  are  only  against  and  fits,  against  holds  between  two  plane  faces,  or 
between  a  plane  face  and  a  cylindrical  shaft,  fits  holds  between  a  cylindrical  shaft  and  a 
hole.  Matrix  equations  are  derived  from  these  constraints  and  can  be  solved 
simuhaneously  to  obtain  the  transformation  matrices.  Their  method  uses  an  algebraic 
approach  to  reduce  the  computation  time,  but  it  requires  a  large  amount  of  computation  to 
rewrite  related  equations  in  a  solvable  form.  The  another  problem  of  this  method  is  that  it 
does  not  guarantee  the  solution  each  time,  especially  when  the  equations  can  not  be 
rewritten  in  solvable  form. 
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Another  similar  method  was  proposed  by  (Lee,  K.  and  Gossard,  D.  C,  1985,  Lee, 
K.  and  Andrews,  G.,  1985)  for  assembly  modeling.  They  firstly  developed  a  hierarchical 
data  structure  for  representing  assemblies  and  then  they  used  spatial  relations  to  position 
and  orient  the  parts  involved  in  the  assembly.  The  spatial  relations  or  constraints  they 
used  are  also  against  and  fit.  After  creating  equations,  Newton-Raphson  iteration  method 
is  used  to  solve  these  equations  so  that  the  needed  translation  and  rotation  matrices  can  be 
obtained. 

Both  of  the  above  methods  need  to  solve  non-linear  simultaneous  equations.  An 
alternate  approach  to  solve  this  problem  is  to  sequentially  impose  each  constraint  such 
that  previously  imposed  constraints  are  not  violated.  Kramer  (1991,  1992)  proposed  a 
Degree  of  Freedom  Analysis  method  to  solve  geometric  constraint  system  in  kinematic 
analysis  of  mechanisms.  The  purpose  of  this  method  is  to  sequentially  position  a  rigid 
body  with  respect  to  another  such  that  the  given  geometric  constraints  are  satisfied.  In 
this  method,  Kramer  developed  semantics  for  defining  constraint  type,  recording  status  or 
rigid  body's  current  degrees  of  freedom  and  describing  the  necessary  actions  for 
transforming  the  body  to  satisfy  the  constraints.  Based  on  the  semantics,  a  Plan  Fragment 
Table  (PFT)  is  created.  Each  entry  in  PFT  has  information  about  the  type  of  the 
constraint,  initial  status  of  DOF,  new  status  of  DOF  and  sequential  actions  needed  for 
transforming  the  rigid  body.  If  a  geometric  constraint  is  given,  the  software  based  on  this 
method  will  analyze  the  current  status  of  DOF  and  the  new  status  of  DOF,  and  then  find 
the  corresponding  entry  in  the  PFT.  Once  the  entry  is  found,  sequential  actions  specified 
by  PFT  will  be  taken  to  do  the  transformation.  The  advantage  of  this  method  is  that  it 
does  not  require  solving  complex  nonlinear  equations.  However,  since  this  method  is 
developed  for  constraint  solving  in  mechanisms,  the  constraint  types  considered  in  this 
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method  are  not  enough  for  solid  modeling.  For  example,  only  "point  to  point  coincident", 
"point  to  axis  coincident"  and  "point  to  plane  coincident"  are  considered  as  translation 
constraints  in  this  method,  but  distance  constraints,  such  as  "point  to  plane  distance", 
which  the  method  does  not  consider,  is  a  very  usual  constraint  type  used  in  dimension- 
driven  solid  modeling.  Another  disadvantage  of  this  method  is  that  a  rather  large  number 
of  entries  would  be  necessary  in  the  Plan  Fragment  Table  to  handle  all  the  constraints 
required  to  fully  constrain  typical  geometric  models.  Moreover,  the  PFT  does  not  have  a 
solution  for  some  constraints  that  often  occur  in  solid  modeling.  For  example,  when  two 
rotational  DOFs  are  left  and  an  "axis  to  axis  angle"  constraint  is  to  be  imposed,  PFT 
suggests  using  numerical  method  to  solve  the  constraint. 


CHAPTER  3 
REPRESENTATION  OF  DESIGN  INTENT 

Overview 

Design  intent  representation  is  the  basis  of  dimension-driven  solid  modeling.  As 
we  reviewed  before,  modem  commercial  CAD  systems  use  history-based  parametric 
modeling  technique  to  do  dimension-driven  editing  of  solid  model,  in  which  the 
geometry  of  a  primitive  is  parametrized  as  a  high-level  data  structure  called  feature  and 
the  design  intent  is  represented  by  a  sequence  of  feature  operations.  In  our  research,  an 
extended  CSG  tree  is  used  to  represent  design  intent.  In  this  chapter,  we  define  extended 
CSG  tree  in  first  section.  In  rest  of  this  chapter,  we  explain  the  concept  of  attach  operator, 
complete  constraint  set  and  primitive  creation  function,  which  are  information  stored  in 
extended  CSG  tree. 

Extended  CSG  Tree 

A  complex  solid  model  can  be  constructed  fi"om  many  simple  geometric 
primitives  by  a  series  of  attach  operations.  These  attach  operations  could  be  various 
regularized  Boolean  operations.  Figure  3-1  is  a  3D  example  showing  the  generation 
process  and  the  CSG  tree  representing  this  process.  Construction  process  for  a  complex 
2D  geometry  can  also  be  represented  by  a  CSG  tree,  in  which  all  primitives  and 
intermediate  geometries  are  2D. 

In  this  example,  there  are  five  3D  primitives  PI,  P2,  P3,  P4  and  P5.  Intermediate 
geometry  Gl  is  obtained  by  subtracting  P2  from  PI.  Intermediate  geometry  G2  is 
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obtained  by  adding  P3  to  Gl.  Intermediate  geometry  G3  is  obtained  by  adding  P5  to  P4. 
The  final  geometry  G4  is  generated  by  subtracting  G3  fi-om  G2.  Each  internal  node 
represents  a  regularized  Boolean  operation.  Before  performing  a  Boolean  operation,  it  is 
possible  that  one  of  primitives  or  intermediate  geometries  involved  in  this  Boolean 
operation  needs  rigid  body  motions  to  be  correctly  positioned  with  respect  to  another.  For 
simplicity,  these  rigid  body  motions  are  omitted  in  this  CSG  tree. 


PI         P2         P3        P4  P5 
Figure  3-1 :  CSG  tree  recording  model  generation  history 

The  traditional  CSG  tree  shown  in  figure  3-1  is  a  good  representation  of  the  final 
solid  model.  The  main  information  stored  in  an  internal  node  of  traditional  CSG  tree 
includes  the  operation  applied  on  its  children,  which  could  be  a  certain  type  of  Boolean 
operation  or  a  rigid  motion  (translation  or  rotation).  However,  it  does  not  uniquely 
represent  the  construction  history  or  design  intent.  The  rigid  body  motions  in  the 
traditional  CSG  tree  only  represent  the  transformations  in  the  current  model  generation 
process.  When  dimensions  are  changed,  the  transformations  needed  could  be  different. 
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However,  the  traditional  CSG  tree  does  not  record  the  information  regarding  how  to 
change  these  rigid  body  motions  in  case  of  dimension  changes. 

Traditional  CSG  tree  does  not  record  the  unique  sequence  of  operations  for 
generating  the  model.  For  an  example,  two  feasible  sequences  corresponding  to  CSG  tree 
in  figure  3-1  are: 

{Cpi ,  Cp2 ,  P1-*P2  ,  CP3 ,  Glu*P3,  Cp4 ,  CP5,  P4^'P5,  G2-*G3}  (3.1) 

and 

{  Cp5 ,  Cpi ,  Cp2 ,  Cp3 ,  Cp4 ,  P4u  P5,  P1-'P2,  Glu  P3,  G2-*G3}  (3.2) 
where  Ci  is  a  primitive  creation  operation  for  primitive  /. 

For  the  updating  algorithm  introduced  in  chapter  5,  however,  it  is  not  necessary 
that  the  original  sequence  of  operations  for  generating  the  model  is  recorded.  As  long  as  a 
sequence  of  operations  can  generate  the  solid  model  corresponding  to  the  CSG  tree,  it  is  a 
feasible  sequence  and  can  be  used  as  the  part  of  model  generation  history  for  geometric 
updating  (regeneration).  It  is  easy  to  see  that  traditional  CSG  tree  records  a  feasible 
sequence,  which  can  be  obtained  by  traversing  the  tree  in  post-order.  Sequence  (3.1)  is 
such  a  sequence. 

To  record  the  model  generation  history  or  design  intent,  we  extend  the  CSG  tree 
by  adding  additional  information,  which  includes  Complete  Constraint  Set  and  primitive 
creation  function,  to  the  internal  nodes  and  leaf  nodes  respectively.  If  there  is  a  node 
representing  rigid  body  motion  in  CSG  tree,  this  node  will  be  deleted  in  the 
corresponding  extended  CSG  tree.  Moreover,  in  the  extended  CSG  tree,  the  Boolean 
operation  is  replaced  by  attach  operator,  which  is  stored  in  internal  nodes. 
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Extended  CSG  tree  can  be  used  to  represent  the  design  intent  of  a  2D  or  3D 
geometry.  If  the  geometry  is  a  3D  solid,  a  2D  profile,  represented  by  an  extended  CSG 
tree,  will  be  stored  in  each  leaf  node  and  used  as  one  of  the  inputs  for  the  primitive 
creation  function. 

Attach  Operator  . 

Attach  Operator  could  be  a  class  or  a  function  stored  in  internal  nodes  that 
attaches  the  geometry  of  right  child  to  geometry  of  left  child,  which  could  be  primitive 
geometry  or  intermediate  geometry.  In  general  case,  the  attach  process  is  performed  by 
doing  regularized  Boolean  union,  difference  or  intersection  based  on  the  type  of  attach 
operator.  The  inputs  of  each  attach  operator  are  the  B-reps  of  the  two  children.  The  result 
is  the  new  geometry  of  the  attach  operation. 

Complete  Constraint  Set 

Before  explaining  the  Complete  Constraint  Set,  we  need  to  make  an  assumption 
for  the  extended  CSG  tree  creation.  We  will  assume  that  the  left  child  of  an  internal  node 
of  the  tree  represents  the  geometry,  which  does  not  need  any  repositioning  and  is  defined 
as  base  geometry  involved  in  the  Boolean  operation.  On  the  other  hand,  the  right  child 
represents  the  geometry,  which  needs  to  be  repositioned  with  respect  to  base  geometry 
and  is  defined  as  the  dependent  geometry  in  the  Boolean  operation.  The  Complete 
Constraint  Set  is  a  set  of  position  constraints  that  completely  fixes  the  six  degrees  of 
freedom  of  the  dependent  geometry  with  respect  to  the  base  geometry  and  correctly 
positions  the  dependent  geometry  based  on  the  design  intent. 
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Each  position  constraint  is  a  spatial  relation  specified  between  two  entities,  which 
belong  to  the  base  geometry  and  the  dependent  geometry  respectively.  For  example,  we 
can  specify  a  distance  between  a  vertex  on  the  dependent  geometry  and  a  face  on  the  base 
geometry.  The  advantage  of  specifying  a  Complete  Constraint  Set  for  the  tree  node 
instead  of  a  rigid  body  motion  is  that  the  transformation  required  to  reposition  the 
dependent  geometry  can  be  automatically  computed  such  that  the  position  constraints  are 
satisfied,  when  any  dimensions  are  changed.  The  method  for  imposition  of  this  Complete 
Constraint  Set  will  be  discussed  chapter  4. 


Primitive  Creation  Function 

A  primitive  geometry  is  initially  created  by  a  primitive  creation  function,  which  is 
stored  in  the  corresponding  leaf  node.  For  a  2D  primitive,  the  fiinction  takes  as  argument 
the  parameters  or  dimensions  of  the  geometry: 

s,  =G,(x„,...,x^,...,x,^^)  (3.3) 

where  Si  is  the  generated  2D  profile,  which  is  a  planar  face,  G,  is  the  primitive  creation 
function  ,  Xik  is  the  kth  dimension  of  this  primitive  i  and  ngi  is  the  number  of  dimensions 
or  independent  variables  of  primitive  i.  Many  primitive  creation  functions  can  be  pre- 
defined for  2D  models.  When  a  new  primitive  is  created  its  geometry  will  depend  on  its 
dimensions  and  the  particular  primitive  creation  function  assigned  to  it.  For  example,  the 
primitive  creation  function  for  a  square  takes  the  side  length  as  input  and  creates  a  square 
face. 
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For  3D  primitive,  the  primitive  creation  function  creates  geometry  by  applying  a 
sweep  operation  using  a  2D  profile  stored  in  the  leaf  node. 

=G,{s^^,x,,,...,x,,,...,x,„^^)  (3.4) 

where  Si  is  resuhed  solid  model,  S2d  is  the  2D  profile  and  other  inputs  are 
parameters  or  dimensions  needed  for  building  this  3D  geometry.  In  a  typical 
implementation  one  may  provide  functions  for  extruding,  revolving,  blending  and 
sweeping  a  2D  profile.  A  3D  block,  for  example,  can  be  created  by  using  an  extrusion 
function.  The  input  2D  profile  can  be  a  square  or  rectangular  face.  Other  inputs  the 
extrusion  fimction  takes  include  the  depth  of  extrusion  and  the  parameters  to  define  the 
extrusion  direction. 


CHAPTER  4 
SEQUENTIAL  CONSTRAINT  IMPOSITION 

Overview 

In  this  chapter,  a  method  for  sequentially  imposing  constraints  in  Complete 
Constraint  Set  is  presented  that  is  particularly  suited  for  imposing  position  constraints  to 
locate  geometric  primitives  or  parts  with  respect  to  each  other.  It  is  assumed  that  the 
geometric  primitives  are  solid  models  or  2D  profiles  that  have  been  created  using  pre- 
defined procedures.  These  primitive  are  "parameterized"  such  that  if  key  parameters 
(dimensions)  are  varied  the  solid  model  or  profile  can  be  either  reconstructed 
automatically  by  following  the  pre-defined  procedure  or  updated  by  using  variational 
geometry  technique  or  the  updating  algorithm  proposed  in  chapter  5  in  this  dissertation. 
The  sequential  constraint  imposition  method  described  here  can  be  used  to  position  these 
parameterized  primitives  with  respect  to  each  other.  In  the  next  section,  the  sequential 
imposition  of  position  constraints  is  introduced  and  its  application  to  dimension-driven 
solid  modeling  is  explained.  Thereafter,  the  position  constraints  are  classified  into 
angular  and  distance  constraints  and  the  method  for  imposing  these  two  types  of 
constraints  sequentially  are  discussed  in  section  3  and  4  respectively.  Application  of  this 
constraint  imposition  technique  to  solid  modeling  is  illustrated  using  an  example  in 
section  5  followed  by  discussion  in  section  6. 
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Dimension-driven  Editing  Using  Position  Constraints 

Geometric  constraints  can  be  sequentially  imposed  to  position  one  geometry  with 
respect  to  another  if  each  constraint  is  imposed  such  that  previously  imposed  constraints 
are  not  violated.  While  this  approach  could  be  used  to  position  vertices,  edges  and  faces 
relative  to  each  other,  it  is  difficult  to  do  so  in  practice  due  to  the  large  number  of 
possible  types  of  constraints,  curves,  surfaces  and  combinations  of  these  that  need  to  be 
considered.  Instead,  one  could  use  the  sequential  approach  to  position  geometric 
primitives  (or  simple  geometry)  with  respect  to  one  another  so  that  they  can  be  combined 
using  Boolean  operations  or  other  attach  operators  to  define  the  final  geometry.  The 
process  of  positioning  primitives  and  combining  them  using  attach  operations  is 
illustrated  in  figure  4-1  using  a  simple  2D  example  and  compared  with  variational 
geometry. 

The  2D  geometry  shown  in  figure  4-1  can  be  created  using  Constructive  Solid 
Geometry  (CSG),  where  simple  primitives  (A,  B  and  C)  are  combined  using  regularized 
Boolean  operations  to  create  the  final  geometry  E.  The  figure  also  shows  the  Boolean 
tree  structure  where  the  leaves  are  the  primitives  and  the  nodes  represent  regularized 
Boolean  operations  and  rigid  body  transformations.  If  the  primitives  are  generated  by 
pre-defined  procedures,  it  can  be  easily  reconstructed  for  any  value  of  the  dimensions  that 
defines  the  primitive.  If  primitives  are  positioned  relative  to  each  other  by  specifying 
position  constraints,  the  rigid  body  transformation  can  also  be  automatically  recomputed 
when  any  dimensions  are  modified.  After  positioning  the  primitives  the  Boolean 
operations  can  be  applied  to  update  the  final  geometry  in  response  to  dimension  changes. 
In  contrast,  for  the  variational  geometry  approach,  one  would  specify  geometric 
constraints  directly  between  vertices  and  edges  of  the  geometry.  In  the  final  geometry  (E) 
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in  figure  4-1  such  constraints  are  indicated  using  symbols  such  as  H  (line  horizontal),  V 
(line  vertical),  1  (perpendicular)  etc.  The  equations  for  these  geometric  and  dimensional 
constraints  are  then  derived  and  solved  simultaneously  to  obtain  the  final  coordinates  of 
the  vertices  and  the  geometry  of  the  edges. 


A  d6 


Figure  4-1:  Variational  geometry  versus  sequential  constraint  imposition 


The  concept  illustrated  in  figure  4-1  can  be  extended  to  3D  solid  models  by  using 
3D  primitives  that  are  constructed  by  sweeping  2D  profiles  along  curves  or  by  revolving 
the  profiles  about  predefined  axes.  Position  constraints  can  be  used  to  position  the 
primitives  with  respect  to  each  other  before  applying  the  Boolean  operations.  The  2D 
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profiles  used  to  construct  the  3D  primitives  can  themselves  be  constructed  by  combining 
2D  primitives  that  are  also  positioned  relative  to  each  other  using  position  constraints.  In 
following  sections,  we  present  a  generalized  method  for  sequentially  imposing  such 
position  constraints  and  discuss  some  of  the  advantages  and  limitations  of  the  approach. 

We  will  use  the  term  object  to  refer  to  both  2D  and  3D  primitives  in  this  chapter. 
To  position  objects  with  respect  to  one  another,  we  assume  that  the  position  constraints 
are  specified  between  points,  lines  or  planes  on  the  two  objects.  The  constraints  are 
imposed  by  rigid  body  transformation  of  one  object  v^th  respect  to  the  other.  We  refer  to 
the  object  that  remains  stationary  as  the  base  (or  parent)  and  the  transformed  object  as  the 
dependent  (or  child).  To  fully  constrain  a  3D  dependent  object  with  respect  to  its  base, 
enough  constraints  must  be  specified  to  remove  all  six  degrees  of  freedom  of  the 
dependent  object.  We  can  classify  position  constraints  into  angular  constraints,  which 
specify  angle  between  entities  and  distance  constraints,  which  specify  distance  (or 
coincidence)  between  entities.  Each  position  constraint  removes  one  or  more  degrees  of 
fi"eedom  of  the  object  to  be  positioned.  A  set  of  position  constraints  C  that  uniquely 
positions  the  dependent  object  and  fully  constrains  it  will  be  referred  to  here  as  a 
complete  constraint  set.  It  can  be  represented  using  a  list  of  constraint  functions, 

C  =  {  Cal,  •  • .,  Can,  Ctl,  ■  ■  ■,  Ctm  }  =  {  Ca,  Q  }  (41) 

where  Cai  are  the  angular  constraint  functions  and  Cdi  are  the  distance  constraint 
functions.  Ca  is  a  complete  angular  constraint  set  if  it  fully  constrains  all  rotational 
degrees  of  freedom  of  the  dependent  object.  Similarly,  Ct  is  a  complete  translation 
constraint  set  if  it  eliminates  all  the  translation  degrees  of  freedom.  Each  type  of 
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constraint  Cj  (angular  or  distance)  can  be  implemented  as  a  function  that  computes  the 
transformation  matrix  [Ti]  that  imposes  the  required  constraint. 

[Ti]  =  Ci(e;,e2,d,Cp)  (4.2) 
In  the  above  equation,  ei  and  e:  are  generic  names  of  geometric  entities  of  the  dependent 
and  base  objects  respectively  between  which  the  constraint  is  specified.  Any  dimensions 
or  parameters  needed  to  specify  the  constraint  is  denoted  by  the  array  d,  which  could  be 
angles,  distances  or  just  flags  used  to  distinguish  between  muhiple  solutions.  Cp  is  the 
array  of  previously  imposed  constraints.  The  transformation  matrix  [Ti]  is  computed  by 
the  function  Ci  such  that  it  transforms  the  object  to  satisfy  the  constraint  without  violating 
any  constraints  in  Cp. 

It  is  possible  to  fully  constrain  one  object  with  respect  to  another  without  using  a 
complete  angular  constraint  set.  For  example,  one  could  fully  constrain  an  object  using 
only  distance  constraints.  However,  in  this  research,  we  will  consider  only  the  case  where 
Ca  is  a  complete  angular  constraint  set  so  that  all  rotational  degrees  of  freedom  are 
removed  using  angular  constraints.  This  is  not  a  severe  limitation  for  application  to  solid 
modeling  because  it  is  often  more  intuitive  to  use  angle  constraints  to  orient  an  object. 

For  a  given  complete  set  of  constraints,  it  is  possible  to  impose  the  constraints  in 
many  different  sequences.  However,  it  is  not  practical  to  develop  an  algorithm  that  allows 
any  arbitrary  order  of  constraint  imposition  since  for  every  type  of  constraint  the  method 
for  deriving  the  transformation  matrix  depends  on  what  constraints  were  previously 
imposed.  Instead,  we  select  a  sequence  that  simplifies  the  procedure  for  deriving  the 
transformation  matrix  and  ensures  that  constraints  can  be  imposed  without  violating 
previous  constraints.  Firstly,  we  make  use  of  the  fact  that  translation  does  not  violate 
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previously  imposed  angular  constraints  while  the  rotation  can  violate  distance  constraints. 
Therefore,  we  impose  the  angular  constraints  first  to  remove  the  rotational  degrees  of 
freedom.  Up  to  three  constraints  may  be  required  to  fix  all  the  angular  degrees  of  freedom 
and  each  of  them  must  be  imposed  without  violating  any  previous  angular  constraints. 
Similarly,  each  distance  constraint  must  be  imposed  without  violating  any  previously 
imposed  distance  constraints.  The  following  sections  describe  methods  for  deriving  the 
necessary  transformation  matrices. 

Sequential  Imposition  of  Angular  Constraints 

Entities  that  could  be  used  to  define  angular  constraints  are  lines  (edges  or  axes) 
or  planes  (faces  or  reference  planes).  Since  a  vector  can  be  used  to  represent  both  the 
orientation  of  a  line  and  the  normal  of  a  plane,  there  are  in  fact  only  two  types  of  angular 
constraints.  The  first  type  is  a  constraint  specifying  the  angle  between  two  vectors  that 
could  be  used  to  specify  angle  between  lines,  line  and  plane  or  plane  and  plane.  This 
constraint,  which  we  will  refer  to  as  the  "vector-to-vector  angle"  (VVA)  constraint,  fixes 
one  rotational  degree  of  freedom.  The  second  constraint  type  is  the  "parallel  vectors" 
(PV)  constraint  that  specifies  that  two  vectors  are  parallel  and  thereby  fixes  two  rotational 
degrees  of  freedom.  This  constraint  can  be  used  to  make  two  planes  parallel  or  two 
edges/axes  parallel  etc.  Using  these  two  constraints,  one  can  define  two  types  of 
Complete  Angular  Constraint  Sets  (CACS).  One  is  a  set  consisting  of  a  PV  constraint 
(fixing  2  d.o.f )  and  a  VVA  constraint  (fixing  1  d.o.f),  which  we  will  refer  to  as  CACS-I. 
The  other  set  consists  of  three  WA  constraints  (fixing  one  d.o.f  each)  that  we  will  refer 
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to  as  CACS-II.  Sequential  imposition  of  these  two  types  of  complete  angular  constraint 
sets  is  described  in  the  next  two  sub-sections. 


CACS-I:  Parallel  vectors  and  vector-to-vector  angle  constraints 

If  the  angular  constraint  set  consists  of  a  PV  and  a  WA  constraint,  they  can  be 
imposed  in  either  sequence.  However,  since  PV  fixes  two  degrees  of  fi-eedom  it  is 
advantageous  to  impose  it  first  as  shown  in  the  figure  4-2.  Let  b,  and  bj  be  the  unit 
vectors  on  the  base  object  and  d,  and  d^  be  the  unit  vectors  on  the  dependent  object 
between  which  the  angles  are  specified.  The  WA  constraint  specifies  that  the  angle 
between  b,  and  d,  is  61,  0<9i<90°,  and  the  PV  constraint  specifies  that  b^  is  parallel  to 

dj .  The  range  of  9i  is  not  a  significant  limitation  because  one  can  always  use  opposite 
direction  of  one  of  vectors  if  9i>90°  is  needed.  The  angle  between  d,  and  dj  (an)  and 
the  angle  between  b,  and  bj  (P12)  are  fixed  since  these  are  vectors  attached  to  the  same 


object. 


Cone-II 


Figure  4-2:  Sequentially  imposing  CACS-I,  Cr-  {PV,  WA) 
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In  figure  4-2,  the  original  orientation  of  the  two  vectors  d,  and  is  shown  using 
dotted  lines.  The  constraints  can  be  imposed  by  two  rigid  body  rotations  of  the  dependent 
object  by  angles     and  ^2  as  shown  in  the  figure.  The  first  rotation  is  about  the  vector  n 

normal  to  both  dj  and  bj  which  can  be  computed  as  n  =  d^  x  bj .  The  rotation  angle  ^\ 
is  in  the  range  from  0  to  180°  due  to  the  property  of  cross-product.  Therefore,  the  angle 
4)1  can  be  computed  as  following: 

^,=cos-\d,\)  (4.3) 
After  rotating  the  dependent  object  by  the  angle  ^i,  d[  and  d'^  are  the  transformed 

positions  of  dj  and  d^  respectively. 

To  impose  the  WA  constraint  without  violating  the  PV  constraint,  the  dependent 
object  must  be  rotated  about  the  d'^  (or  bj )  vector.  The  angle  of  rotation  is  computed 
such  that  after  the  rotation  d\  is  positioned  at  the  intersection  of  cone-I  and  cone-II  as 
shown  in  the  figure.  Cone-I  has  bj  as  its  axis  and  the  angle  between  the  cone  surface  and 
the  axis  (cone-angle)  is  ai2.  Cone-II  has  b,  as  its  axis  and  a  cone-angle  of  9i.  The 
intersection  of  cone-I  and  cone-II  represents  a  direction  that  makes  an  angle  9i  with  b, 

and  an  angle  an  with  b^  (and  d^ ).  Notice  that  depending  on  the  relative  position  and 
angles  of  the  cones,  these  cones  may  intersect  along  two  lines,  touch  each  other  along  one 
line  or  not  intersect  at  all.  Therefore  the  constraint  set  CACS-I  could  have  two  solutions, 
one  solution  or  no  solution.  From  figure  4-2,  it  is  clear  that  the  cones  will  intersect  and 
the  constraint  set  will  have  solutions  if  the  following  condition  is  satisfied. 

|Pi2-ai2|  <Gi<pi2  +  ai2  (4.4) 
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where  P12  is  the  angle  between  b,  and  b^,  0<  P12  <180°,  an  is  the  angle  between  d,  and 
d,,  0<ai2<180°. 

The  angle  of  the  second  rotation  ^2  about     that  would  position  d\  at  an  angle  61 
from  b,  can  be  computed  as 

(t),=/(d;,b„b„0,)  (4.5) 

Equation  (4.5)  can  be  abstracted  as  following  problem:  compute  the  angle  (j)  by 
which  to  rotate  a  unit  vector  a  about  another  unit  vector  b  so  that  it  makes  an  angle  0 
with  a  given  unit  vector  c  .  The  final  position  of  the  vector  a  after  the  rotation  is  shown 
as  a  . 


(a)  (b) 

Figure  4-3 :  Rotation  angle  computation 

If  a  overlaps  b  (a  =  0),  we  will  check  whether  a  •  c  =  cos^ .  If  a  •  c  =  cos^ ,  the  function 
(4.5)  returns  an  angle  0,  otherwise,  there  is  no  solution  for  (j).  If  c  overlaps  b,  we  will 
check  whether  6  =  a.  If  0  =  a,  the  function  (4.5)  returns  an  angle  0,  otherwise,  there  is  no 
solution  for  (J).  As  shown  in  Figure  4-3,  in  general  case  other  than  above  two  extreme 
cases,  the  final  position  of  a'  can  be  computed  as: 

a'  =  OP  +  PB  =  cosab  +  rcos(t)n  +  rsin  (4.6) 
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where,  r  =  sina  is  the  length  of  PA  and  PB  and 


PA  a-cosab 


PA 


a  -  cosab 


.  OPxPA  £  . 
m  =  =z — =r  =  bxa 

OPxPA 


(4.7) 


The  angle  of  rotation  (j)  must  be  computed  such  that  the  angle  between  a'  and  c  is  9  so 


that  a'  •  c  =  cosG .  Substituting  (4.7)  into  (4.6)  we  get, 
pcos^  +  qsm^  =  k 


(4.8) 


where. 


p  =  nc- 


(a  -  cosab)  c 


a-cosab 


^  =  mc  =  (bxa)c 
cos^-bccosa 


(4.9) 


sm  a 


Dividing  two  sides  of  (4.8)  by  ^Jp^        and  let 


cosy 


sm  / 


(4.10) 


The  value  of  y  is  uniquely  determined  by  (4. 10).  We  have  equivalent  equation  of  (4.8)  a; 

k 


cos(y  -  <p) 


(4.11) 


Therefore 


(p  =  /±cos~ 


(4.12) 
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Since  the  values  of  p  and  q  are  not  0  in  general  case,  ^  has  two  solutions  from  (4.12) 


when  >     ■  When  k  =  —\Jp^        ,  </> has  two  repeated  solutions  /  ± ;r which  are 

actually  same.  When  k  =  ■\lp^  +q^  ,  ^ has  one  solution  y . 

After  (|)i  and  ^2  in  figure  4-2  are  obtained,  the  transformation  matrix  for  imposing 

this  constraint  set  can  be  expressed  as 

[TJ  =  [T,(b„(z),)][T,(n,<^,)]  (4.13) 
[T^(n,(t))] ,  the  transformation  matrix  that  rotates  an  object  by  an  angle  ^  about  the 

vector  n  =      riy,  rizY ,  is  defined  as  following  (Zeid,  1991): 


«^(l-cos</J)  +  cos(Z>  /7^«^(l-cos(^)-A7^  sin^  n^/7^(l-cos(^)  +  «^sin^ 
n^ny{\-cos</>)  +  b^sm^  nl(\-cos(p)  +  cos(p  n^/7^(l  -  cos(Z>) -n,  sin  ^ 
n^n^(\ -cos(^)- by  sin  ^   /7_^«,(1  -  cos^) +     sin  <^      «^(l-cos(Z>)  +  cos<^ 

(4.14) 


CACS-II:  Three  vector-to-vector  angle  constraints 

In  the  most  general  case,  three  distinct  vectors  from  each  object  are  used  to  define 
the  three  VVA  constraints.  However,  this  is  not  a  convenient  or  suitable  way  to  orient 
one  object  with  respect  to  the  other  because  a  solution  exists  only  in  special  cases  for  a 
small  range  of  angles.  Furthermore,  one  cannot  intuitively  determine  whether  a  solution 
exists.  For  these  reasons,  we  only  consider  the  case  where  two  constraints  use  the  same 
vector  on  the  dependent  object.  Without  loss  of  generality,  we  can  assume  that  the 

constraints  Cai  and  Ca2  refer  to  the  same  vector  d, .  Cai  specifies  the  angle  9i  between  d, 
and  bj ,  Cai  specifies  the  angle  62  between  d,  and  bj ,  and  Ca3  specifies  the  angle  83 
between  dj  and  bj .  All  61,  62  and  63  are  in  the  range  from  0  to  90°  . 
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Figure  4-4:  CACS-II  with  two  distinct  vectors  on  the  dependent  object 

The  first  constraint,  Cai,  can  be  imposed  by  positioning  d;  on  the  surface  of  cone- 
I  (axis  b, ,  angle  9i).  This  can  be  achieved  by  rotating  the  object  about  a  vector  n  normal 
to  both  di  and  bi  by  an  angle     where  n  =  dj  x  bj  and  (j)i  can  be  computed  using  (4.15). 

cos(<z), +^)  =  d,  b,.  (4.15) 

Ca2  can  be  imposed  by  positioning  d,  on  the  surface  of  cone-II  (axis  bj ,  angle  Qi) 
by  rotating  the  dependent  object  about  by  an  angle  (j)2  as  shown  in  Figure  4-4  (a).  The 
angle  can  be  computed  as  (j).  =  /(d,,b,,b2,92).  The  new  position  of  d, ,  which  satisfies 
both  Cai  and  Ca2,  is  at  the  intersection  of  cone-I  and  cone-II  and  is  denoted  as  dj  in  the 
figure.  A  solution  exists,  if  62  satisfies  the  following  equation: 

|pi2-ei|<02<pi2  +  ei  (4.16) 
where  Pn  is  the  angle  between  b,  and  bj,  0<  Pn  <180°. 

To  avoid  violating  Cai  and  Ca2,  Ca3  can  be  imposed  only  by  rotating  the  object 
about  dj .  The  final  position  of  dj  is  at  the  intersection  of  cone-III  (axis  bj ,  angle  63)  and 

cone-IV  (axis  d[ ,  angle  an)  and  is  denoted  as  d"  in  the  Figure  4-4  (b).  If  the  two  cones 
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intersect  the  angle  of  rotation  about  d|  can  be  computed  as  /(d3,d[,b3,93) .  The 
cones  intersect  and  the  solution  exists  only  if  63  satisfies  the  following  condition. 

|Gi3-ai3|<e3<ei3  +  ai3  (4.17) 
In  the  above  equation,  813  is  the  angle  between  dj  and     ,  0<  0i3  <180°,  au  is  the  angle 

between  d,and  83,  0<  an  <180°.  The  transformation  matrix  required  to  impose  this 
complete  angular  constraint  set  can  be  expressed  as 

[Tn]  =  [T,(b3,<t)3)][T,(d;,(t),)][T,(n,(t),)]  (4.18) 

Sequential  Imposition  of  Distance  Constraints 

Distance  constraints  specify  distance  between  or  coincidence  of  points  (vertex  and 
reference  points),  lines  (edge  or  axis)  and  planes  (faces  or  reference  planes).  We  do  not 
consider  curves  and  surfaces  since  such  entities  can  be  positioned  with  respect  to  each 
other  by  specifying  constraints  between  reference  points,  lines  or  planes  that  are  attached 
to  the  curve  or  surface.  The  possible  distance  constraints  between  points,  lines  and  planes 
are  listed  in  table  4-1. 
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Table  4-1 :  Distance  constraints 


Distance  Constraint 

Acronym 

fixed 

Translation 
Space 

1. 

 -■  -■ — TJ  

Vertex  (point)-  Vertex  (point)  Coincident 

(VVC) 

-1 
J 

point 

2. 

Line  -  Line  Coincident 

L 

line 

J. 

Plane  -  Plane  Coincident 

(PPC) 

1 
1 

plane 

A 

4. 

Vertex  (point)  -  Line  Coincident 

(VLC) 

L 

line 

C 

J, 

Vertex  (point)  -  Plane  Coincident 

(VFC) 

1 
1 

plane 

O. 

Line  -  Plane  Coincident 

(LFC) 

1 
1 

plane 

7 

vertex  \point)  to  venex  ypoini)  uisiance 

(  V  \  U) 

1 
1 

spnere 

8. 

Vertex  (point)  to  Line  Distance 

(VLD) 

1 

cylinder 

9. 

Vertex  (point)  to  Plane  Distance 

(VPD) 

1 

plane 

10. 

Line  to  Line  Distance 

(LLD) 

1 

cylinder 

11. 

Line  to  Plane  Distance 

(LPD) 

1 

plane 

12. 

Plane  to  Plane  Distance 

(PPD) 

1 

plane 

The  constraints  that  are  highlighted  in  bold-italics  in  table  4-1  are  the  five 


fundamental  distance  constraints  since  all  the  others  can  be  converted  to  one  of  these 
five.  The  plane-plane  coincident  (PPC)  and  distance  (PPD)  constraints  assume  that  the 
two  planes  are  parallel  therefore  they  are  identical  to  the  vertex  plane  coincident  (VPC) 
and  distance  (VPD)  constraints  respectively  and  can  be  imposed  by  using  any  point  on 
one  of  the  planes  as  the  vertex.  Similarly,  the  line-plane  coincident  (LPC)  and  distance 
(LPD)  constraints  assume  that  the  line  is  parallel  to  the  plane  and  hence  they  are 
equivalent  to  the  VPC  and  VPD  respectively.  There  are  two  cases  for  LLC  and  LLD.  If 
the  two  lines  are  parallel  to  each  other,  LLC  and  LLD  are  equivalent  to  VLC  and  VLD 
respectively  and  can  be  imposed  by  using  any  point  on  one  of  line  as  the  vertex.  If  the 
two  lines  are  not  parallel,  LLC  and  LLD  can  be  imposed  by  constructing  a  plane  parallel 
to  one  line  that  passes  through  the  other  line  and  then  making  the  first  line  coincident 
with  the  plane  for  LLC  or  placing  it  at  the  specified  distance  in  case  of  LLD.  Therefore, 
LLC  and  LLD  are  equivalent  to  VPC  and  VPD  respectively.  In  this  case,  LLC  ensures 
that  they  intersect  and  LLD  makes  the  shortest  distance  between  the  lines  equal  to  the 
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specified  distance.  Finally,  since  VPC  is  actually  a  special  case  for  VPD  it  does  not  need 
to  be  considered  separately.  Thus  we  have  only  five  independent  constraints  WC,  VLC, 
VVD,  VLD  and  VPD.  Among  these  five  fundamental  constraints,  the  vertex- vertex 
coincident  (VVC)  is  the  only  constraint  that  fixes  all  the  translation  degrees  of  fi-eedom. 
This  constraint  is  easy  to  impose  since  the  required  translation  is  the  vector  fi-om  one 
vertex  to  the  other.  The  remaining  four  constraints  could  be  imposed  in  many  different 
sequences  but  we  will  use  the  preference  order  VLC^  VPD^  VVD^  VLD.  Notice  that 
in  all  four  constraints,  the  first  entity  is  a  vertex  (point)  while  the  second  entity  may  be 
point,  line  or  plane.  We  will  refer  to  the  point  corresponding  to  this  vertex  as  the 
constraint  point.  The  entity  that  belongs  to  the  dependent  object  will  be  referred  to  as  the 
dependent  entity  and  the  entity  attached  to  the  based  object  is  the  base  entity.  The 
possible  combinations  in  which  these  constraints  could  be  applied  are  listed  in  table  4-2. 

Table  4-2:  Possible  complete  distance  constraint  sets 


VLC  VPD 

VLC  ^  WD 

VLC  ^  VLD 

VPD     VPD  VPD 

VPD  ^  VPD  ^  VVD 

VPD     VPD  ->  VLD 

VPD  ^  VVD  ^  VVD 

VPD     VVD  ->  VLD 

VPD     VLD  VLD 

WD     WD  WD 

WD     VVD  VLD 

WD     VLD  ^  VLD 

VLD     VLD  ->  VLD 

Translation  Space 

If  previously  imposed  angular  constraints  have  removed  all  the  rotational  degrees 
of  freedom  of  the  object,  it  can  only  be  translated  to  impose  the  distance  constraints. 
After  a  distance  constraint  has  been  imposed  the  object  can  still  be  translated  without 
violating  this  constraint  if  it  has  any  degrees  of  fi-eedom  remaining. 
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Deflnition-I.  The  translation  space  of  a  point  on  the  dependent  object  with  respect 
to  a  distance  constraint  is  the  set  of  points  to  which  the  point  can  be  moved  to  by 
translating  the  dependent  object  without  violating  the  constraint. 

Definition-II:  The  translation  space  of  a  distance  constraint  is  the  set  of  points  to 
which  points  on  the  dependent  entity  can  be  moved  to  by  translating  the  dependent  object 
without  violating  the  constraint. 

Alternatively,  we  could  say  that  the  translation  space  of  the  constraint  is  the 
translation  space  of  any  point  on  the  dependent  entity.  The  geometry  of  the  translation 
space  of  various  constraints  is  described  in  table-Ill  as  point,  line,  plane,  sphere  or 
cylinder.  It  is  convenient  to  represent  the  translation  space  (curve/surface)  using 
parametric  equations  of  the  form  R{pi)  (curve)  or  R(pi,p2)  (surface),  where  pi  and  p2  are 
parameters  and  R  is  a  position  vector  of  points  on  the  translation  space  curve  or  surface. 
Below  we  list  the  translation  space  of  the  VLC,  VPD,  VVD  and  VLD  constraints. 


Table  4-3 :  Translation  space  expressed  using  parametric  equations 


Constraint 

Type 

Parametric  equation  for  translation  space 

1. 

VLC 

Line 

R'(M)  =  P„+Md 

2. 

VPD 

Plane 

R''(s,/)  =  P„  +5a  +  /b 

3. 

WD 

Sphere 

R''(<t),9)  =  C„  +/;(sin(t)COs6i+sin(t)sin6j  +  cos(t)k) 

4. 

VLD 

Cylinder 

R'^  (z,  (p)  =  P,  +  zl  +  r.  (cos  9111  +  sin  cpii) 

In  table  4-3,  R\u),  R^(5,/),  R^((p,0)  and  R%,^i9)  represent  the  translation  space  of 
the  constraints  VLC  (line),  VPD  (plane),  VVD  (sphere)  and  VLD  (cylinder)  respectively. 
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Po  is  a  point  on  the  translation  space.  For  the  VLC  constraint,  the  translation  space  is  a 
line  since  points  on  the  dependent  entity  can  only  move  along  a  line  once  this  constraint 
has  been  imposed.  This  is  true  regardless  of  whether  the  vertex  or  the  line  is  the 
dependent  entity.  In  table  4-3,  d  is  a  unit  vector  parallel  to  the  line  in  the  VLC  constraint. 
Similarly,  for  the  VPD  constraint,  the  translation  space  is  a  plane  since  the  dependent 
entity  (whether  it  is  the  vertex  or  the  plane)  can  only  move  along  a  plane  that  is  parallel  to 
and  at  the  desired  distance  from  the  plane  referred  to  in  the  constraint  definition,  a  and 

b  are  mutually  perpendicular  vectors  parallel  to  this  plane.  The  translation  space  is  a 
sphere  for  the  VVD  constraint  whose  center  Co  is  the  vertex  on  the  base  object  and  radius 
is  equal  to  the  distance  between  the  points  specified  by  the  constraint.  Finally,  in  the 
VLD  constraint  has  cylindrical  translation  space  whose  axis  i  is  a  unit  vector  parallel  to 
the  line  and  radius  rc  is  equal  to  the  specified  distance  between  point  and  edge.  lii  and  n 

are  mutually  perpendicular  unit  vectors  that  are  also  perpendicular  to  1 .  P,  is  the  point 
corresponding  to  the  vertex  if  the  vertex  is  the  base  entity  otherwise  if  the  line  is  the  base 
entity  then  Pa  is  a  point  on  this  line. 

Theorem:  If  R^(u)  is  the  translation  space  of  a  point  Po  (on  the  dependent  object) 
with  respect  to  a  constraint  X,  then  the  translation  space  of  any  other  point  P  on  the 
dependent  object  with  respect  to  this  constraint  is  R^(u)  =  R^(u)  +  (P-Po),  assuming 

that  previously  imposed  angular  constraints  have  eliminated  all  the  rotational  degrees  of 
freedom  of  the  dependent  object. 

The  above  theorem  is  based  on  the  observation,  that  since  the  angular  degrees  of 
freedom  have  been  removed  when  one  point  on  the  object  is  translated  along  a  curve  or 
surface  every  other  point  on  the  object  moves  along  an  offset  curve  or  surface  of  identical 
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shape.  More  formally,  after  the  constraint  X  has  been  imposed,  the  point  Po  is  on  its 
translation  space  R^.  Therefore,  there  exists  a  value  of  the  parameters  uo  for  which 
R^(uo)  =  Po.  Substituting  u  =  uo  in  the  definition  of  Rp  (u),  we  get      (u„)  =  P  showing 

that  P  is  on  Rp  .  When  P  is  moved  to  another  point  on  Rp ,  if  Po  stays  on  its  translation 

space  R^  (so  that  constraint  X  is  not  violated),  then  Rp  is  the  translation  space  of  P.  Let 

the  dependent  object  be  translated  such  that  P  is  moved  to  another  point  P'  =  Rp(u,) 

which  is  on  R^ .  During  this  translation  Po  is  moved  to  P^'  =  Pq  +  (P'  -  P)  We  can  show 

that  Pq  is  on  R^  as  follows: 

P'  =  R; (u,)  =  R^(u,)  +  (P  -  Po)  ^  R^(u,)  =  Po  +  (P'  -  P)  =  Po'  . 

Imposing  the  First  Distance  Constraint 

Let  R^^^(pi,p2)  (surface)  or  R^^V^)  (curve)  represent  the  translation  space  of  the 
first  constraint  and  Pi  be  the  position  of  a  point  on  the  dependent  entity.  To  impose  the 
constraint,  the  dependent  object  must  be  translated  such  that  the  point  Pi  lies  on  the 
translation  space  of  the  constraint.  The  translation  vector  required  is  a  vector  fi^om  Pi  to 
any  point  on  the  translation  space.  Arbitrary  values  could  be  assigned  to  the  parameters 
pi  and  p2  to  compute  a  point  on  the  translation  space.  Alternately,  the  point  Pi  can  be 
projected  on  to  the  translation  space  to  minimize  the  translation  distance.  This  can  be 
achieved  by  solving  the  following  optimization  problem  numerically. 

mn\r,-R''\p„p,)f  (4.19) 
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If  p[  and  p[  are  the  optimal  values  of  the  parameters  computed  from  the  above 
minimization,  the  translation  vector  (t)  can  be  computed  as 

t  =  (P;  -  P, )  where,  P.*  =  R^"  {p, ,  p\ )  (4.20) 

Imposing  the  Intermediate  Distance  Constraint 

The  second  constraint  can  be  considered  an  intermediate  distance  constraint  if  the 
first  and  second  constraint  fix  only  one  degree  of  freedom  each  so  that  one  more 
constraint  is  required  to  fully  constraint  the  dependent  object.  Let  R^^\qi,q2)  (surface)  be 
the  translation  space  of  the  intermediate  distance  constraint  and  P2  be  a  point  on  the 
dependent  entity  of  this  constraint.  Note  that  since  the  first  constraint  has  already  been 
imposed,  P2  is  a  point  on  the  transformed  position  of  the  dependent  entity.  To  impose  the 
intermediate  constraint,  Pj  "^'^^t  be  translated  to  a  point  on  the  intersection  of  its 
translation  spaces  with  respect  to  the  first  constraint  and  the  translation  space  of  the 
second  constraint.  Using  the  theorem  stated  earlier,  the  translation  space  of  P2  with 

respect  to  the  first  constraint  is  R^"(/?,,/72)  +  (P2  -P,*),  because  P,'  is  the  transformed 
position  of  P,  after  imposing  the  first  constraint.  The  intersection  curve  can  be 
represented  by  the  equation, 

R''\p„p,)  +  iV,-V:)  =  R''\q„q,)  (4.21) 
In  equation  (4-21)  we  have  three  simultaneous  equations  in  four  unknown  parameters  and 
therefore  it  represents  a  curve.  A  point  on  this  curve  could  be  found  by  projecting  the 
constraint  point  Pj  on  to  this  curve.  The  projected  point  can  be  found  by  the  following 
optimization: 
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|2 


P3-R^^H^„^.)1  (4.22) 


Min 

subject  to:R''\p,,p,)  +  (p,-F;)=R'''  {q, ,  q, ) 


If  p[,P2,q\  and  q\  are  optimal  values  of  the  parameters  for  the  above  optimization 
problem,  then  the  projected  point  is 

p;  =  R^'^(/>>;)+(p,-p;)=R^^H<,^;)  (4.23) 

The  translation  required  to  impose  the  second  constraint  is  t  =  (Pj  -  P2 ) 


Imposing  the  Last  Distance  Constraint 

Let  Rp'^ij-^j^)  be  the  translation  space  of  the  last  constraint  and  P3  be  a  point  on 
the  transformed  position  of  the  dependent  entity  after  imposing  the  earlier  constraints. 
Translating  P3  to  the  intersection  of  the  translation  spaces  of  P3  with  respect  to  all  the 
distance  constraints  will  impose  the  final  constraint  without  violating  the  previous 
constraints.  The  translation  spaces  of  P3  with  respect  to  the  first  two  distance  constraints 

are  R^'^(/?,,/?2)  +  (P3 -P;)andR''^(^,,^2)  +  (P3 -^2)  where  P/and  P*  are  transformed 
positions  of  Pi  and  P2  after  imposing  the  first  and  second  constraints.  The  new  position  of 
P3  can  be  computed  by  solving  simultaneously  the  following  six  equations 

R^'^(/.„/>,)  +  (P3  -P;)  =  R^^>(^„^,)  +  (P3  -P;)  =  R''\r,,r,)  (4.24) 

If  /?*,     ql,  q\,  r*  and  r*  are  the  solutions  for  the  six  parameters  obtained  by  solving  the 
above  equations,  then  new  position  of  P3  can  be  computed  as 

Yl=R^'\r:,r;)  or  P;  =K'\plp\)^{V,-V:)^K^'\qlq\)  +  iV,-V;) 

(4.25) 
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Translating  the  dependent  object  along  the  translation  vector  t  =  (P3 -Pj)  will 
impose  the  final  constraint.  The  method  described  above  is  very  general  and  applicable  to 
any  kind  of  constraint  if  implemented  using  numerical  procedures  for  solving  equations 
or  optimization.  However,  in  practice,  the  implementation  can  be  made  more  efficient  by 
using  analytical  solutions  for  intersection  of  translation  spaces  whenever  possible. 

Analytical  Methods  for  Sequentially  Imposing  Distance  Constraints 

In  most  cases,  the  procedure  presented  above  can  be  performed  analytically 
because  we  can  find  analytical  solution  for  intersection  of  translation  spaces.  Intersection 
of  two  translation  spaces  that  are  both  two-dimensional  yields  a  one  dimensional 
translation  space.  Intersection  of  a  one-dimensional  translation  space  and  a  two- 
dimensional  translation  space  yields  a  zero  dimensional  translation  space,  which  could  be 
one  or  two  points.  Other  than  applying  numerical  techniques  to  find  a  point  on  this 
translation  space,  one  can  often  find  an  analytical  solution  that  can  be  expressed  as  a 
parametric  equation  of  a  curve.  Below  we  list  a  few  such  cases. 

(1)  Plane  -  Plane  intersection 

Let  P, (5,, /))  =  /?, +5ja, and  Pj{s^,t^)  =  +  5^012 +^2^2  planar 
translation  spaces  whose  normal  vectors  are  «,  and  .  The  intersection  of  these  two 
planes  is  a  line  whose  parametric  equation  is: 

L=Lo+sd  (4.26) 
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^1  and     =  a, ,  if  we  select  the  vectors  on 


the  plane  as  a,  =  «,  x     and  b^=a^xh^,  a^=     and  b^=a^y.fi^. 


(2)  Plane  and  Sphere  intersection 

Intersection  of  plane  and  sphere  leads  to  a  circular  translation  space.  Let 
P  {s  ,t  )-R  +  s  a  +t  b   be  the  planar  translation  space  whose  normal  vector  is 

axb  =  h  and  base  point  R  and  the  spherical  translation  space  be  S  with  center  C  and 
radius  r.  The  intersection  circle  can  be  expressed  as: 


where,  7?^  =C -C)»«]/j  is  the  center  of  the  circle  and  its  radius  is 


(3)  Plane  and  Cylinder  intersection; 

A  circular  translation  space  is  obtained  by  the  intersection  of  a  planar  and 
cylindrical  translation  spaces.  Let  the  planar  space  be  P(s  ,t  )  =  R  +s  a  +t  b  with 


CR  (9)  =  R,+r^  {a  cos  9  +  *  sin  9) 


(4.27) 


r.  =  ^|r'-[(R-C)•hf 


normal  n-axb  and  the  cylindrical  space  be  C(u,0)  =  R^  +u  d  +  r^{m cosO  + 1  sin  0) . 


The  intersection  ellipse  is: 


E{e)  =  C,  +m^Acos0  +  l^Bsm0 


(4.28) 
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The  center  of  the  elHpse  is  C^=R^+ud ,  the  direction  of  the  major  and  minor 
axes  are  l^  -hxd  and  m^-l^xh  while  the  magnitudes  are  A  =  ^    ^  and  B  =  r^  where 

d»h 


(4)  Sphere  and  Sphere  intersection: 

Intersection  of  two  spherical  translation  spaces  whose  centers  are  Q  and     ,  and 
radii  are  ri  and  r^,  yields  the  following  circular  translation  space. 

CR(e)  =  ^,+r,(acose  +  *%ine)  (4.29) 

^2  _  ^  2  J 

The  center  of  the  circle  is     =    '  ?  ^  (C;  -  C,)  + — (Q  +  Q)  and  the  radius  of 

2 1  (^^C^2 1  2 


the  circle  is  n.  =  Jr,^  --^ 


'l'^2 


a  and  b  are  two  unit 


vectors  on  the  circle  plane  such  that  axb  =  C^Cj  and  bxC^C^  =  a . 
(5)  Line  and  Plane  intersection: 

Intersection  of  a  line  translation  space  and  a  planar  translation  space  yields  an 

intersection  point.  Let  the  line  translation  space  be  Z=Zo+M^and  the  planar 

translation  space  he  P  (s  ,t  )  =  R  +s  a  +  /  ^  ,  the  normal  of  the  plane  is  h  =  axb.  The 
intersection  point  is: 


(4,30) 


(6)  Line  and  Sphere  intersection 
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Intersection  of  a  line  translation  space  and  a  spherical  translation  space  consists  of 
one  or  two  points.  Let  the  line  translation  space  be  L  =  L^+u  d and  the  spherical 

translation  space  be  with  center  C  and  radius  r,  the  intersection  points  are: 


{C-L,)*d±^r'-{\C-L,  \y-[{C-L,).df 


(4.31) 


(7)  Circle  and  Sphere  intersection 

Intersection  of  a  circular  translation  space  and  a  spherical  translation  space 

consists  of  one  or  two  points.  Let  the  circular  translation  space  be  with  center  C^,  normal 

h  and  radius  Vc  and  the  spherical  translation  space  be  with  center  C^^  and  radius  r,.  A 

plane  with  normal  h  and  base  point  has  an  intersection  circle  with  Sphere.  Let  C  and 
r  be  the  center  and  the  radius  of  the  intersection  circle  respectively,  we  have 

c=c:^+[(c:-c^).A2]/5 

C,C  =C  -Q 

cc       .  - 

Let  /  =  —  —  and  /w  =  /  x  « ,  the  intersection  point  is: 

IQC| 

P=C^+C^P  (4.32) 

Where 


C^P  =rcose  i  +r  sine 


c 


m 


where 
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6'  =  ±cos"'  ^     '  ^ 


2/-JQCI 

The  ±  means  there  are  possibly  two  intersection  points. 


Application  to  Solid  modeling 

The  appHcation  of  the  sequential  constraint  imposition  method  to  dimension- 
driven  solid  modeling  is  illustrated  in  this  section  using  an  example.  We  assume  that  the 
solid  is  constructed  by  combining  primitives  or  simple  solids  using  regularized  Boolean 
operations.  The  solids  to  be  combined  are  positioned  with  respect  to  each  other  using 
position  constraints  between  vertices,  edges  and  faces  of  the  solids.  Figure  4-5  shows  the 
partial  CSG  tree  of  the  final  solid  model  S?.  For  brevity,  the  details  of  how  the  solids  S2, 
S4  and  Se  were  created  by  combining  simple  primitives  are  not  shown  in  the  figure.  At 
each  node  of  the  CSG  tree  we  assume  that  the  solid  on  the  left  is  the  base  object  and  the 
solid  on  the  right  is  the  dependent  object.  Only  the  Boolean  operation  at  each  node  is 
shown  in  the  CSG  tree,  however,  it  is  assumed  that  a  constraint  set  is  also  defined  for 
each  node,  The  dependent  object  is  positioned  with  respect  to  the  base  by  sequenfially 
imposing  the  constraints  before  the  Boolean  operation  is  applied  to  combine  them. 
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At  node  1,  the  union  operation  is  applied  to  create  S3  by  combining  the  soHds  Si 
and  S2.  The  dependent  object  S2  was  oriented  with  respect  to  Si  by  making  the  faces  fi 
and  fz  parallel  (using  a  PV  constraint  that  makes  the  normal  vectors  of  the  faces  parallel). 
Notice  that  the  normal  vectors  of  these  faces  are  in  opposite  directions  in  this  figure.  This 
can  be  specified  using  an  appropriate  flag  in  the  array  denoted  as  d  in  equation  (4.2).  Due 
to  symmetry  of  Si,  rotation  of  S2  about  the  axis  Al  does  not  make  any  difference  in  the 
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geometry  of  S3.  Therefore,  it  is  not  necessary  to  fully  constrain  the  rotational  degrees  of 
freedom  of  S2.  The  three  translation  constraints  used  at  node  1  are: 

(i)  Line  to  Plane  Coincident  (LPC)  to  make  plane  fi  coincident  with  axis  Al, 

(ii)  Plane  to  Plane  Distance  (PPD)  to  specify  distance  between  the  faces  f2  and  fl  and 

(iii)  Line  to  Plane  Distance  (LPD)  to  specify  distance  between  face  f4  and  axis  Al. 
Notice  that  the  translation  set  used  at  node  1  is  VPD->VPD^VPD,  since  all  three 

constraints  used  above  (LPC,  PPD  and  LPD)  reduce  to  VPD  as  discussed  in  table  4-1.  At 
node  2,  Ss  is  created  by  the  union  of  S3  and  S4  after  positioning  S4  with  respect  to  S3 
using  the  following  angular  and  translation  constraint  sets. 
Complete  angular  constraint  set  (CACS-II): 

(i)  Vector  to  vector  angle  (WA)  constraint  to  specify  an  angle  of  135°  between  the 
normal  vectors  of  face  f?  and  face  ft, 

(ii)  Vector  to  vector  angle  (VVA)  constraint  to  specify  that  the  normal  vector  of  face 
f?  is  perpendicular  to  the  normal  of  face  fi. 

(iii)  Vector  to  vector  angle  (VVA)  constraint  to  specify  that  the  faces  fl  and  f6  (or 
their  normal  vectors)  are  at  an  angle  of  30°. 

Complete  translation  constraint  set  (VPD->VPD^VPD): 

(i)  Line  to  plane  coincident  (LPC)  constraint  to  make  the  axis  Al  coincident  with  the 
face  f7, 

(ii)  Line  to  line  coincident  (LLC)  to  make  the  axis  Al  intersect  the  axis  A4  and 

(iii)  Line  to  plane  distance  (LPD)  to  specify  the  distance  between  the  axis  A4  and  the 
face  fl . 
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Finally,  the  solid  S?  is  created  by  subtracting  the  solid  Se  from  S5  after  position  Se 
with  respect  to  S5  using  the  following  complete  constraint  sets. 
Complete  angular  constraint  set  (CACS-I): 

(i)  Parallel  vectors  (PV)  constraint  to  make  the  faces  (or  the  normal  vectors  of  the 
faces)  fi  and  fg  parallel. 

(ii)  Vector  to  vector  angle  (WA)  constraint  to  make  the  angle  between  the  faces  f4 
and  fg  equal  to  zero. 

Complete  translation  constraint  set  (VEC^VPD): 

(i)  Line  to  line  coincident  (LLC)  constraint  to  make  the  axes  Al  and  A5  coincident, 

(ii)  Plane  to  plane  coincident  (PPC)  to  make  the  planes  fl  and  f8  coincident. 

When  defining  the  solid  interactively  using  a  CAD  system,  each  constraint  can  be 
imposed  as  soon  as  they  are  defined  using  the  approach  described  here.  To  enable 
dimension-driven  editing  the  constraints  used  at  each  node  must  be  stored  as  introduced 
in  chapter  3  so  that  when  any  dimensions  are  changed,  the  solid  can  be  updated  by 
recreating  the  primitives,  if  necessary,  re-imposing  the  constraints  sequentially  to  position 
them  and  then  applying  the  Boolean  operations.  For  the  example  shown  in  figure  4-5, 
many  alternate  sets  of  constraints  could  be  used  to  position  the  primitives  and 
intermediate  solids  relative  to  each  other  exactly  as  shown.  However,  when  a  dimension 
needs  to  be  changed  and  the  solid  needs  to  be  updated,  different  sets  of  constraints  may 
lead  to  different  final  solids.  Therefore,  the  constraints  must  be  carefully  selected  to 
ensure  that  they  represent  the  design  intent. 


66 

Discussion 

A  method  for  sequentially  imposing  position  constraints  to  position  rigid  bodies 
with  respect  to  each  other  is  described.  This  method  can  be  used  in  solid  modeling  and 
assembly  modeling  to  position  geometric  primitives  or  parts  with  respect  to  each  other.  In 
comparison  to  variational  geometry,  the  main  advantage  of  this  approach  is  that  this 
approach  does  not  require  solving  all  the  constraints  used  in  the  model  simultaneously. 
The  constraints  between  each  pair  of  solids  in  a  Boolean  operation  can  be  imposed  one  at 
a  time  both  during  construction  and  when  the  model  is  updated  after  changing  any 
dimension.  The  sequential  constraint  imposition  can  therefore  be  thought  of  as  the 
decomposition  of  this  large  set  of  constraints  into  simpler  sets  that  are  solved  to  position 
rigid  bodies  or  primitives.  In  solid  modeling  application,  the  geometry  of  the  primitive 
itself  may  be  generated  either  by  pre-defined  procedures  or  by  variational  geometry. 
However,  even  when  variational  geometry  technique  is  used  to  update  the  primitives,  by 
sequentially  positioning  the  primitives,  we  can  fully  de-couple  the  position  constraints 
from  the  constraints  controlling  the  geometry  of  the  primitive. 

Availability  of  analytical  solutions  for  most  constraints  enables  fast  imposition  of 
constraints  both  during  construction  and  updating  of  the  solid  model.  Another  advantage 
of  using  analytic  solutions  is  that  all  solutions  are  available  when  there  are  multiple 
solutions.  Therefore,  the  user  can  interactively  select  the  desired  solution  based  on  the 
design  intent.  During  updating,  this  same  solution  can  be  automatically  reselected  by 
establishing  a  method  to  distinguish  between  these  solutions.  For  example,  if  the  designer 
selected  the  solution  corresponding  to  the  smallest  value  of  a  parameter,  then  this  solution 
can  be  used  during  regeneration  also.  In  most  cases,  we  are  also  able  to  derive  conditions 
for  existence  of  a  solution. 
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The  sequential  approach  described  here  is  only  applicable  when  all  rotational 
degrees  of  freedom  are  eliminated  using  angular  constraints.  As  mentioned  earlier,  this  is 
not  a  serious  limitation  for  applications  in  solid  and  assembly  modeling.  We  are  currently 
investigating  the  possibility  of  extending  the  method  to  handle  constraint  sets  that  do  not 
include  a  complete  angular  constraint  set. 

A  preferred  sequence  was  adopted  for  imposing  the  position  constraints  in  order 
to  minimize  the  number  of  different  cases  that  need  to  be  considered.  In  addition,  the 
preferred  sequence  was  selected  with  the  goal  of  simplifying  the  procedure  for  imposing 
the  constraints  without  violating  previous  constraints.  However,  this  does  not  imply  that 
one  has  to  specify  constraints  in  this  sequence  when  building  a  solid  model  interactively. 
For  example,  the  designer  may  prefer  to  specify  the  distance  constraints  before  the  angle 
constraints.  This  does  not  pose  a  problem  since  the  object  can  be  positioned  (or 
repositioned,  if  a  wrong  sequence  was  initially  used)  in  the  right  sequence  after  the  object 
is  fully  constrained. 


CHAPTER  5 

NON-REGENERATIVE  UPDATING  ALGORITHM 
Overview 

In  this  chapter,  a  detailed  updating  algorithm  for  dimension-driven  solid  modeling 
is  presented,  which  can  directly  modify  the  B-rep  of  the  original  solid  in  response  to 
dimension  change.  The  organization  of  the  rest  of  this  chapter  is  as  follows.  In  first 
section,  we  discuss  problems  of  current  updating  algorithm  for  history-based  parametric 
modeling.  In  second  section,  the  concept  of  reusing  geometric  entity  in  Boolean  operation 
is  introduced.  In  third  and  fourth  sections,  we  introduce  how  to  modify  primitive  and 
reposition  geometry  in  updating  process.  In  fifth  section,  an  algorithm  for  updating  solid 
models  without  topology  change  is  presented  first,  followed  by  description  of  a  more 
robust  updating  algorithm,  which  can  handle  topology  change  during  updating.  In  final 
section,  We  discuss  how  to  take  advantage  of  local  modification  to  further  improve  the 
proposed  algorithm. 


Analysis  of  Existing  Algorithm 

Based  on  the  history  representation  presented  in  chapter  3,  the  "rollback" 
algorithm  in  parametric  modeling  can  be  used  to  regenerate  the  final  geometric  model  in 
case  of  dimension  change.  In  the  example  shown  in  figure  3-1,  the  model  generation 
history  can  be  represented  by  the  following  feasible  sequence: 
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{Cp,  ,  Cp2  ,  Pos(P2,  PI),  P1-'P2  ,  Cp3  ,  Pos(P3,  Gl),  Glu'P3,  Cp4  ,  Cps,  Pos(P5,  P4),  P4u*P5, 
Pos(G3,G2),G2-'G3}  (5.1) 

where  Cx  represents  an  operation  that  create  primitive  x  and  Pos(X,  Y)  represents 
an  operation  that  position  geometry  X  with  respect  to  Y  by  imposing  a  given  complete 
constraint  set.  If  any  dimensions  of  the  primitive  P2  are  changed,  the  algorithm  will 
rollback  the  state  after  operation  Cpi  and  redo  the  following  operations: 

{Cp2 ,  Pos(P2,  PI),  PI-*P2  ,  Cp3 ,  Pos(P3,  Gl),  Glu*P3,  Cp4  ,  Cpj,  Pos(P5,  P4),  P4u*P5,  Pos(G3, 
G2),G2-*G3}  (5.2) 

which  is  almost  the  same  as  recreating  the  whole  model.  Obviously,  this  kind  of  updating 
algorithm  based  on  reconstruction  is  inefficient.  By  carefully  studying  this  example,  we 
find  that  P2  is  a  relatively  independent  primitive  in  the  model.  The  only  geometry 
involved  in  the  Boolean  operation  and  position  constraints  related  to  P2  is  primitive  PI. 
Intuitively,  local  modification  related  to  PI  and  P2  should  be  sufficient  to  update  the 
solid.  The  other  disadvantage  of  the  rollback  algorithm  is  that  the  primitives  have  to  be 
recreated  and  Boolean  operations  have  to  be  re-applied  to  the  model.  Moreover, 
persistent  naming  is  needed  for  this  kind  of  updating  to  identify  the  topological  entities 
referred  to  in  the  history.  Persistent  naming  is  slow  and  unreliable.  Three  possible  ways 
to  improve  the  updating  efficiency  are;  (i)  making  modification  more  locally,  (ii) 
avoiding  unnecessary  recreation  of  geometric  entities  and  (iii)  trying  to  avoid  using 
persistent  naming.  Based  on  these  ideas,  a  non-regenerative  updating  algorithm,  which 
can  directly  modify  the  B-rep,  is  designed  and  presented  in  this  chapter. 
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Geometric  Entity  Reuse  in  Boolean  Attach  Operation 

Topological  entities  in  B-rep  include  vertices,  edges,  faces  and  other  optional  high 
level  entities  such  as  shell  and  body.  Among  these  topological  entity  types,  vertex,  edge 
and  face  have  their  own  corresponding  geometric  entities,  which  are  point  (coordinates), 
curve  and  surface  respectively.  For  a  set  of  topological  entities  E,  the  notation  G(E)  v^all 
be  used  to  represent  the  geometric  entities  corresponding  to  E.  During  the  Boolean 
operation,  we  can  reuse  geometric  entities  of  the  two  original  B-reps  involved  in  the 
Boolean  operation  whenever  possible  instead  of  creating  new  geometric  entities.  In  other 
words,  some  topological  entities  in  the  resulting  B-rep  can  share  geometric  entities  with 
topological  entities  in  the  two  original  B-reps.  The  purpose  of  geometric  entity  reuse  is 
that  any  changes  in  geometric  models  of  the  two  children  will  directly  cause  the 
corresponding  changes  in  their  reused  entities  in  the  solid  model  of  the  parent.  Figure  5-1 
shows  a  sub-tree  representing  a  Boolean  operation  in  a  CSG  tree.  /  is  an  internal  node 
with  left  child  L,  and  right  child  R,  of  the  CSG  tree.  Let  Sr,  and  Su  be  the  geometric 
model  (B-rep)  represented  by  nodes  Rt  and  Z,  respectively  and  5",  be  the  geometric  model 
(B-rep)  generated  by  attaching  Sm  to  Su  with  a  Boolean  operation  in  node  /.  Sm  and  Su 
can  be  primitive  or  intermediate  geometry. 


Figure  5-1;  Boolean  operation  with  geometry  entity  reuse 
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Si  can  be  further  expressed  as  a  set  of  B-rep  entities  as  following: 

S,  =  {  Efi  ,  Ee. ,  Evi}  (5-3) 

where  Ef„  Eei  and  Evi  are  face,  edge  and  vertex  sets  respectively.  Similarly,  Sr  and  Sl  can 

be  expressed  as: 

SRi={EfRi,EeR,,EvR,}  (5.4) 

SLi={EfLi,EeL,,EvL,}  (5-5) 

All  topological  entities  in  Si  are  newly  created  during  the  Boolean  operation.  The 
reusable  entities  are  the  corresponding  geometric  entities  of  faces,  edges  and  vertices. 
Considering  entity  reuse,  Si  can  also  be  written  as: 

Si={  Ei"^Ei^"^Ei^Ei^}  (5.6) 

where  Ei*"*  is  a  set  of  intersection  entities  including  intersection  vertex  set  Evi""  and 
intersection  edge  set  Eei""*: 

E,'"*={Evr,Ee*"'}  (5.7) 
An  intersection  vertex  is  the  intersection  of  one  edge  from  SRi  or  Su  and  one  face  from 
SLi  and  SRi.  The  edge  and  face  involved  in  the  intersection  are  usually  called  intersecting 
edge  or  intersecting  face  respectively.  An  intersection  edge  is  the  intersection  of  one 
intersecting  face  from  SRi  and  one  intersecting  face  from  Su. 
Ei'"*"  is  a  set  of  subdivision  entities: 

Ei^""  =  {  EeR^"\  EeL^"^  EfR^"^  En."'  }  (5.8) 
where  Eer'"''  is  set  of  subdivision  edges  generated  by  subdividing  intersecting  edges  of 
SRi ,  EeL^"''  is  set  of  subdivision  edges  generated  by  subdividing  intersecting  edges  of  Sy  , 
EfR^"''  is  set  of  subdivision  faces  generated  by  subdividing  intersecting  faces  of  SRi  and 
EfL^''  is  set  of  subdivision  faces  generated  by  subdividing  intersecting  faces  of  Su  . 
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Ei^  is  a  set  of  topological  entities  sharing  geometric  entities  with  a  subset  (E^')  of 
the  topological  entities  in  Sri: 


E.^-  {Ef,^  Eei^Evi^} 

(5.9) 

E-Ri  -  (tfRi  ,  t-eRi  ,  t-vRi  / 

[p.  IV) 

G(Ei'')  =  G(ER,') 

(5.11) 

G(Efi'')  =  G(EfR,') 

(5.12) 

G{Ee^)  =  G(EeRi') 

(5.13) 

G(Evi'')  =  G(EvRi') 

(5.14) 

where  Eri  e  Sri  and 

EfRi'  c  EfRi ,  EeRi'  c  EeRi  and  EvRi'  Q  EvRi  (5.15) 
Similarly,  Ei^  is  a  set  of  topological  entities  sharing  geometric  entities  with  a  subset  (Eu) 
of  the  topological  entities  in  Su: 

Ei^={Ef,^Ee^Ev,^}  (5.16) 

EL'={EfL',EeL',EvL'}  (5-17) 


Figure  5-2;  3D  example  of  geometry  entity  reuse 
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Figure  5-2  shows  an  3D  example  of  Boolean  union  operation,  the  subsets  of 
topological  entities  in  two  children,  whose  corresponding  geometric  entities  are  shared  by 
resuhed  body,  are: 

Eri'  =  {all  vertices,  edges  and  faces  except  face  Fr2  and  vertices,  edges  on  this 


face}  (5.18) 

Eli'  =  {all  vertices,  edges  and  faces}  (5  19) 
and  the  topological  entities  with  newly  generated  geometric  entities  are: 

Ev'"'={Vi,,V,2,  V,3,Vi4}  (5.20) 

Ee"^={en,e,2,ei3,ei4}  (5.21) 

EeR^"*'={e.5,e,6,ei7,ei8}  (5.22) 

EeL"''={}  (5.23) 

Em"''={F,2,Fi3,Fi4,Fi5}  (5.24) 

Efl.^"''={Fn}  (5.25) 


where  Vii,  Vi2  ,  Via  and  Vi4  are  intersection  vertices,  Cii  ,  ei2  ,  ei3  and  ei4  are 
intersection  edges;  eis ,  eie ,  and  eig  are  subdivision  edges;  Fii ,  Fi2 ,  Fi3 ,  Fi4  and  Fij  are 
subdivision  faces.  Because  only  geometric  entities  are  reused  and  the  topological  entities 
of  child  nodes  keep  unchanged,  the  B-reps  Sri  and  Su  of  both  children  are  not  affected 
and  can  be  stored  in  their  own  nodes  after  the  Boolean  attach  operation. 


Modification  of  Primitive 

Dimension  changes  of  any  primitive  will  cause  the  modification  of  the  geometry 
of  all  entities  of  the  primitive  and  this  modification  will  propagate  in  bottom-up  direction 
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along  the  CSG  tree  and  cause  the  modification  of  the  final  geometric  model  due  to  the 
geometric  entity  reuse.  Therefore,  modification  of  primitives  is  the  starting  point  of 
updating  process.  As  mention  before,  each  type  of  primitive  can  be  parameterized  by  a 
set  of  independent  dimensions  (parameters)  and  has  its  own  predefined  geometric 
creation  procedure.  Similarly,  we  can  predefine  geometric  modification  procedure 
Mi(Axii,  Axin)  for  each  type  of  primitive,  where  i  is  the  type  and  Axy  is  the  change  of  j- 
th  dimension  and  n  is  the  number  of  independent  dimensions  of  this  type  of  primitive. 
The  procedure  Mi  modifies  the  geometry  of  all  the  entities  of  the  primitive  in  its  local 
coordinate  system  in  response  to  changes  in  the  independent  dimensions.  Figure  5-3 
shows  how  a  geometric  modification  procedure  for  rectangular  block  works. 


(a)  (b) 

Figure  5-3:  Modification  of  primitive 

There  are  three  independent  dimensions  for  rectangular  block  and  the  origin  of  the 
local  coordinate  system  is  at  vertex  1.  When  dimension  da  is  changed,  geometric 
modification  procedure  will  move  corresponding  points  (coordinates)  of  vertices  5,  6,  7 
and  8  a  distance  Ada  along  the  direction  1^5.  After  moving  the  vertices,  the  geometry  of 
affected  faces  and  edges  are  also  moved  to  their  new  positions.  In  this  modification 


75 

procedure,  only  the  positions  of  the  geometric  entities  are  changed  so  that  no  new  entities 
are  created.  For  a  primitive  soUd,  this  kind  of  geometric  modification  is  actually  a  non- 
regenerative  updating  process  In  a  complex  final  geometric  model,  a  primitive  could  be 
in  its  original  position  (figure  5-3  (a))  or  any  transformed  position  (Figure  5-3  (b)). 
However,  the  primitive  modification  always  happens  in  primitive  geometry's  local 
coordinate  system. 

Geometry  Repositioning  in  Updating  process 

For  a  sub-tree  shown  in  figure  5-1,  Sri  needs  to  be  repositioned  with  respect  to  Su 
in  two  cases.  The  first  case  is  that  the  value  of  at  least  one  position  constraint  in  node  i's 
constraint  set  is  changed.  For  example,  if  "face  to  face  angle  "  is  a  position  constraint  of 
the  constraint  set,  the  change  of  the  angle  value  will  cause  the  repositioning  of  right  child 
Sr].  The  second  case  is  that  the  geometric  entity  of  at  least  one  of  the  topological  entities 
involved  in  the  constraint  set  is  changed  or  relocated.  Any  modification  of  the  primitive 
or  a  change  in  constraints  of  the  sub-trees  Li  and  Ri  could  cause  the  second  case.  In  both 
cases,  a  series  of  transformation  matrices  can  be  obtained  using  the  method  for  imposing 
constraints  described  in  previous  chapter. 

Each  obtained  transformation  matrix  is  applied  to  the  sub-tree  Ri  in  a  recursive 
way  such  that  geometric  entities  of  all  primitives  (represented  by  leaf  nodes)  under  Sni 
and  all  entities  generated  in  the  construction  process  of  Sri  are  repositioned.  Suppose  Sj 
needs  to  be  repositioned  with  respect  to  its  left  sibling  and  the  transformation  matrix  T 
has  already  been  obtained  for  imposing  one  of  the  position  constraints  in  the  constraint 


76 

set  of  node  j's  parent,  the  recursive  algorithm  for  applying  this  transformation  to  node  j 
can  be  written  as; 

if  node  j  is  a  leaf  node,  G(Sj)  =  T  G(Sj) ; 
else 

GCEj*"')  =  T  G(Ej"") 
apply  T  to  node  Lj  and  Rj; 

Figure  5-4  Algorithm  for  geometry  repositioning 

If  node  j  is  a  leaf  node  and  represents  a  primitive,  the  transformation  matrix  is 
directly  applied  to  all  geometric  entities  G(Sj).  If  node  j  is  an  internal  node,  the 
transformation  matrix  is  directly  applied  to  geometric  entities  generated  in  the  Boolean 
operation  stored  in  this  node,  which  include  geometric  entities  G(Ej'"*)  of  all  intersection 
entities  Ej'"'  ,  and  then  the  transformation  matrix  is  recursively  applied  to  left  and  right 
children  of  node  j. 


Updating  Algorithm 

The  purpose  for  updating  Si  is  to  update  all  geometric  entities  in  G(Si)  and  add  to 
or  delete  some  topological  entities  from  Si  to  match  the  dimension  change.  In  case  there 
is  no  topology  change,  no  topological  entities  are  added  to  or  deleted  fi"om  set  Si  and  only 
the  geometric  entities  in  G(Si)  need  to  be  updated.  In  case  of  topology  change,  the  size 
and  members  of  all  subset  Ei"' ,  Ei^"**,  Ei^  ,  Ei^  of  Si  could  change.  In  this  case,  the  first 
thing  to  do  for  updating  is  to  delete  unnecessary  members  from  Si  and  create  and  add  new 
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topological  entities  to  Si  if  any.  The  second  thing  to  do  is  to  update  the  geometric  entities 
of  undeleted  entities  in  Si. 

Updating  Algorithm  (No  Topology  Change) 

Assuming  that  no  topology  change  occurs,  the  updating  algorithm  for  modifying 

Si  in  figure  5-1  can  be  expressed  as: 

1.  update  Sl,; 

2.  update  Sr,; 

3.  reposition  Sri; 

4.  update  GCE,"' )  and  G(Er'); 

Figure  5-5  Updating  algorithm  when  there  is  no  topology  change 

In  step  1,  if  subtree  Li  represents  a  primitive,  the  primitive  geometry  is  modified, 
otherwise,  the  updating  algorithm  will  be  recursively  applied  to  Li.  In  this  step,  G(Ei^)  is 
updated.  Similarly,  in  step  2,  G(Ei^)  is  updated.  In  step  3,  G(E^)  is  transformed  to  new 
position  if  necessary. 

Before  discussing  step  4,  a  little  bit  discussion  about  regularized  Boolean 
operation  is  necessary.  Many  algorithms  (Mortenson,  1985;  Mantyla,  1988;  Hoffmann, 
1989;  Lee  1999)  have  been  introduced  for  regularized  Boolean  operation.  Basically, 
Boolean  operation  involves  the  following  three  steps:  a)  intersection  computation  for 
creating  Ei'"^  b)  set  membership  classification  for  non-intersecting  entities  for  creating 
Ei^  and  Ei*^  and  c)  subdivision  computation  for  creating  Ei^"''.  To  create  Ei*"*,  we  have  to 
compare  each  pair  of  faces  from  EfRi  and  EfLi  and  then  compute  the  geometry  of 
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intersection  vertices  and  edges  for  tlie  intersecting  face  pairs.  To  create  entities  in  Ej'"'', 
we  have  to  perform  set  membership  classification  to  create  subdivision  edges  and 
construct  the  connectivity  among  intersection  edges,  subdivision  edges  and  edges  m  Ei 
and  Ei^  to  create  subdivision  faces.  In  the  initial  construction  process  of  the  soUd,  this 
process  is  inevitable.  However,  during  the  updating  process  a  large  amount  of  work  can 
be  saved  by  recording  some  information  during  initial  Boolean  operation. 

The  most  important  information  is  how  an  intersection  entity  was  computed.  This 
information  includes  the  entities  involved  in  the  intersection  and  can  be  easily  recorded 
as  the  attribute  of  an  intersection  entity.  For  example,  one  way  to  record  this  information 
is  to  record  two  faces  as  the  attribute  of  an  intersection  edge  generated  from  the 
intersection  of  these  two  faces.  Based  on  this  information,  we  can  update  GCEi*"'  )  and 
G(E,'"'')  in  step  4  of  updating  process  instead  of  redoing  the  Boolean  operation. 

For  an  intersection  entity,  if  the  entities  involved  in  the  intersection  are  known  we 
do  not  have  to  find  the  intersecting  face  pairs.  Instead,  we  can  directly  re-compute 
G(Ei'"Y  In  the  example  shown  in  figure  5-2,  36  (6x6)  face  intersecfion  computafions  are 
needed  when  redoing  Boolean  operation,  whereas  only  4  face  intersection  computations 
are  needed  when  updating  the  geometry  of  the  intersection  entities  based  on  the  recorded 
intersection  information. 

A  subdivision  entity  in  Ei'"**  usually  includes  lower  level  intersection  entity  in  Ej*"' 
,  lower  level  subdivision  entity  and  shares  part  of  its  geometry  with  subdivided  entity  in 
Sri  or  Su.  For  example,  a  subdivision  edge  usually  includes  one  or  two  intersection 
vertices  in  Ei'"'  and  shares  a  geometry  entity  (curve)  with  the  subdivided  edge  in  Sri  or 
Su-  A  subdivision  face  usually  includes  intersection  edges,  subdivision  edges  and  edges 
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in  Ei^  and  Ei^.  Subdivision  face  also  shares  geometry  entity  (surface)  with  the  subdivided 
face.  Based  on  this  discussion,  it  is  easy  to  see  that  G(Ei^"'')  depends  on  G(Ei'"' ),  G(SRi), 
G(SLi),  G(E^)  and  G(Ei^)  which  are  all  updated  after  the  first  three  steps  and  re- 
computing G(Ei'"').  Therefore,  G(Er^)  is  updated  automatically. 

Compared  to  the  rollback  algorithm  or  other  regenerative  updating  algorithms  that 
involve  redoing  Boolean  algorithm,  the  advantages  of  proposed  algorithm  is  obvious 
when  there  is  no  topology  change.  Firstly,  much  fewer  face  intersection  computations  are 
needed.  Secondly,  set  membership  classification  for  non-intersected  entities  is  not 
necessary  because  the  size  of  Ei^  and  Ei^  and  their  members  are  not  changed.  Finally, 
subdivision  computation,  which  involves  set  membership  classification  and  connectivity 
construction,  is  not  needed  any  more.  Moreover,  persistent  naming  can  be  avoided  since 
no  topological  entity  is  recreated. 

In  most  engineering  design  application,  dimension  change  does  not  change 
topology  of  the  solid  model.  In  the  literature  (Shapiro,  V.  and  Vossler,  D.  L.,  1995; 
Raghothama,  S.  and  Shapiro,  V.,  1998),  the  parametric  solid  models  with  different 
dimensions  but  same  topology  are  classified  as  in  same  "parametric  family".  Based  on  the 
no-topology-change  assumption,  a  user  needs  to  recreate  the  solid  model  manually  if 
topology  change  is  expected.  The  updating  algorithm  presented  in  this  section  works 
efficiently  in  the  case  of  no-topology-change.  However,  the  no-topology-change 
assumption  is  a  limitation.  To  improve  the  designer's  flexibility,  an  extension  of  the 
above  updating  algorithm  is  presented  in  the  next  section  to  handle  topology  change. 
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Updating  Algorithm  Considering  Topology  Change 

Considering  the  case  of  topology  change,  the  recursive  updating  algorithm  for 
updating  Si  in  figure  5-1  can  be  expressed  as: 

1.  update  Su; 

2.  update  Siu; 

3.  reposition  Srj; 

4.  detect  topology  change; 

5.  if  there  is  no  topology  change 

a)  update  GCE,"' )  and  G(Er^); 

else 

b)  update  E,"'; 

c)  update  E,^  and  Ei^; 

d)  update  E.'"''; 

Figure  5-6:  Updating  algorithm  considering  topology  change 

The  first  three  steps  are  the  same  as  those  in  figure  5-5.  Step  4  is  for  detecting 
possible  topology  change.  Topology  change  can  be  detected  by  re-doing  intersection 
computation  to  check  if  any  member  or  size  of  E,""  changes.  Figure  5-7  shows  some 
examples  of  topology  change  when  changing  constraint  dimensions  of  Siu  in  Figure  5-2. 

From  figure  5-7  (a)  to  (b),  the  number  of  intersection  edges  in  Ei'"'  changes  from 
six  to  four.  Three  of  six  original  intersection  edges  do  not  exist  any  more  and  one  new 
intersection  edge  is  generated  and  added  to  Ei"".  From  figure  5-7  (a)  to  (c),  the  number  of 
intersection  edges  in  Ei*"'  does  not  change  and  is  still  six.  However,  the  members  of  Ei"* 
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have  changed.  Four  of  the  original  intersection  edges  in  the  set  Ei  do  not  exist  any  more 
while  four  new  intersection  edges  have  been  generated  and  added  to.  From  figure  5-7  (b) 
to  (d),  the  number  of  intersection  edges  in  Ei"*  does  not  change  and  is  still  four.  However, 
none  of  the  original  intersection  edges  exist  any  more  while  four  new  intersection  edges 
have  been  generated  and  added  to  Ei"".  In  this  example,  36  (6x6)  face  to  face  intersection 
computations  are  needed  to  detect  topology  change.  Therefore  step  4  is  the  same  as 
intersection  computation  required  during  Boolean  operation. 


(a) 
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Figure  5-7:  Examples  of  topology  change  (a)  case  1;  (b)  case  2;  (c)  case  3;  (d)  case  4. 


If  no  topology  change  is  found  in  step  4,  step  5(a)  will  be  executed.  Otherwise, 
step  5(b)  will  be  executed.  The  procedure  for  step  4,  5(a)  and  5(b)  are  described  below: 
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for  any  face  pair:  fjS  EfR,  and  fys  EfLi  (j  =  1,...,  size  of  EfR,  and  k  =  1,...,  size  of  Em) 
compute  the  intersection  of  this  face  pair; 
if  these  two  faces  intersect 

if  this  intersection  existed  before 

update  the  geometry  of  corresponding  intersection  edge  and  vertices; 

else 

create  new  intersection  entities  based  on  the  intersection  result 
and  add  them  to  E,"'; 

return  value  indicating  topology  change  happens; 

else 

if  this  intersection  existed  before 

delete  the  corresponding  intersection  entities  from  E"'; 
return  value  indicating  topology  change  happens; 

Figure  5-8;  Topology  change  detection 

During  this  procedure,  we  not  only  checked  if  topology  change  occured,  but  also 
updated  G(Ei'"')  when  there  was  no  topology  change  and  updated  Ei""*  when  there  was 
topology  change.  When  there  is  no  topology  change,  G(Ei^"'')  gets  automatically  updated 
after  updating  G(Ei"").  When  there  is  topology  change,  we  update  Hi'"*  instead  of  GCEi"*). 
Updating  Ei'"'  involves  three  steps:  updating  the  geometry  of  the  intersection  entities  that 
existed  before  and  still  exists  after  dimension  change,  deleting  intersection  entities  which 
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do  not  exist  any  more  from  Ei'"^  and  adding  new  intersection  entities  into  E-^.  All  deleted 
intersection  entities  are  temporarily  stored  as  a  set  Ei*^^'"'"'  and  all  new  intersection  entities 
are  temporarily  stored  as  a  sub-set  E;"'"'-'^  where  Ei"'''"'"'^  updated  Ei*^.  These  two  are 
needed  in  step  5(d). 

In  step  5(c),  we  need  to  re-classify  all  non-intersecting  entities  of  Su  and  Sri..  The 
purpose  of  the  reclassification  is  to  delete  topological  entities,  which  do  not  belong  to  the 
final  solid  any  more,  from  Ei^  and  Ei^,  and  to  create  and  add  new  topological  entities  into 
Ei'^  and  Ei^  that  share  geometry  with  the  corresponding  entities  in  Su  or  Sm.  For  example, 
in  Figure  5-7,  when  Sm  is  moved  down  with  respect  to  Su  from  (b)  to  (d),  the  face  of  Si 
that  shares  geometry  with  top  face  of  Sri  will  be  deleted  from  Ei'*  and  a  new  face  sharing 
geometry  with  the  bottom  face  of  Sri  will  be  created  and  added  to  Ei'^. 

The  subdivision  entity  set  Ei'"''  is  updated  in  step  5(d),  which  can  be  divided  into 
two  parts:  1)  updating  subdivision  edge  sets  Esr""''  and  Ecl"'';  2)  updating  subdivision 
face  sets  Em""''  and  EfL*"*'.  The  procedure  for  updating  subdivision  entity  is  similar  to  that 
in  Boolean  operation  and  requires  set  membership  classification.  The  difference  is  that 
subdivision  entities  that  existed  before  and  exist  after  dimension  changes  can  be  modified 
but  not  recreated  in  Ei'"''.  Similar  to  step  5(b)  and  step  5(c),  the  subdivision  entities  that 
do  not  belong  to  Si  any  more  will  be  deleted  from  Ei'"''  and  new  subdivision  entities  will 
be  added  to  Ei'"''. 

Let  e'"^  be  a  subdivision  edge  of  the  updated  solid  model  Si  whose  end  vertices 
are  vi  and  V2 ,  there  are  five  possible  cases  : 

1 .  vi  or  V2  G  Ei"'  -  Ei"""-'"'  and  the  other  g  Ei'^  u  Ei^; 
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2.  V,  or  V2  €  Er     and  the  other  e  E;''  kj  E,^ 

3.  both  vi  and  V2eEr-Ei~'; 

4.  both  vi  and  V2  g  E^^-'"'  ; 

5.  vi  or  V2  G       -  Er''"'  and  the  other  e  Ei"'""'"'; 

The  first  two  cases  can  be  easily  found  in  figure  5-7.  Assuming  the  soUd  is 
changed  from  figure  5-7  (a)  to  (b),  subdivision  edge  ei  is  case  1  and  e2  is  case  2.  The  last 
two  cases  are  shown  in  figure  5-9  (b)  and  (c).  Obviously,  if  any  intersection  vertex  e 
gdei-int^  the  subdivision  edge  using  this  vertex  as  one  end  does  not  exist  any  more  and 
should  be  deleted  fi-om  EeR^"''  or  EeL*"^. 

In  case  1  and  3,  nothing  needs  to  be  done  because  these  subdivision  edges  existed 
previously.  Examples  of  case  2,  4  and  5  are  shown  in  Figure  5-9.  In  these  cases,  new 
subdivision  edge  will  be  created  and  added  to  Esr'"''  or  Ebl'"''.  In  case  2  (figure  5-9(a)) 
and  5  (figure  5-9(c)),  we  need  to  use  the  new  intersection  vertex  to  do  set  membership 
classification  to  subdivide  the  edge  it  lies  on  and  find  another  existing  vertex  along  the 
classified  direction.  In  case  4  (figure  5-9(b)),  we  need  to  do  set  membership  classification 
two  fimes  to  create  the  new  subdivision  edge.  In  these  three  cases,  if  any  edge  e  Ei'^ ,  Ei^ 
,  EeR^'"''  or  EeL^"**  pass  through  both  end  vertices,  this  edge  should  be  deleted  from 
corresponding  subset. 
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Figure  5-9  Creation  of  new  subdivision  edge 
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Based  on  the  above  analysis,  we  do  not  need  to  use  any  vertexe  Ei""  -  Ej"^"'''"*  to 
do  classification.  Therefore,  for  updating  subdivision  edge  sets,  we  only  need  to  go 
through  all  intersection  vertices  in  Ei"*^'''"'  and  E i"*"^"'"'  to  do  following  procedure  for  each 
vertex: 

if  vertex  v  e  E,"*''  "* 

delete  all  subdivision  edges,  of  which  v  is  one  of  end  vertex ,  from 
EeR'""  or 

else  if  vertex  V6  E,""*"'"' 

use  methods  for  case  2,  4  and  5  to  subdivide  edges  which  belong  to  E^iu  u  EeLi 

and  pass  through  v  to  create  new  subdivision  edges  and  add  them  to 

EeR'"''  or  E.l'"''  ;  if  necessary,  delete  certain  edges  in  E,"^ ,  E,^ ,  Ecr''*  or  Ecl'"" 

Figure  5-10:  Procedure  for  updating  subdivision  edge  sets 

Updating  the  subdivision  face  sets  Efs""''  and  EfL'"''  consists  of  three  steps.  The 
first  step  is  to  delete  subdivision  faces,  which  do  not  belong  to  Si  any  more,  from  E{r"^ 
and  Enf"^.  The  procedure  to  determine  and  delete  such  a  subdivision  face  can  be  simply 
described  as: 

if  all  intersection  edges  in  this  face  e  Ei'*"'"'"' 
delete  this  face  from  Ek'"''  and  Ea,'"''; 

Figure  5-11:  Procedure  for  deleting  a  subdivision  face 

The  second  step  is  to  create  the  new  subdivision  faces  and  add  them  to  the  sets 
EfR^"''  or  EfL*"''.  The  third  step  is  to  modify  subdivision  faces,  which  still  exist  but  have 
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different  topology.  Let  f"^  be  a  subdivision  face  of  updated  Si  ,  the  following  possible 
cases  for  edges  forming  this  face  need  to  be  considered: 


1 .  at  least  one  edge  g  Ei""  -  Ei"*^"*"*  and  no  edge  belongs  to  Ei"~; 

2.  at  least  one  edge  g  Er""'"*  and  no  edge  belongs  to  Ej*"'  -  E^^"'"'; 

3.  at  least  one  edge  e  E;"'  -  Ei"'""*"'  and  at  least  one  edge  e  Ei"""'"'"*; 

Examples  of  all  these  three  cases  can  be  found  in  Figure  5-7.  When  solid  changes  from 
(a)  to  (b),  face  fl  and  f2  are  in  case  1,  face  D  is  in  case  2.  When  solid  changes  from  (b)  to 
(c),  fi  in  (c)  is  in  case  3. 

There  are  two  sub-cases  for  case  1  based  on  whether  the  number  of  edges  e  Ei"*  - 
gnew-int  ^hangcs  and  does  not  change.  Face  fl  and  f2  in  (b)  are  these  two  sub-cases 
respectively.  In  the  second  sub-case,  the  topology  of  this  subdivision  face  does  not 
change  and  nothing  needs  to  be  done.  In  the  first  sub-case,  we  need  to  modify  the  edge 
connectivity  of  this  subdivision  face  by  deleting  intersection  edges  g  Ei"*^'''"'  and 
subdivision  edges,  which  do  not  exist  anymore,  and  adding  new  subdivision  edges.  But 
the  geometry  of  this  face  does  not  change  and  no  new  face  is  created. 

In  case  2,  the  subdivision  face  is  a  brand  new  subdivision  face  and  needs  to  share 
geometry  with  the  corresponding  subdivided  face  in  Emi  or  EfLi.  For  this  face  edge 
connectivity  must  be  constructed  among  edges  from  Ei"^"''",  EeR*"''  or  EeL*"'' ,  and  Ei^  or 
Ei^.  In  this  case,  if  a  face  in  Ei"^  or  Ei^  shares  geometry  with  the  subdivided  face,  the  face 
must  be  deleted  from  Ei*^  or  Ei^. 
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In  case  3,  the  subdivision  face  still  exists  and  its  geometry  does  not  change. 
However,  we  need  to  reconstruct  edge  connectivity  among  edges  from  Ei'"',  EeR^"''  or 
EeL'"\and  E;^  orE^ 

The  algorithm  (shown  in  figure  5-6)  presented  in  this  section  can  handle  both 
cases  with  topology  changes  and  without  topology  change.  In  case  of  no  topology 
change,  the  efficiency  is  less  than  the  algorithm  in  figure  5-5  since  we  need  to  redo 
intersection  computation  to  determine  if  topology  change  happens.  But  this  updating 
algorithm  is  still  much  more  efficient  than  other  updating  algorithms  that  involve  redoing 
the  Boolean  operations  since  it  has  all  other  advantages  of  algorithm  in  figure  5-5. 

In  case  topology  change  happens,  this  algorithm  never  recreates  topology  entities 
that  existed  in  the  original  solid.  From  the  algorithm  described  above,  it  is  easy  to  see  that 
it  is  very  similar  to  redoing  Boolean  operation  except  that  it  does  not  recreate  topology 
entity.  In  the  worst  case,  the  efficiency  of  this  algorithm  is  similar  to  that  of  redoing 
Boolean  operations.  However,  it  still  has  the  advantage  of  not  using  persistent  naming  for 
topology  entities.  The  efficiency  of  this  algorithm  is  obvious  since  in  most  engineering 
design  modification  there  is  no  topology  change. 

Localization  of  the  Geometric  Modification 

The  algorithm  in  figure  5-6  can  be  further  improved  by  limiting  the  geometric 
modification  to  be  as  local  as  possible.  In  the  updating  algorithm,  if  any  of  Su  or  Sri  is 
not  modified,  one  or  both  of  step  for  updating  Su  and  step  for  updating  Sr]  can  be 
skipped.  In  figure  3-1,  for  example,  if  only  dimensions  of  P2  are  changed,  B-rep  stored  in 
leaf  nodes  PI  and  P3  and  internal  node  3  (G3)  are  not  affected.  Therefore,  when  the 


S8 

updating  algorithm  is  applied  to  nodes  1,  2  and  4,  we  can  skip  the  step  for  updating  PI, 
P3  and  G3  respectively.  Moreover,  step  3  for  repositioning  Sri  can  also  be  skipped  if  the 
constraint  value  and  geometric  entities  corresponding  to  all  topological  entities  involved 
in  complete  constraint  set  are  not  changed.  In  figure  3-1,  repositioning  in  nodes  2  and  4 
can  be  skipped  because  the  complete  constraint  sets  in  these  two  nodes  are  not  affected 
by  changing  dimensions  of  P2. 

There  are  many  different  ways  to  localize  the  geometric  modification  when 
implementing  the  updating  algorithm.  Two  possible  ways  are  given  below. 

For  skipping  steps  1,  2  for  updating  Su  or  Siu  in  certain  cases,  an  integer  flag 
called  modify_flag  is  stored  in  each  node  of  the  extended  CSG  tree.  Originally,  this  flag 
is  set  to  be  0  for  any  node,  which  means  the  geometry  represented  by  the  node  does  not 
need  to  be  updated.  If  any  dimension  of  a  primitive  is  modified,  the  modify_flags  of  the 
leaf  node  corresponding  to  this  primitive  and  internal  nodes  on  the  path  from  this  leaf 
node  to  root  node  will  be  set  to  be  1.  If  any  constraint  dimension  in  a  complete  constraint 
set  stored  in  an  internal  node  is  modified,  the  modify  flags  of  this  node  and  all  internal 
nodes  on  the  path  from  this  node  to  root  node  is  also  set  to  be  1.  The  modify_flag  of  one 
node  is  1  means  the  geometry  represented  by  this  node  needs  to  be  updated. 

To  avoid  repositioning  the  dependent  geometry  in  step  3  in  certain  cases,  an 
integer  flag  called  moved  flag  can  be  attached  to  each  geometric  entity  as  one  of  its 
attributes.  Originally,  moved  flags  of  all  geometric  entity  are  set  to  be  0,  which  means 
that  no  change  happened  to  any  geometric  entity.  A  moved  flag  can  be  set  to  be  1,  which 
means  change  happened  to  corresponding  geometric  entity.  Three  types  of  changes  are 
possible.  In  first  case,  when  primitive  geometry  is  modified  due  to  dimension  change,  the 
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moved  flags  corresponding  to  some  changed  geometric  entities  is  set  to  1.  For  example, 
in  figure  5-3,  moved_flag  of  line  of  edge  5  6  is  set  to  1  while  move_flag  of  line  of  edge 
1-2  is  still  0  after  primitive  modification.  In  second  case,  when  a  geometric  entity  is 
translated  or  rotated  by  a  transformation  matrix  in  updating  process,  the  moved  flag  of 
this  entity  is  set  to  1.  In  last  case,  if  the  geometric  entity  of  an  intersection  entity  is 
updated  by  re-computing  the  intersection,  the  corresponding  moved  flag  is  set  to  1.  A 
complete  constraint  set  of  an  internal  node  includes  topological  entities  involved  and  each 
of  these  topological  entities  has  its  geometric  entity.  If  moved  flag  of  any  one  of  these 
geometric  entities  is  1,  repositioning  of  the  node's  right  child  with  respect  to  left  child  is 
necessary.  Otherwise,  the  repositioning  (step  3)  can  be  skipped  when  updating  algorithm 
is  applied  to  this  node. 

At  the  end  of  updating  process  of  the  solid  model,  the  modify  flag  of  each  node 
in  extended  CSG  tree  and  moved  flag  of  all  geometric  entities  is  reset  to  0.  Considering 
localization  of  geometric  modification,  the  updating  algorithm  in  figure  5-6  can  be 
rewritten  as: 

1 .  if  modify_flag  of  Sl,  =  1,  update  Su, 

2.  if  modify_flag  of  Sr,  =  1,  update  Siu; 

3.  if  necessary,  reposition  Sri; 

4.  detect  topology  change; 

5.  if  no  topology 

a)  update  GCE,""' )  and  G(Er\ 

else 

b)  update  Ei"' ; 

c)  update  E,"^  and  E,^; 

d)  update  E,'""; 
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Figure  5-12:  Revised  updating  algorithm 

Step  4  and  5  in  the  algorithm  can  be  thought  of  as  a  procedure  to  modify  the 
corresponding  Boolean  operation.  Assuming  this  procedure  can  be  represented  by  M(- ) 
or  M(yj*)  which  correspond  to  -*  and  u'  respectively,  the  updating  process  for  solid 
model  in  figure  3-1  when  dimension  of  primitive  P2  is  changed  can  be  represented  as  a 
list  of  operations  as  following: 

{  Mp2,  Pos(P2,  PI),  PI  M(-*)  P2,  Gl  M(u*)  P3,  G2  M(-*)  03 }  (5.26) 
where  Mp2  represents  the  procedure  of  modifying  primitive  P2.  Since  no  new 
topological  entity  is  created,  Mp2  is  more  efficient  than  Cp2  in  (5.2).  As  discussed  before, 
even  in  worst  case  (topology  change  happened)  efficiency  of  step  4  and  5  is  still  better 
than  redoing  Boolean  operation  since  it  eliminates  the  need  for  persistent  naming.  In  most 
engineering  design  modification  (no-topology-change)  the  efficiency  of  step  3  and  4  is 
much  better  than  redoing  Boolean  operation.  Therefore,  operations  M(-*)  and  M(u')  are 
more  efficient  than  -*  and  u*.  Comparing  (5.26)  with  (5.2),  there  is  no  doubt  that  the 
algorithm  in  figure  5-12  is  much  more  efficient  that  "rollback"  algorithm. 


CHAPTER  6 

IMPLEMENTATION  OF  A  DIMENSION-DRIVEN  SOLID  MODELING  SYSTEM 

Overview 

To  verify  the  proposed  algorithms  for  constraint  imposition  and  updating,  an 
experimental  Feature-Based  solid  Modeling  system  (FBM),  which  can  perform 
dimension-driven  solid  modeling,  is  implemented  using  object-oriented  programming 
language  C++  (Deitel,  H.  M.  and  Deitel,  P.  J.,  1994)  on  IRIX  6.2  operation  system  on 
SGI  workstation.  The  system  uses  ACIS  (ACIS,  1996)  as  the  underlying  geometric 
modeler  to  create  and  display  the  geometry  and  uses  Motif  window  manager  to  create 
user  interface.  A  stand-alone  constraint  imposition  package  CONSPACK  was 
implemented  without  using  any  data  structure  from  ACIS.  The  FBM  system  can  only  run 
in  SGI  workstation  that  has  ACIS  and  Motif  available  while  CONSPACK  is  written  in 
pure  standard  C++  and  thus  can  be  used  on  different  platforms  for  other  applications. 
FBM  and  CONSPACK  together  consist  of  about  21000  lines  of  C++  source  code.  In  this 
chapter,  some  implementation  details  are  introduced  and  discussed. 


System  Structure 

The  dimension  driven  solid  modeling  system  architecture  can  be  illustrated  as  in 
figure  6-1.  The  FBM  system  consists  of  a  User  Interface  and  a  Dimension-driven 
Modeler.  The  User  Interface  creates  a  graphical  window,  menu  bars  and  dialogue  boxes 
to  let  user  issue  commands  to  Dimension-driven  Modeler.  The  appearance  of  the 
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graphical  user  interface  and  the  window,  menu  bars  will  be  introduced  in  chapter  7.  The 
Dimension-driven  Modeler  is  the  core  of  the  FBM  system  and  consists  of  many  C++ 
classes  and  some  utility  functions.  Responding  to  commands  from  the  User  Interface,  the 
Dimension-driven  Modeler  can  perform  operations  such  as  solid  model  creation, 
dimension  modification  and  solid  model  updating,  etc.  During  these  operations,  the 
Dimension-driven  Modeler  needs  to  communicate  with  ACIS  and  CONSPACK  and 
display  the  resuhs  in  the  graphical  window  of  the  User  Interface.  In  the  following 
sections,  we  will  firstly  introduce  solid  representation  scheme  in  ACIS.  After  that,  we 
will  focus  our  discussion  on  the  Dimension-driven  Modeler  and  CONSPACK. 


User  Interface 


Dimension-driven 
Modeler 


FBM 


T 


CONSPACK 


ACIS  geometric 
modeler 


Figure  6-1:  Structure  of  Dimension-driven  solid  modeling  system 


Solid  Representation  in  ACIS 

ACIS  represents  a  solid  using  B-rep  and  its  data  structure  is  very  similar  to  half- 
edge  data  structure  introduced  in  chapter  2.  The  entities  used  to  construct  a  solid  model  in 
ACIS  can  be  classified  into  topological  entities  and  geometric  entities.  Figure  6-2  shows 
the  topological  entity  hierarchy  in  ACIS  representation  and  the  corresponding  geometric 
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entity  types.  This  figure  is  modified  from  ACIS  manual  (ACIS,  1996)  by  deleting  some 
topological  entity  types  not  used  in  our  system,  such  as  SUB  SHELL  and  WIRE. 


BODY 
LUMP 

m 

SHELL 


I 


FACE 

SURFACE 

LOOP 


COEDGE 

HZ 


EDGE 


CURVE 


VERTEX 


^  APOINT 


Figure  6-2:  Topological  and  geometric  entities  in  ACIS 


BODY  is  the  highest  level  topological  entity  and  all  solid  models  including  2D  or 
3D  geometry  are  represented  as  BODY  in  ACIS.  A  BODY  is  composed  of  LUMPs,  each 
of  which  is  a  ID,  2D  or  3D  set  of  points  in  space  that  is  disjoint  from  all  other  LUMPs.  A 
LUMP  is  composed  of  SHELLs,  each  of  which  is  a  set  of  connected  FACEs  that  are  used 
to  bound  the  outside  of  a  solid  or  an  internal  void  (hollow).  A  FACE  is  a  surface  bounded 
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by  a  LOOP  which  is  a  list  of  COEDGEs  containing  EDGEs  in  a  certain  direction.  A 
EDGE  is  a  curve  bounded  by  VERTEXs.  Among  these  topological  entities,  only  FACE, 
EDGE  and  VERTEX  are  entity  types  in  general  B-rep  and  have  their  own  geometric 
entities,  which  are  SURFACE,  CURVE  and  APOINT  respectively.  All  the  above 
topological  and  geometric  entity  types  are  defined  using  corresponding  classes  in  ACIS 
that  are  derived  from  a  base  class  named  ENTITY. 

Feature  Deflnition 

A  feature  defined  in  FBM  is  a  high  level  entity  which  contains  information  of 
geometric  definition,  attach  operator  and  parent-child  relation  with  respect  to  other 
features.  An  object  modeling  diagram  is  shown  in  figure  6-3  to  illustrate  the  FEATURE 
class  definition. 

An  attach  operator  of  a  feature  is  a  class  containing  a  set  of  position  constraints, 
which  can  be  evaluated  by  using  CONSPACK  to  position  the  feature  with  respect  to  a 
prior  feature  it  will  be  attached  to.  The  attach  operator  also  contains  methods  for  attach 
operation  and  modifying  attach  operation.  Attach  operation  creates  new  B-rep  of  the  prior 
feature,  which  we  refer  to  as  the  parent  feature,  by  attaching  geometry  of  the  feature  onto 
the  geometry  of  its  parent  feature.  Pointers  of  all  child  features  of  the  feature  are  stored  as 
FEATURE  list  in  the  feature  class  and  so  is  the  pointer  to  its  parent  feature.  Therefore,  a 
feature  has  pointers  pointing  to  its  parent  feature  and  to  all  its  child  features. 
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FEATURE 


Geometry  Creation  function 
Geometrymodify  function 
regenerate  function 


0 


Attach 

Operator 

I 


parent 


Independent 
Variables 


FEATURE  3D 


Profile 


child  features 


FEATURE  list 


Figure  6-3:  Object  modeling  diagram  of  feature  definition 


The  primitive  geometry  of  a  feature  is  defined  by  a  particular  Geometry  creation 
function,  which  takes  independent  variables  (dimensions)  of  the  primitive  as  input  and 
stores  the  output  B-rep  as  the  first  BODY  in  the  BODY_list  of  the  feature.  Independent 
variables  and  Geometry  creation  function  of  a  feature  are  defined  when  initializing  the 
object  of  FEATURE  class.  For  example,  a  feature  with  rectangle  geometry  can  be 
initialized  by  given  lengths  of  two  sides  and  a  rectangle  creation  function,  which  uses  the 
side  lengths  to  create  the  B-rep  of  the  rectangle.  The  feature  with  a  2D  primitive 
geometry  is  called  2D  feature  and  can  be  directly  initialized  from  the  FEATURE  class. 
Most  commercial  systems  allow  the  user  to  sketch  2D  profile  and  express  various 
geometric  constraints  and  dimensions  so  that  the  profile  can  be  updated  using  variational 
geometry  techniques.  In  the  implementation  of  FBM,  we  always  use  feature  attach 
operation  to  combine  simple  2D  profiles  to  create  2D  feature  with  complex  shape.  In  this 
way,  variational  geometry  method  can  be  avoided,  Geometry  modify  function  modifies 
the  geometry  of  the  feature's  primitive  in  response  to  value  changes  of  independent 
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variables.  Regeneration  function  of  a  feature  is  an  implementation  of  the  updating 
algorithm  introduced  in  chapter  5  and  will  be  further  discussed  late. 

A  class  FEATURESD  is  defined  as  a  derived  class  of  FEATURE  and  inherits  all 
protected  and  public  members  from  FEATURE  class.  Since  3D  primitive  geometry  is 
created  by  extruding,  sweeping  or  revolving  a  two-dimensional  profile,  there  is  an 
additional  data  member  called  Profile  Feature,  which  is  a  2D  feature  representing  the  2D 
profile,  in  the  definition  of  FEATURE_3D  shown  in  figure  6-3.  Independent  variables  of 
3D  feature  include  necessary  dimension(s),  such  as  extrusion  distance  or  revolving  angle. 
Geometry  creation  function  in  3D  feature  takes  Profile  Feature  and  independent 
variables  as  its  input  parameters. 

A  BODY  list  is  maintained  in  feature  definition  that  includes  a  list  of  B-reps  of 
the  feature  in  different  stage  of  model  construction  process.  Inifially,  this  list  only  has  one 
BODY,  which  is  the  B-rep  of  primitive  geometry  of  this  feature.  When  a  child  feature  is 
attached  to  this  feature,  the  child  feature  is  added  to  the  FEATURE  list  of  this  feature 
while  the  resuhing  B-rep  of  the  attach  operation  is  added  to  BODY  list.  There  is  a 
mapping  between  FEATURE  list  and  BODY  list  in  a  feature,  as  shown  in  figure  6-4. 


BODY] 

B0DY2 

BODY/ 

BODYn 

primitive 

I 

child  1 

child  /■-; 

child  n-1 

Figure  6-4:  Mapping  between  BODY  list  and  FEATURE  list 
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The  first  BODY  on  BODY  list  is  the  primitive  geometry  of  the  feature  and  the 
last  one  on  the  Ust  is  the  current  BODY  of  the  feature.  BODYJist  has  a  mapping  with 
child  feature  Hst:  /-th  BODY  on  BODY  list  is  the  result  of  attaching  the  current  BODY 
of  f/-y>th  child  feature  to  (i-l)-th  BODY  on  BODYJist. 


Parent-child  Feature  Tree 

In  chapter  3,  we  described  using  extended  CSG  tree  to  represent  design  intent. 
The  extended  CSG  tree  is  implemented  as  Parent-child  Feature  Tree  in  FBM  based  on  the 
feature  definition.  Using  F;  to  represent  the  feature  with  primitive  Pi,  the  Parent-child 
Feature  Tree  of  the  example  in  Figure  3-1  is  shown  in  Figure  6-5. 

Fl 


F2 


F3 


F4 


F5 


Figure  6-5:  Parent-child  Feature  Tree 


It  is  easy  to  see  that  Parent-child  Feature  Tree  is  a  complete  implementation  of  the 
extended  CSG  tree.  The  intermediate  geometry  stored  in  an  internal  node  of  the  extended 
CSG  tree  is  a  BODY  stored  in  the  BODY_list  of  a  feature,  which  is  the  one 
corresponding  to  the  left  most  leaf  node  of  the  sub-tree  rooted  at  this  internal  node.  The 
design  intent  represented  by  the  extended  CSG  tree  is  composed  of  three  parts:  primitive 
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creation  function  in  each  leaf  node,  complete  constraint  set  and  attach  operator  in  each 
internal  node.  The  primitive  creation  function  is  the  Geometry  creation  function  in  the 
feature.  The  complete  constraint  set  and  attach  operator  in  an  internal  node  of  extended 
CSG  tree  are  combined  and  implemented  as  attach  operator  of  a  feature.  If  the  right  child 
of  this  internal  node  is  a  leaf  node,  the  feature  is  the  one  corresponding  to  this  leaf  node. 
Otherwise,  the  feature  is  the  left  most  leaf  node  of  the  sub-tree  rooted  at  this  internal 
node. 

Functions  for  Creating  and  modifying  Primitives 
Primitive  Creation 

In  feature  definition,  Geometry  creation  function  is  actually  a  pointer  to  one  of 
functions  for  creating  primitives,  which  are  pre-defined  functions  in  the  FBM  system.  We 
have  implemented  Sanctions  for  creating  the  2D  primitives  shown  in  figure  6-6  and 
function  for  extruding  2D  profile  to  create  3D  primitive.  Functions  for  creating  more  2D 
primitives  and  function  for  creating  revolved  3D  primitive  can  be  easily  added  to  the 
system  in  the  fliture. 

The  independent  variables  of  each  primitive  are  labeled  in  the  figure.  A  square 
needs  only  one  parameter,  which  is  the  side  length,  to  define  its  geometry  while  a 
rectangle  needs  two  independent  variables,  which  are  the  lengths  of  its  two  sides.  Three 
triangle  creation  functions  are  defined  in  the  system.  Three  sides  of  triangle  l  have  same 
length,  so  it  only  needs  one  parameter  to  define  its  shape  while  triangle_2  needs  three 
parameters.  Triangle_3  has  a  constant  90°  angle  and  thus  needs  only  two  side  lengths  as 
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independent  variables.  Circle  has  one  independent  variable,  which  is  its  radius.  Both 
arc  l  and  arc_2  need  radius  and  an  angle  to  define  their  geometry. 


b 


square  rectangle         triangle_l         triangle_2  triangle_3 


circle  arc_l  arc  2 


Figure  6-6:  2D  Primitives 

Details  of  functions  for  creating  above  2D  primitives  are  trivial  since  only  2D 
geometry  computation  and  basic  ACIS  geometry  construction  operations  are  involved. 
However,  it  is  worth  mentioning  that  each  primitive  has  an  implied  local  coordinate 
system.  One  example  of  this  local  coordinate  system  is  shown  for  triangle_2  in  figure  6- 
6.  For  primitives  in  first  row  of  figure  6-6,  the  origin  of  the  local  coordinate  is  always 
located  at  left-bottom  vertex  of  the  primitive  and  the  x  axis  is  attached  to  the  bottom  edge 
of  the  primitive  along  the  direction  of  corresponding  COEDGE.  For  primitives  in  second 
row  of  figure  6-6,  the  origin  is  located  at  center  of  the  circle  or  arc.  The  x  axis  is  parallel 
to  that  of  the  global  coordinate  system  for  the  circle,  attached  to  right  side  edge  for  arc  l 
and  parallel  to  the  only  straight  edge  for  arc_2  respectively.  The  local  coordinate  system 
for  a  2D  primitive  is  important  in  the  primitive  modification  process. 
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End  face 


End  edge 


End  vertex 


•Extrusion  edge 


Extrusion  face 


copy  of  2D  profile  face 


Figure  6-7:  Example  of  extrusion 


A  function  called  create  extrusion  is  implemented  for  extruding  a  2D  profile  to 
create  a  3D  primitive.  The  extrusion  direction  is  the  normal  vector  of  the  planar  face  and 
the  distance  is  an  input  parameter  of  this  function.  This  function  firstly  reuses  the 
geometry  of  the  2D  profile  face  (including  the  geometry  of  edges  and  vertices  of  this 
face)  to  create  a  copy  of  the  profile  face.  Figure  6-7  shows  an  example,  which  is 
extrusion  of  a  square.  The  resulting  block  contains  the  profile  face  copy  with  the  edges 
and  vertices  on  it,  4  extrusion  edges,  4  extrusion  faces,  4  end  vertices,  4  end  edges  and  1 
end  face.  For  brevity,  we  call  the  edges  and  vertices  on  profile  face  copy  as  profile  edges 
and  profile  vertices  respectively. 

The  posifion  of  an  end  vertex  can  be  found  by: 


where  Ve  and  Vp  represent  the  positions  of  the  end  vertex  and  its  corresponding  profile 
vertex;  d  is  the  extrusion  distance  and  V  is  the  extrusion  direction  vector.  An  extrusion 
edge  is  a  straight  edge,  which  can  be  easily  created  using  corresponding  Ve  and  Vp.  An 
end  edge  is  an  offset  of  the  corresponding  profile  edge.  Each  extrusion  face  is  a  four-edge 


Ve- Vp  +  C/V 


(6-1) 
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face  bounded  by  the  corresponding  profile  edge,  end  edge  and  two  extrusion  edges 
corresponding  to  two  vertices  of  the  profile  edge.  The  geometry  of  the  extrusion  face  is 
determined  by  the  geometry  of  the  corresponding  profile  edge.  The  face  is  a  planar  face  if 
the  profile  edge  is  a  straight  edge  and  is  a  cylindrical  face  when  the  profile  edge  is  an  arc. 
The  end  face  is  a  planar  face  with  a  normal  same  as  the  extrusion  direction  and  is 
bounded  by  all  end  edges. 

Primitive  Modification 

The  need  for  an  algorithm  for  modifying  a  2D  primitive  has  been  described  in 

chapter  5.  For  each  type  of  2D  primitive,  we  implemented  a  separate  primitive 
modification  function,  which  takes  a  pointer  of  FEATURE  as  input.  The  modified 
independent  variables  and  the  B-rep  of  the  primitive  in  current  position  can  be  obtained 
using  this  pointer.  The  primitive  modification  function  changes  the  geometry  of  the  2D 
primitive  of  the  given  feature,  while  the  orientation  and  position  of  the  local  coordinate 
system  is  kept  unchanged 

A  function  called  modify  extrusion  is  implemented  to  modify  the  3D  extrusion 
primitive.  This  function  takes  a  pointer  of  FEATURE_3D,  pt  3D,  as  input.  The 
implementation  of  this  function  can  be  illustrated  as  flowchart  showing  in  figure  6-8.  A 
integer  flag,  modify  flag  as  introduced  in  last  section  of  chapter  5,  is  used  to  specify  if 
dimension  change  occurred  in  the  2D  profile. 
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using  pt  3D  to  get  profile 
feature  and  extrusion  distance, 
get  the  modified  flag  of 
profile  feature 


profile  feature 
^regeneraUon  fimction 
to  updating  its  geometry 

\ 

re-compute  the  geometry 
of  end  vertices,  extrusion 
edges,  end  edges  and 
extrusion  faces 

no 


no 


move  end  vertices  a  distance 
new  distance  -  old  distance 
along  ex'trusion  direction 


Figure  6-8:  Flowchart  of  fUnction  modify  extrusion 
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Definition  of  Attach  Operator 

Attach  operator  is  implemented  as  a  class  called  Attach  operator.  The  definition 
of  the  Attach  operator  is  illustrated  in  figure  6-9  (a)  using  an  object  modeling  diagram. 
The  Attach_operator  is  just  an  abstract  class,  which  is  used  to  derive  more  specific  type 
of  attach  operators  including  Protrusion_2D,  Cut_2D,  Hole_2D,  Fillet_2D,  Chamfer_2D 
and  Boolean  union  for  3D.  The  inheritance  hierarchy  is  shown  in  figure  6-9  (b). 

The  abstract  class  Attach  operator  contains  a  pointer  to  the  feature  it  belongs  to,  a 
pointer  to  an  object  of  class  ConstraintSet,  which  is  defined  in  CONSPACK,  and  a  list  of 
pointers  to  objects  of  class  RELATION.  The  feature  is  the  child  feature  involved  in  this 
attach  operation  and  its  parent  feature  can  be  obtained  through  the  child  feature.  The 
object  of  ConstraintSet  is  firstly  initialized  when  the  Attach_operator  is  initialized.  In 
each  updating  process,  if  the  child  feature  needs  to  be  repositioned  with  respect  to  its 
parent,  the  existing  ConstraintSet  object  is  deleted  and  a  new  object  is  initialized  again. 
The  ConstraintSet  object  is  the  tool  for  evaluating  all  user-defined  position  constraints  to 
position  the  child  feature.  RELATION  is  an  abstract  class  for  user-defined  position 
constraint  and  can  be  used  to  derive  specific  position  constraint  classes.  RELATION 
class  has  a  member  function  called  evaluate  relation  for  evaluating  its  position  constraint 
to  get  the  corresponding  transformation  matrix.  The  function  evaluate  constraints  in 
Attach  operator  is  used  to  fix  all  degrees  of  freedom  of  the  child  feature  with  respect  to 
its  parent.  The  implementation  of  this  function  is  quite  simple.  The  function  loops 
through  all  user-defined  position  constraints  on  RELATION_list  and  calls  its 
evaluation  relation  function  to  impose  the  constraint. 


104 


FEATURE 

Attach  operator 

evaluateconstraints  function 
virtual  attach  function 

ConstraintSet 


RELATION  list 


(a) 


Attachoperator 


Protrusion  2D 


Cut  2D 


Hole  2D 


Fillet  2D 


Chamfer  2D 


Boolean  union 


(b) 


Figure  6-9:  Object  modeling  diagram  and  class  hierarchy  of  class  Attach  operator:  (a) 
definition  of  Attach_operator  class;  (b)  derived  classes  of  Attach_operator 


Generally,  the  procedure  for  Boolean  attach  operation  can  also  be  used  for  2D 
attach  operation.  However,  using  some  specially  defined  attach  operations,  such  as  those 
list  in  figure  6-9  (b),  in  2D  case  can  provide  higher  efficiency  than  using  general  Boolean 
operation.  In  fact,  these  2D  attach  operations  are  simplifications  of  general  Boolean 
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operation  in  3D  case.  Since  each  class  derived  from  Attach  operator  has  different  attach 
operation,  the  attach  Junction  in  base  class  is  a  virtual  fUnction  and  in  each  derived  class 
there  is  a  fianction  with  the  same  name  as  the  virtual  function. 

Definition  of  RELATION  Class 

User-defined  position  constraints  are  implemented  as  derived  classes  of  the 
abstract  class  RELATION.  The  definition  of  class  RELATION  can  be  described  in  figure 
6-10  using  object  modeling  diagram.  RELATION  has  a  pointer  to  child  feature.  When 
the  position  constraint  corresponding  to  this  RELATION  is  evaluated,  a  transformation 
matrix  is  obtained  and  is  applied  to  the  current  BODY  of  the  child  feature.  Other  data 
stored  in  RELATION  class  include  constraint  type,  parent  entity,  child  entity  and 
constraint  value.  Constraint  type  is  an  integer  specifying  the  type  of  user-defined  position 
constraint  represented  by  the  object  of  this  class.  Parent  entity  and  child  entity  are  two 
topological  entities  involved  in  this  position  constraint,  which  are  from  current  parent 
BODY  and  child  BODY  respectively.  Constraint  value  is  an  angle,  distance  or  just  a  flag, 
which  is  needed  for  defining  the  position  constraint.  Function  evaluate  relation  is  virtual 
fianction  and  there  is  a  corresponding  function  with  the  same  name  in  each  derived  class 
of  RELATION. 

We  have  implemented  18  user-defined  position  constraints,  each  of  which  is 
derived  from  the  RELATION  class.  Figure  6-11  shows  the  names  of  these  derived  classes 
and  fiill  names  of  the  corresponding  position  constraints.  The  evaluation  relation  in 
each  of  these  derived  classes  converts  the  user-defined  position  constraint  into  a 
constraint  type  defined  in  chapter  4  and  then  uses  ConstraintSet  to  impose  this  constraint. 
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The  converting  process  is  quite  straightforward  for  user-defined  3D  constraint,  with  the 
assumption  that  all  faces  involved  in  these  position  constraints  are  planar  faces. 


FEATURE 

RELATION 

constraint  type 

Parent  entity 

Child  entity 

constraint  value 

virtual  evaluate  relation 

function 

Figure  6-10:  Object  modeling  diagram  for  class  RELATION 

Some  2D  constraints  are  also  defined  to  position  one  2D  profile  with  respect  to 
another.  Two  2D  angular  constraints  EEP_2D  and  EEA_2D  can  be  directly  converted  to 
PV  and  WA  constraints  in  chapter  4  respectively.  Distance  constraints  EEC_2D, 
VVC_2D,  CEN_2D  and  VVD_2D  can  be  directly  converted  to  VLC,  WC,  WC  and 
WD  constraints  defined  in  chapter  4  respectively.  Other  three  distance  constraints 
EED_2D,  EVD  2D  and  ECD_2D  can  be  converted  to  VPD  constraint.  However,  for  any 
edge  in  these  three  constraints,  a  plane  passing  through  it  and  perpendicular  to  the  2D 
profile  will  be  constructed  and  used  to  replace  the  edge  in  the  constraints. 
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WD  2D 


VVC  2D 


EED  2D 


EEA  2D 
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CEN  2D 
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edge  to  edge  distance  in  2D 


edge  to  edge  angle  in  2D 


edge  to  vertex  distance  in  2D 


edge  to  center  distance  in  2D 


center  to  center  coincident  in  2D 


face  to  face  parallel  in  3D 


face  to  face  angle  in  3D 


face  to  face  merge  (coincident)  in 


face  to  face  distance  in  3D 


face  to  vertex  distance  in  3D 


vertex  to  vertex  coincident  in  3D 


vertex  to  edge  coincident  in  3D 


vertex  to  vertex  distance  in  3D 


vertex  to  edge  distance  in  3D 


Figure 


6-11:  Inheritance  hierarchy  of  RELATION  classes 
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Definitions  of  Protrusion_2D,  Cut_2D  and  Hole_2D  Classes 

Protmsion_2D  class  inherits  all  public  and  protected  data  of  Attach  operator  and 
has  its  own  attach  function.  The  input  parameters  of  this  attach  function  are  two  2D 
geometries  of  parent  and  child  features.  Before  performing  this  attach  operation,  one 
child  feature's  edge  should  overlap  one  parent  feature's  edge.  These  two  edges  are  called 
merged  edges  in  Protrusion_2D  attach  operation.  Figure  6-12  shows  an  example  of  this 
kind  of  attach  operation.  The  situations  before  positioning,  before  attaching  and  after 
attach  operation  are  shown  in  (a),  (b)  and  (c)  respectively. 


(a)  (b)  (c) 

Figure  6-12:  Example  of  Protrusion_2D  attach  operation  (a)before  positioning;  (b)before 

attaching;  (c)  after  attach  operation 


Cut_2D  class  is  very  similar  to  Protrusion_2D  and  also  inherits  data  from 
Attach  operator  while  it  has  different  attach  fiinction.  Figure  6-13  shows  the  procedure 
of  a  Cut_2D  attach  operation. 


(a)  (b)  (c) 

Figure  6-13:  Example  of  Cut_2D  attach  operation  (a)before  positioning;  (b)before 

attaching;  (c)  after  attach  operation 
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Hole_2D  class  inherits  data  from  the  Attach  operator  and  also  has  its  own  attach 
function.  The  2D  geometry  of  the  child  feature  must  be  a  single  loop  profile.  Before 
performing  Hole_2D  attach  operation,  the  child  feature  should  be  positioned  totally 
inside  the  profile  of  parent  feature.  It  is  assumed  that  there  is  no  intersection  between  the 
profiles  of  child  feature  and  parent  feature.  Figure  6-14  shows  the  example  of  Hole_2D 
attach  operation. 


(a)  (b)  (c) 

Figure  6-14:  Example  of  Hole_2D  attach  operation  (a)before  positioning;  (b)before 

attaching;  (c)  after  attach  operation 

It  is  notice  that  no  new  vertex  is  created  during  above-mentioned  2D  attach 
operations.  The  intersection  entity  set  is  empty.  All  geometric  entities  in  the  resulting  2D 
geometry  are  entities  reused  from  two  2D  geometries.  Therefore,  during  the  updating 
process,  step  4  and  5  of  algorithms  presented  in  figure  5-6  and  5-12  can  be  omitted.  This 
is  a  major  advantage  for  defining  2D  attach  operators. 

Definitions  of  Filiet_2D  and  Chanifer_2D  Classes 

Fillet  2D  and  Chamfer_2D  classes  have  some  characteristics,  which  are  different 
from  other  2D  attach  operators.  The  child  feature  involved  in  this  attach  operation  has 
independent  variables  for  defining  the  size  of  the  fillet  or  chamfer,  but  it  does  not  have  a 
primitive  geometry.  Therefore,  in  the  attach  function  the  input  for  child  feature's  BODY 
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is  always  NULL.  The  object  modeling  diagrams  for  these  two  classes  are  shown  in  figure 
6-15  (a)  and  (b). 

There  are  three  new  data  members,  edge!,  edge2  and  fillet  or  chamfer  edge,  in 
these  two  classes.  Fillet  or  chamfer  is  buih  on  the  first  two  edges.  Since  the  child  feature 
representing  fillet  or  chamfer  does  not  have  its  geometry,  there  is  no  need  to  define  any 
position  constraints  for  these  attach  operators.  However,  for  consistence  purpose,  we 
store  fillet  edge  and  chamfer  edge  in  these  two  classes  and  then  define  two  dummy 
position  constraints  classes:  FLT_2D  for  fillet  and  CHF_2D  for  chamfer.  Instead  of 
computing  transformation  matrix,  the  evaluate  relation  in  these  two  classes  updates  the 
geometries  of  the  fillet  and  chamfer  edge.  The  attach  operations  for  creating  fillet  and 
chamfer  are  shown  in  figure  6-16. 


Fillet  2D 


edgel 
edge2 
Fillet  edge 
attach  function 

0 


Public  and  protected 
data  member  of 
Attachoperator 


Chamfer  2D 


edgel 
edge2 

chamfer  edge 
attach  fianction 


Public  and  protected 
data  member  of 
Attachoperator 


(a) 


(b) 


Figure  6-15:  Object  modeling  diagram  for  Fillet_2D  and  Chamfer_2D  (a)  Fillet_2D 

class;  (b)  Chamfer_2D. 


Ill 


dl/\d2 


one  independent 
variable:  radius 


two  independent 
variables:  dl  and  d2 


(a) 


(b) 


Figure  6-16:  Examples  of  Fillet_2D  and  Chamfer_2D  attach  operations  (a)  Fillet_2D 
attach  operation;  (b)  Chamfer_2D  attach  operation. 

Deflnition  of  Boolean_union  Class 

The  definition  of  Boolean  union  class  can  be  described  using  the  object  modeling 
diagram  shown  in  figure  6-17.  Boolean  union  class  also  inherits  all  public  and  protected 
member  of  class  Attach  operator.  The  attach  function  of  this  class  takes  the  current 
BODYs  of  the  parent  and  the  child  feature  as  inputs  and  returns  a  new  BODY,  which  is 
the  resuh  of  the  Boolean  union  operation  applied  on  these  two  BODYs.  The  function 
modify  bin  union  is  used  to  modify  the  resultant  BODY  when  dimensions  of  the  two 
BODYs  involved  in  the  attach  operation  or  their  relative  positions  are  changed.  This 
function  is  used  in  the  updating  process  to  perform  steps  4  and  5  in  the  updating 
algorithm  presented  in  figure  5-6  and  figure  5-12.  Assuming  the  child  feature  is  the  /-th 
child  on  the  child  FEATURE  list,  the  inputs  of  this  function  are  the  current  BODY  of  the 
child  feature,  /-th  and  (i+l)-Xh.  BODY  on  parent  feature's  BODYJist.  The  resuh  of  this 
function  is  that  the  (i+l)-X\y  BODY  on  parent  feature's  BODYJist  is  modified. 
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Booleanunion 

attach  function 

modify  bln  union  function 

5 


i 

i 

Public  and  protected  data 
member  of  Attach  operator 

int s  edge  li St ;  source  face  l i st ; 
ints_vtx_list 

Figure  6-17:  Object  modeling  diagram  of  class  Boolean  union 


All  intersection  edges  and  intersection  vertices,  which  are  generated  in  Boolean 
union  operation,  are  stored  in  Boolean  union  class  as  ints  edge  list  and  ints  vtx  list.  A 
list  of  faces  are  stored  in  this  class  as  source  face  list  to  specify  intersected  face  pair  for 
each  intersection  edge.  The  length  of  source  face  list  is  two  times  the  length  of 
ints  edge  list.  There  exists  a  mapping  between  these  two  lists  as  shown  in  Figure  6-18 
that  /-th  edge  on  ints  edge  list  is  the  intersection  edge  of  (2i-J)-th  face  and  2;'-th  face  on 
source  face  list. 


ints_edge_list 

1 

n 

ir  IP 

source_face_list 

1 

2 

2/-1 

2i 

2n-l 

2n 

Figure  6-18:  Mapping  between  ints  edge  list  and  source  face  list 


The  process  of  Boolean  union  attach  operation  can  be  illustrated  by  flowchart  in 
Figure  6-19.  Function  modify_bln_union  is  the  implementation  of  step  4  and  5  of  the 


113 

updating  algorithm  described  in  figure  5-6  and  5-12.  The  implementation  can  be 
illustrated  by  flowchart  in  figure  6-20.  As  mentioned  before,  there  are  three  BODYs 
involved  in  this  function.  In  following  flowchart,  we  use  A  and  B  to  represent  parent  and 
child  BODYs  involved  in  Boolean  union  operation.  The  BODY  resulting  from  the 
Boolean  operation  is  represented  as  C. 


fori=l  to  number  of 
faces  in  parent  BODY 


for  j=l  to  number  of 
faces  in  child  BODY 


compute  intersection  between  face  i  and  j, 
if  intersection  exists,  add  intersection  edge, 
vertices  to  ints_edge_list  and  ints_vtx_list. 
add  face  i,  j  to  source  face  list 


classify  all  unintersected  faces  and  add 
qualified  faces  to  resulted  BODY 


I 


for  i=  1  to  length  of 
ints_edg_list 


get  face  1,  face2  corresponding 
to  edge  i  from  source  face  list 


start  from  intersection  edge  i  to  create  subdivision 
faces  of  facel  and  face2.  subdivision  edges  are 
created  in  this  step 

add  created  subdivision  faces  to  resulted  BODY 


Figure  6-19:  Flowchart  of  Boolean  union  attach  function 
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topologychange  =0 


for  i=l  to  number  of 
faces  in  BODY^ 


forj=l  to  number  of 
faces  in  BODY  B 


compute  intersection  between  face  i  and  j 


add  new  intersection 
edge  and  vertices  to 
intsedgelist  and 
ints  vtx  list,  add 
pair(i  j)  to 

sourcefacelist  and 
mark  it  NEW, 
topology  change  =  1 


T 


mark  pair(ij)  ON, 
update  geometry  of 
corresponding 
intersection  edge 
and  vertices 


mark  pair  (ij)  OFF  from  source  face  list 
mark  corresponding  intersection  edge  OFF 
from  ints  edge  list;  topology_change=l 


Figure  6-20:  Flowchart  of  function  modify  bln  union 
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The  procedure  face jjrocessor  is  used  to  perform  step  5(c)  and  5(d)  of  the 
algorithms  in  figure  5-6  and  figure  5-12.  This  procedure  can  be  illustrated  by  flowchart  in 
figure  6-21.  In  this  procedure,  a  small  function  face  status  test  will  test  the  given  face 
based  on  the  information  stored  in  source  face  list  and  the  marks  generated  and  return  an 
integer  k.  There  are  following  cases  for  value  of  k. 

•  k=  -1 :  this  face  is  not  on  source  face  list 

•  k=  0:  this  face  was  originally  on  source  face  list,  but  is  not  anymore 

•  k=  1 :  this  face  was  not  on  source  face  list,  but  is  now 

•  k=  2:  this  face  was  an  intersected  face  and  is  still  one  without  any  new  intersection 

•  k=  -2:  this  face  was  an  intersected  face  and  is  still  one,  but  with  new  intersection 
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k=  facestatustest  (face) 

kl  =  classifyface  (face,  body) 

(if  the  face  is  a  parent  face,  the  body  is 

child  BODY  B:  if  the  face  is  a  child 

face,  the  body  is  parent  BODY  A  ) 


create  copy 
of  face  and 
add  it  to 
BODYC 


remove  copy  of 
this  face  from 
BODYC 


create  copy 
of  face  and 
add  it  to 
BODYC 


Yes 


remove  copy  of 
this  face  from 
BODYC 


find  the  new  intersection 
edge  on  this  face  and 
subdivide  this  face  to  create 
new  subdivision  face  and  add 
it  to  BODY  C 


get  the  subdivision  face 
correspondmg  to  this  face 
and  modify  this  face  using 
new  intersection  edge  on 
this  face 


Figure  6-21:  Flowchart  of  procedure  face_processor 
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Another  small  function  classify Jace  is  used  to  test  if  a  non-intersected  face 
belongs  to  resultant  BODY.  Function  classify  face  returns  an  integer  kl,  which  indicates 
the  given  face  belongs  to  resuhed  BODY  when  its  value  is  1  and  otherwise  if  the  value  is 
0. 

Implementation  of  Updating  Algorithm 

Updating  algorithm  introduced  in  figure  5-12  is  implemented  as  a  recursive 
member  function  regenerate( )  of  FEATURE  class.  The  implementation  can  be  described 
using  flowchart  shown  in  figure  6-22.  To  take  advantage  of  localizing  geometry 
modification  introduced  in  chapter  5,  a  modify  flag  is  attached  to  each  BODY  on 
BODY  list  of  a  feature.  Initially,  modify  flags  of  one  feature  are  all  0.  In  the  process  of 
dimension  change,  modify  flags  of  a  feature  is  set  to  1  based  on  following  rules: 

1 .  if  the  dimensions  of  the  primitive  are  changed,  all  modify  flags  must  be  set 
tol; 

2.  if  any  constraint  dimension  is  changed,  modify  flag  of  current  BODY  must 
set  to  1 ; 

3.  if  modify  flag  of  current  BODY  is  1,  and  this  feature  is  the  /-th  child  of  its 

parent  feature,  the  modify  flags  corresponding  to  ij+l),  (i+2)  (n+J)-th 

BODYs  on  the  parent  feature's  BODY  list  must  set  to  1,  where  n  is  the  total 
number  of  children  of  parent  feature. 

A  moved  flag  is  attached  to  each  geometric  entity  involved  in  the  solid  construction 
process  and  is  initially  0.  The  moved  flag  could  be  set  to  1  in  the  process  of  primitive 
modification,  constraint  imposition  or  intersection  re-computation  as  long  as  the 
geometric  entity  is  changed  or  moved.  Based  on  moved  flag,  a  small  flmction 
reposition_need( )  is  used  to  detect  if  the  position  constraints  of  this  feature  need  to  be  re- 
evaluated. 
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delete  Attach  operator's  current  ConstraintSet  and  initialize  a  new  one 


Attach  operator     evaluate  constraints 


Attach  operator  calls  modify_bln_union( )  or  other  modifying  boolean  function 


Figure  6-22:  Implementation  of  updating  algorithm 
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Stand-alone  Constraint  Imposition  Package 

The  sequential  constraint  imposition  algorithm  introduced  in  chapter  4  is  partially 
implemented  as  a  stand-alone  constraint  imposition  package  CONSPACK.  This  package 
can  be  used  to  position  one  rigid  body  with  respect  to  another  by  specifying  spatial 
relations  between  entities  of  these  two  rigid  bodies.  The  primary  purpose  of  developing 
this  package  is  for  positioning  the  geometry  of  one  feature  with  respect  to  its  parent  in  the 
FBM  system.  However,  the  package  is  written  in  such  a  way  that  it  can  be  easily 
integrated  to  other  applications  such  as  assembly  modeling.  The  whole  package  is 
composed  of  two  C++  classes,  which  are  PositionConstraint  and  ConstraintSet,  including 
about  1300  lines  of  standard  C++  source  code. 


Definition  of  PositionConstraint  Class 

PositionConstraint  is  a  base  class  used  to  derive  different  position  constraint 
types.  This  class  can  be  illustrated  using  object  modeling  diagram  shown  in  figure  6-23. 


ConstraintSet 


O- 


PositionConstraint 


pconstrainttype 

basevector 

basejpoint 

baseentitytype 

dependentvector 

dependent_point 

dependententitytype 

value 

solution  select  flag 

applytransform  function 

virtual  constraint  evaluation  function 


Figure  6-23:  Object  modeling  diagram  of  PositionConstraint  class 
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Position  constraint  type  represented  by  PositionConstraint  class  is  specified  using 
an  integer  pconstraint  type.  The  position  constraint  types  available  in  CONSPACK  are 
listed  in  table  6-1. 


Table  6-1 :  Position  constraint  types 


Position  constraint  name 

Abbreviation 

pconstrainttype 

vector  to  vector  angle 

CWA 

1 

vector  to  vector  parallel 

CWP 

2 

vertex  to  vertex  distance 

CWD 

3 

vertex  to  line  coincident 

CVLC 

4 

vertex  to  plane  distance 

CVPD 

5 

vertex  to  line  distance 

CVLD 

6 

vertex  to  vertex  coincident 

CVVC 

7 

As  mentioned  in  chapter  4,  we  use  base  object  and  dependent  object  to  refer  to 
two  rigid  bodies  to  be  positioned  by  a  position  constraint.  Entities  involved  in  the 
position  constraint  could  be  vertex  (point),  line,  plane  or  vector  attached  on  these  two 
objects.  Both  line  and  plane  can  be  represented  by  combination  of  a  point  and  a  vector. 
Therefore,  we  use  base  vector  and  base  joint,  which  are  arrays  containing  three  double 
values,  to  represent  entity  from  base  object  and  use  base  entity  type,  which  is  an  integer, 
to  specify  the  type  of  entity.  Similarly,  dependent  vector,  dependent _point  and 
dependent  entity  type  are  used  to  represent  entity  from  dependent  object  and  specify  the 
type  of  the  entity.  The  entity  types  and  the  corresponding  integer  value  of 
base  entity  type  or  dependent  entity  type  are  listed  in  table  6-2.  The  data  member  value 
represents  the  constraint  value  such  as  angle  or  distance.  An  integer  solution  select  Jlag 
is  used  to  record  information  regarding  which  solution  is  needed  in  muUiple  solution 
case.  Function  apply  transform  is  used  to  transform  the  dependent_vector  and  dependent 
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point  with  a  given  homogeneous  transformation  matrix.  The  virtual  function 
constraint  evaluation  is  defined  for  the  derived  classes  to  create  their  constraint 
evaluation  function. 


Table  6-2:  Entity  types 


Entity  type 

integer  value 

vertex  (point) 

1 

line 

2 

plane 

3 

vector 

4 

Seven  classes  corresponding  to  the  constraint  types  in  table  6-1,  are  derived  from 
the  PositionConstraint  class.  These  classes  inherit  all  the  members  of  PositionConstraint, 
but  have  their  own  constraint_evaluation  function.  The  inheritance  hierarchy  is  shown  in 
figure  6-24. 


PositionConstraint 


VV  A 


VV  P 


W  D 


VL  C 


VP  D 


VL  D 


W  C 


Figure  6-24:  Hierarchy  of  PositionConstraint  class 


The  PositionConstraint  class  has  a  pointer  to  an  object  of  ConstraintSet  it  belongs 
to.  Therefore,  one  individual  position  constraint  knows  information  about  previously 
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imposed  position  constraint.  Based  on  this  information  and  using  methods  introduced  in 
chapter  4,  the  particular  constraint  evaluation  flinction  of  a  position  constraint  computes 
and  returns  a  transformation  matrix,  which  can  be  used  to  impose  this  constraint  without 
violating  previously  imposed  position  constraints. 


Definition  of  ConstraintSet  Class 

ConstraintSet  is  an  interface  class,  by  which  user  communicates  with 
CONSPACK  and  defines  and  imposes  position  constraints.  In  this  process,  user  does  not 
have  to  be  aware  of  PositionConstraint  classes.  The  definition  of  ConstraintSet  class  can 
be  illustrated  using  object  modeling  diagram  shown  in  figure  6-25. 


 ConstraintSet  

fixedrotationdofs 

fixedtranslationdofs 

numrotcons 

niuntrscons 

rotationorder 

translationorder 

applytransform  function 

define  constraint  function 

computeconstrainttransf  function 

1 

PositionConstraint 


Figure  6-25:  Object  modeling  diagram  of  ConstraintSet  class 


ConstraintSet  contains  two  arrays  of  objects  of  PositionConstraint  class.  These 
two  arrays  are  rotation_cons[3],  which  contains  no  more  than  3  angular  constraints,  and 
translation_cons[3],  which  contains  no  more  than  3  distance  constraints.  The  data 
member  fixed  rotation  dofs  and  fixed Jranslation  jiofs  are  integers  recording  how  many 
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degrees  of  freedom  of  the  dependent  object  has  been  fixed.  The  data  members 
num  rot  cons  and  num  trs  cons  count  the  number  of  angular  constraints  and  distance 
constraints  this  ConstraintSet  contains.  The  position  constraints  in  ConstraintSet  can  be 
defined  in  any  order.  However,  CONSPACK  always  imposes  defined  constraints  in  the 
preferred  order  described  in  chapter  4.  This  preference  order  is  generated  and  stored  in 
ConstraintSet  as  integer  arrays  rotation_order[3]  and  translation_order[3].  For  example,  if 
the  preference  order  for  imposing  angular  constraints  is  {rotation_cons[l], 
rotation_cons[2],  rotation_cons[0]},  then  integers  stored  in  the  array  rotation  order 
should  be  1,  2,  0  respectively. 

The  procedure  of  using  CONSPACK  can  be  described  as: 


An  entity  stored  in  PositionConstraint  is  actually  a  simulation  of  the  entity  on  real 
object.  After  defining  a  position  constraint,  the  position  and  orientation  of  the  simulated 
dependent  entity  stored  in  PositionConstraint  is  exactly  same  as  the  entity  on  real  object. 
However,  the  dependent  object  will  be  transformed  to  new  position  or  orientation  by 
using  the  computed  transformation  matrix.  At  this  time,  the  simulated  dependent  entity 
stored  can  not  reflect  the  real  position  and  orientation  of  the  corresponding  entity  on 
dependent  object  anymore.  To  solve  this  problem,  we  defined  a  private  function  called 
apply  transform  to  transform  the  simulated  dependent  entities  such  that  the  consistency 
between  simulated  entity  and  corresponding  entity  on  the  real  object  is  maintained. 


define  a  position 
constraint 


evaluate  this 
constraint  to  get  a 
transformation  matri.x 


transform  the 
^  dependent  object 


with  the  matrix 


Figure  6-26:  Procedure  of  using  CONSPACK 
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The  functions  define  constraint  and  compute  constraint  trans/  are  the  only  two 
functions  in  CONSPACK  that  user  needs  to  know.  By  using  define  constraint  function, 
user  can  define  a  particular  position  constraint  and  add  it  to  ConstraintSet.  The 
transformation  matrix  for  imposing  this  defined  position  constraint  can  be  obtained 
immediately  by  using  compute  constraint  transf  after  the  constraint  definition. 


Directions  for  Using  CONSPACK 

Using  CONSPACK  is  quite  simple  since  user  only  needs  to  know  how  to  use 
function_constraint  and  compute  constraint  transf.  To  position  an  object  with  respect  to 
another,  user  needs  to  first  initialize  an  object  of  class  ConstraintSet,  for  example,  SETl. 
This  can  be  done  by  writting  a  line  of  C++  code  as: 

ConstraintSet*  SETl  =  new  ConstraintSet;  (6.1) 
Now  SETl  can  be  used  as  a  tool  to  impose  position  constraint  one  by  one.  The  function 
declaration  of  define  constraint  is: 


int  define_constraint( 

double*  base_vec, 

double*  base_pos, 

int  basetype, 

double*  dep  vec, 

double*  dep_pos, 

int  deptype, 

double  value 

int       solution  flag 

int  positionconstype 

) 


Figure  6-27:  Declaration  of  function  define_constraint 


To  use  this  function,  user  needs  to  convert  the  real  entity  on  the  base  object  to  the 
first  three  input  parameters  and  the  real  entity  on  dependent  object  to  the  next  three  input 
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parameters.  In  this  parameters,  base  type  and  dep  type  can  be  determined  using  table  6- 
2.  The  constraint  value  is  stored  in  the  variable  "value".  Solution  flag  can  be  set  and 
recorded  interactively  if  the  user  feels  that  one  of  the  multiple  solutions  is  preferred.  The 
type  of  the  constraint,  position  cons  type,  can  be  select  based  on  table  6-1.  If  constraint 
definition  is  successful,  the  function  will  return  1.  Otherwise  it  returns  0.  After  all  input 
parameters  are  defined,  the  following  C++  code  can  be  used  to  obtained  the 
transformation  matrix: 


SETl  ->  defme_constraint(  ); 

double*  rslt; 

result  =  SETl->  computeconstrainttransfO; 


(6.2) 


Array  rslt  is  a  double  array  with  length  of  17.  The  obtained  homogenous 
transformation  matrix  can  be  written  as; 


7  = 


rslt[l]  rslt[2]  rslt[3]  rslt[4] 

rslt[5]  rslt[6]  rslt[7]  rslt[S] 

rslt[9]  rslt[\0]  rslt[\\]  rslt[\2] 

rslt[U]  rslt[\4]  rslt[\5]  rslt[\7] 


(6.3) 


Matrix  T  will  be  used  to  transform  the  dependent  object  and  then  the  above 
procedure  will  be  repeated  to  define  and  impose  next  position  constraint.  CONSPACK 
can  handle  under-constrained  problem  However,  if  user  tries  to  defined  an  extra 
constraint  after  all  degrees  of  freedom  are  fixed,  CONSPACK  will  give  warning 
information  and  define  constraint  function  will  return  0. 
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Table  6-3 :  Available  constraint  combinations 


Angular  constraint  combination 

Distance  constraint  combination 

CVVP  -  C WA 

CWC 

CVVA  --  CVVA  --  CWA 

CVLC  -  CVPD 

CVLC  -  CVVD 

CVLC  --  CVLD 

CVPD  -  CVPD  --  CVPD 

CVPD  -  CVPD  --  CVVD 

CVPD  -  CVPD  -  CVLD 

Different  combinations  of  position  constraints  can  be  used  to  fully  fix  all  six  DOF 
of  one  object  with  respect  to  another.  Currently,  the  implemented  combinations  are  listed 
in  table  6-3.  As  mentioned  before,  the  order  of  the  combinations  is  not  important  and  the 
user  can  define  the  constraints  in  one  combination,  in  any  order,  he  wants. 


CHAPTER  7 
RESULTS 

Overview 

In  this  chapter,  some  implementation  results  are  presented.  The  user  interface  is 
introduced  in  the  first  section  to  explain  the  menus  and  the  process  of  generating  a  solid 
model,  changing  dimension  and  updating  the  solid  model.  An  example  of  a  3D  solid 
model  generated  by  extruding  a  2D  profile  is  presented  in  second  section.  In  final  section, 
we  use  another  example  to  show  the  process  of  constructing  a  simple  3D  solid  by 
Boolean  union  operation  and  updating  the  solid  model  when  a  dimension  is  changed. 

User  Interface 

As  shown  in  figure  7-1,  the  user  interface  of  the  FBM  system  is  a  window  with  7 
pull  down  menus.  The  "File"  menu  allows  user  to  open,  save  file  and  quit  fi^om  the 
system.  "Edit"  menu  has  buttons  for  dimension  change  and  updating.  "View"  menu  is  for 
changing  background  color,  viewpoint  and  redisplaying  the  graphics  in  the  window. 
"2D_Feat"  menu  has  buttons  for  creating  2D  primitives  shown  in  figure  6-6.  By  clicking 
on  different  button  on  this  menu,  user  can  define  and  create  the  corresponding  2D 
primitive.  "3D_Feat"  menu  is  meant  for  similar  commands  for  3D  primitive  creation. 
Currently,  we  only  have  "Extrusion"  on  this  menu.  By  clicking  on  "Extrusion"  button, 
user  can  interactively  pick  an  existing  2D  primitive  as  profile  and  give  extrusion  distance 
to  extrude  it.  "Constraint"  menu  has  18  buttons,  each  of  which  represents  a  position 
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constraint.  By  click  on  a  button,  user  can  interactively  define  position  constraint. 
"Feature-Types"  menu  has  6  buttons  for  defining  and  applying  different  attach  operations 
including  Protrusion_2D,  Cut_2D,  Hole_2D,  Fillet_2D,  Chamfer_2D  and 
Boolean  union. 


The  process  of  creating  a  solid  model  using  this  user  interface  consists  of  two 
steps:  creation  of  2D  profiles  and  creation  of  3D  solid  model.  The  process  for  creating  a 
2D  profile  can  be  summarized  as  follows: 

1 .  Use  the  commands  under  the  menu  "2D_Feat"  to  create  2D  features  with 
different  primitives  interactively; 
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2.  Associate  a  2D  child  feature  with  its  2D  parent  feature  by  defining  an 
particular  attach  operation  for  the  child  feature  using  2D  attach  operators 
listed  under  the  menu  "Feature  Types"; 

3.  Define  2D  position  constraints  available  under  the  menu  "Constraints"  to 
position  the  2D  child  feature  with  respect  to  its  2D  parent  feature; 

4.  Use  "Apply  Feature"  command  under  "Feature  Types"  menu  to  create  a  new 
geometry  for  the  parent  feature  using  the  attach  operator  of  the  child  feature; 

5.  Repeating  the  step  2,  3  and  4  to  create  a  2D  profile  of  the  desired  shape. 

Assuming  that  the  needed  2D  profiles  has  been  created,  the  process  for  creafing  a 
3D  solid  model  can  be  summarized  as  following  steps: 

1 .  Extrude  the  2D  profiles  to  create  3D  primitive  features  using  the  command  in 
the  menu  "3D_Feat"; 

2.  Associate  a  3D  child  feature  with  its  3D  parent  feature  by  defining  a 
Boolean  union  attach  operation  using  "Feature_Types"menu; 

3.  Define  3D  position  constraints  using  the  commands  listed  under  the  menu 
"Constraints"  to  position  the  3D  child  feature  with  respect  to  its  3D  parent 
feature, 

4.  Apply  the  Boolean  union  attach  operation  using  the  command  under 
Feature  types  menu  to  create  a  new  B-rep  for  the  parent  feature; 

5.  Repeating  steps  2,  3  and  4  to  create  the  3D  solid  with  the  desired  shape. 
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The  dimensions  of  a  created  2D  profile  or  3D  solid  model  can  be  modified  and 
the  geometry  can  be  updated.  The  procedure  is  quite  simple: 

1 .  Click  the  "Edit  Dimensions"  button  on  the  "Edit"  menu  to  interactively  select 
dimensions  to  be  modified  and  input  the  new  values  of  these  dimensions; 

2.  Click  the  "Update"  button  on  the  "Edit"  menu  to  update  the  geometry  of  the 
2D  profile  or  3D  solid. 

Example  of  Creating  and  Updating  a  3D  Extrusion  Solid 

The  3D  extrusion  solid  shown  in  figure  7-1  is  generated  by  extruding  the  2D 
profile  shown  in  figure  7-2,  which  is  created  by  using  all  the  five  2D  attach  operators 
available  to  combine  nine  2D  primitive  features  together.  Among  these  features,  Fl  is  a 
base  feature  with  a  rectangle  shape.  F2  is  a  feature  with  Cut_2D  attach  operator  and 
rectangular  shape.  F3  and  F7  are  features  with  Protrusion_2D  attach  operator  and  they 
have  geometry  of  Arc  l  and  Arc_2  respectively.  F4  and  F8  are  features  with  Hole_2D 
attach  operator  and  circular  geometry.  F6  is  also  a  feature  with  Hole_2D  attach  operator 
and  rectangle  geometry.  F5  and  F9  are  features  with  Fillet_2D  and  Chamfer_2D  attach 
operators  respectively.  Some  of  the  dimensions  are  shown  in  the  figure  7-2. 

After  significant  dimension  changes,  the  solid  shown  in  figure  7-1  is  updated  to 
the  one  shown  in  figure  7-3.  Part  of  the  design  intent  of  this  solid  model  is  represented  by 
position  constraints  of  feature  F7,  which  are  "edge  to  edge  coincident"  between  the 
straight  edge  of  Arc_2  of  F7  and  vertical  edge  of  Fl,  and  "vertex  to  vertex  distance" 
shown  as  dl3  in  figure  7-2.  Another  part  of  the  design  intent  is  represented  by  "center  to 
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center  coincident"  constraint  between  F8  and  F7.  It  is  noticed  that  increasing  the  angle  of 
Arc_2  in  feature  F7  does  not  violate  all  these  position  constraints. 


Figure  7-2:  Features  and  dimensions  of  a  2D  profile 


Figure  7-3 :  Updated  3D  extrusion  solid 


Example  of  Creating  and  Updating  a  general  3D  solid 

In  this  section,  an  example  of  a  3D  solid  is  presented  that  was  created  by  applying 
Boolean  union  attach  operation  to  two  3D  primitive  solids.  For  clarity,  the  picture  is 
shown  in  wire-frame  form.  In  figure  7-4,  two  3D  primitives,  a  bearing  and  a  plate  with 
four  holes,  are  shown.  These  solids  are  to  be  combined  together  using  the  Boolean  union 
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operation.  In  figure  7-5,  the  bearing  has  been  correctly  positioned  with  its  parent  feature, 
the  plate.  One  of  the  position  constraints  is  the  distance  d  between  two  planar  faces  of 
these  two  primitives.  Figure  7-6  shows  the  solid  created  by  the  Boolean  union  operation. 


«sa'  fbmcxJefer 


^^^^ 


:Fiie    Edit    View    2D_reat;    3I>_Feat    ConSti-aint;    Featur  e_Types 


Figure  7-4:  Two  3D  primitives  to  be  combined  by  the  Boolean  operation 


Figure  7-6:  Resultant  solid  of  the  Boolean  union  operation 


134 


Figure  7-8:  Updating  of  3D  solid  in  case  of  topology  change. 
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In  figure  7-5  and  7-6,  distance  d  =  -15,  where  the  sign  "-"  means  child  face  is  at  a 
distance  d  in  a  direction  opposite  to  parent  face  normal.  When  changing  this  dimension  to 
-5,  the  updated  solid  is  shown  in  figure  7-7.  In  this  case,  the  topology  of  the  solid  does 
not  change.  If  we  further  change  the  dimension  to  d  =  15,  the  updated  solid  is  shown  in 
figure  7-8.  It  is  clear  that  the  updating  algorithm  successfully  updated  the  solid  with 
topology  change  without  violating  any  constraints. 


CHAPTER  8 
CONCLUSION 

Conclusions 

The  research  presented  in  this  dissertation  is  composed  of  four  parts  including 
design  intent  representation  using  extended  CSG  tree,  sequential  constraint  imposition 
method,  non-regenerative  updating  algorithm  and  implementation  of  FBM  system  and 
CONSPACK.  These  four  parts  cover  almost  all  aspects  of  dimension-driven  solid 
modeling. 

Extended  CSG  tree  is  a  complete  representation  of  design  intent  because  of  the 
additional  information  of  position  constraints  and  primitive  creation  functions.  Compared 
to  design  intent  representation  in  variational  geometry,  storing  complete  constraint  set 
and  primitive  creation  function  in  nodes  of  extended  CSG  tree  can  be  thought  of  as  a  de- 
coupling of  the  simultaneous  equations  in  variational  geometry.  Representing  the  design 
intent  using  a  tree  structure  rather  than  a  sequential  list  of  operations  makes  possible  a 
more  sophisticated  updating  algorithm  than  the  one  based  on  "roll  back". 

Although  extended  CSG  tree  is  a  new  representation  for  design  intent,  the  major 
contributions  of  this  dissertation  are  the  sequential  constraint  imposition  method  and  non- 
regenerative  updating  algorithm.  Sequential  constraint  imposition  method  is  a  method  for 
positioning  one  rigid  object  with  respect  to  another  by  sequentially  imposing  a  set  of 
position  constraint.  The  main  advantage  of  this  approach  is  that  it  does  not  require 
solving  all  the  constraints  used  in  the  model  simuhaneously.  The  constraints  between 
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each  pair  of  object  can  be  imposed  one  at  a  time  both  during  construction  and  when  the 
model  is  updated  after  changing  dimension.  Analytical  solutions  exist  for  most  constraint 
combinations  and  the  availability  of  these  analytical  solutions  enables  fast  imposition  of 
constraints  both  during  construction  and  updating  of  the  solid  model.  Using  analytical 
method  has  another  advantage  that  all  solutions  are  available  when  there  are  muhiple 
solutions.  Therefore,  the  user  can  interactively  select  the  desired  solution  based  on  the 
design  intent.  During  updating,  this  same  solution  can  be  automatically  reselected  by 
storing  information  to  distinguish  among  these  solutions. 

Another  important  contribution  of  this  dissertation  is  the  non-regenerative 
updating  algorithm  in  dimension-driven  solid  modeling.  The  algorithm  we  proposed  is 
the  first  updating  algorithm  in  this  area  that  does  not  need  entity  recreation.  As  we 
described  in  previous  chapters,  the  "rollback"  algorithm  used  in  commercial  CAD  system 
usually  involves  recreating  the  geometric  and  topological  entities,  constraint  solving  and 
Boolean  operation.  Compared  to  this  kind  of  algorithm,  the  updating  algorithm  in  figure 
5-12  is  much  more  efficient  because  of  following  reasons: 

1 .  In  response  to  dimension  changes,  the  algorithm  modifies  the  geometry  of  the 
solid  model  (B-rep)  instead  of  re-construct  it  during  the  updating  process. 
Therefore  recreation  of  topological  and  geometric  entities  can  be  avoided  to 
save  a  lot  of  time  and  memory. 

2.  Since  we  do  not  recreate  the  B-rep,  there  is  no  need  to  create  correspondence 
between  the  entities  before  and  after  the  dimension  change.  Therefore,  time 
consuming  persistent  naming  procedure  can  be  avoided. 
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3.  In  case  of  no-topology-change,  there  is  no  need  to  redo  the  Boolean  operation 
during  the  updating  process.  Even  though  intersection  re-computation  is 
needed  for  detecting  topology  change,  it  takes  much  less  time  than  Boolean 
operation. 

4.  In  case  of  topology  change,  we  modify  the  result  of  Boolean  operation 
instead  of  redoing  it.  This  procedure  is  very  similar  to  redoing  Boolean 
operation  except  that  topological  entity  recreation  is  never  needed.  In  the 
worst  case,  the  efficiency  is  same  as  that  of  Boolean  operation.  However,  the 
advantages  listed  in  1  and  2  are  still  applicable. 

5.  Non-linear  simultaneous  equation  solving  in  primitive  recreation  can  be 
avoided  by  using  primitive  modification  instead  of  primitive  recreation. 

6.  Non-linear  simultaneous  equation  solving  in  primitive  reposition  can  be 
avoided  by  using  sequential  constraint  imposition  method  proposed  in  this 
dissertation. 

7.  The  non-regenerative  updating  algorithm  is  further  improved  by  ensuring  that 
geometric  modification  is  localized.  With  the  localization  approach,  a  lot  of 
steps  in  the  non-regenerative  updating  algorithm  can  be  skipped  to  improve 
updating  efficiency.  The  design  intent  representation  using  extended  CSG  tree 
contributes  a  lot  to  this  localization  approach,  which  is  not  possible  when 
design  intent  is  only  recorded  as  a  sequential  list  of  model  construction 
operations. 
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In  this  dissertation,  the  proposed  approaches  and  algorithms  have  been 
implemented  in  a  feature-based  dimension-driven  solid  modeling  system  (FBM)  and  a 
stand-alone  constraint  imposition  package  (CONSPACK).  The  main  purpose  of 
implementing  the  FBM  system  is  to  test  the  proposed  methods.  This  system  can  go 
through  the  whole  process  of  dimension-driven  solid  modeling  although  it  is  not  as  robust 
as  commercial  system.  The  software  package  CONSPACK  can  be  used  not  only  for  the 
purpose  of  this  dissertation,  but  also  for  other  applications  related  to  rigid  body 
positioning.  It  is  designed  with  object-oriented  data  structure  such  that  it  can  be  easily 
integrated  into  other  systems. 

During  the  implementation  of  FBM  system,  we  experienced  some  difficulties 
when  using  ACIS  as  the  geometric  modeler.  Since  all  existing  CAD  systems  are  using 
regenerative  updating  algorithms,  ACIS  does  not  allow  users  to  easily  modify  the  B-rep 
and  the  resuh  of  Boolean  operation.  Some  information  required  for  implementing  the 
updating  algorithm  presented  in  this  dissertation  is  lost  or  not  recorded  during  ACIS 
Boolean  operations.  Therefore,  we  implement  our  own  Boolean  operation. 

Future  Directions 

Dimension-driven  solid  modeling  is  becoming  an  active  research  topic  in  CAD 
area  because  it  has  direct  application  in  industry.  Further  research  on  this  topic  can  make 
the  contribution  of  this  dissertation  more  valuable  in  industrial  application.  The  possible 
extensions  for  the  research  presented  in  this  dissertation  are: 

1.  Extending  the  sequential  constraint  imposition  method  to  handle  constraint 
sets  that  do  not  include  a  complete  angular  constraint  set.  For  example,  it  may 
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be  necessary  to  use  three  WC  constraints  to  fully  fix  3  rotation  and  3 
translation  DOFs.  How  to  convert  this  type  of  constraint  set  to  the  complete 
constraint  set  defined  in  this  dissertation  would  be  very  valuable  in  practical 
applications 

2.  In  the  updating  algorithm  presented,  a  better  topology  change  detection 
algorithm  would  significantly  improve  the  efficiency  of  this  algorithm.  The 
detection  method  using  intersection  re-computation  is  very  time  consuming. 

3.  Although  the  implementation  is  enough  for  testing  the  methods  and 
algorithms  presented  in  this  dissertation,  the  CONSPACK  and  FBM  are 
incomplete  from  industrial  point  of  view.  Currently  only  analytical  methods 
are  implemented  in  CONSPACK.  Therefore,  some  combinations  of  distance 
constraints  are  not  available  at  this  time.  Implementing  numerical  methods  in 
CONSPACK  is  a  valuable  extension  of  this  research.  We  only  implemented 
the  Boolean  union  operation  in  FBM.  Other  Boolean  operations  and  other 
primitive  creation  function  could  be  added  to  this  system  in  future. 


LIST  OF  REFERENCES 


ACIS  Manuals,  Spatial  Technology  Inc,  1996. 

Aldefeld,  B.,  "Variation  of  geometries  based  on  a  geometric-reasoning  method", 
Computer-Aided  Design,  Vol.  20,  No.  3,  pp.  117-126,  1988. 

Ambler,  A.  P.  and  Popplestone,  R.  J.,  "Inferring  the  positions  of  bodies  from  specified 
spatial  relationships".  Artificial  Intelligence,  Vol.  6,  pp.  157-174,  1975. 

Bettig,  B.  and  Shah,  J.,  "Systematic  derivation  of  a  standard  set  of  declarative  constraints 
for  shape  definition  in  CAD",  Proceedings  of  1999  Design  Engineering 
Technical  Conferences,  September  12-15,  Las  Vegas,  Nevada,  1999. 

Bouma,  W.,  Fudos,  I.,  Hoffmann,  C,  Cai,  J.  and  Paige,  R.,  "Geometric  constraint  solver", 
Computer-Aided  Design,  Vol.  27,  No.  6,  pp.  487-501,  1995. 

Brain,  M.,  "MOTEF  programming:  The  essentials  ...  and  more".  Digital  Press,  Newton, 
MA,  1992. 

Buchanan,  S.  A.  and  De  Pennington  A.,  "Constraint  definition  system:  a  computer- 
algebra  based  approach  to  solving  geometric-constraint  problems", 
Computer-Aided  Design,  Vol.  25,  No.  12,  pp.741-750,  1993. 

Capoyleas  V.,  Chen  X.  and  Hoffman  C.  M.,  "Generic  naming  in  generative  constraint- 
base  design",  Computer-Aided  Design,  Vol.  28,  No.  1,  pp.  17-26,  1996. 

Chen,  X.  and  Hoffmann,  C.  M.,  "Design  compilation  of  feature-based  and  constraint- 
based  CAD",  Proceedings  of  the  symposium  on  Solid  Modeling,  ACM, 
NY,  p  13-19,  1995a. 

Chen,  X.  and  Hoffmann,  C.  M.,  "Towards  feature  attachment",  Computer-Aided  Design, 
Vol.  27,  No.  9,  pp.  695-702,  1995b. 

Chen,  X.  and  Hoffmann,  C.  M.,  "On  editability  of  feature-based  design",  Computer- 
Aided  Design,  Vol.  27,  No.  12,  pp.  905-914,  1995c. 

Chung,  J.  C.  H.,  Patel,  D.  R.  and  Cook,  R.  L.,  "Feature-based  modeling  for  mechanical 
design".  Computer  &  Graphics,  Vol.  14,  No.  2,  pp.  189-199,  1990. 

Deitel,  H.  M.  and  Deitel,  P.  J.,  "C++  How  to  program".  Prentice  Hall,  1994. 


141 


142 


Duan,  W.,  Zhou,  J.  and  Lai,  K.,  "FSMT:  a  feature  solid-modelling  tool  for  feature-based 
design  and  manufacture",  Computer-Aided  Design,  Vol.  25,  No.  1,  pp.  29- 
38,  1993. 

Fudos,  I.  and  Hoffmann,  C,  "A  Graph-Constructive  Approach  to  Solving  Systems  of 
Geometric  Constraints",  ACM  Transactions  on  Graphics,  Vol.  16,  No.  2, 
pp.  179-216,  1997. 

Gao,  X.-S.  and  Chou,  S.-C,  "Solving  geometric  constraint  systems,  part  I.  A  global 
propagation  approach",  Computer-Aided  Design,  Vol.  30,  No.  1,  pp.  47- 
54,  1998a. 

Gao,  X.-S.  and  Chou,  S.-C,  "Solving  geometric  constraint  systems,  part  U.  A  symbolic 
approach  and  decision  of  rc-constructibility",  Computer-Aided  Design, 
Vol.  30,  No.  2,  pp.  47-54,  1998b. 

Ge,  J.-X.,  Chou,  S.-C.  and  Gao,  X.-S.,  "Geometric  constraint  satisfaction  using 
optimization  method",  Computer-Aided  Design,  Vol.  31,  No.  14,  pp.  867- 
879,  1999. 

Gossard,  D.  C,  Zuffante,  R.  P.,  Sakurai,  H.,  "Representing  dimensions,  tolerance,  and 
features  in  MCAE  systems",  IEEE  Computer  Graphics  and  Applications, 
pp.  51-59,  1988. 

Han,  J.  and  Requicha,  A.  A.  G.,  "Integration  of  feature  based  design  and  feature 
recognition",  Computer-Aided  Design,  Vol.  29,  No.  5,  pp.  393-403,  1997. 

Hoffmann,  C.  M.  and  Juan  R.,  "Erep  -  An  editable,  high-level  representation  for 
geometric  design  and  analysis",  IFIP  Transactions,  Geometric  Modeling 
for  Product  Realization,  P.  R.  Wilson,  M.J.  Wozny  and  M.J.  Pratt,  eds., 
Elsevier  Science  Publishers  (North-Holland),  pp.  129-164,  1992. 

Hoffmann,  C.  M,  and  Joan-Arinyo,  R.,  "On  user-defined  features",  Computer-Aided 
Design,  Vol.  30,  No.  5,  pp.  321-332,  1998. 

Kramer,  G  A.,  "Using  degrees  of  freedom  analysis  to  solve  geometric  constraint 
systems",  ACM  Symposium  Foundations  of  Solid  Modeling,  ACM,  New 
York,  371-378,  1991. 

Kramer,  G.  A.,  "Solving  Geometric  Constraints  Systems:  A  Case  Study  in  Kinematics", 
MIT  Press,  Cambridge,  MA,  1992. 

Kumar,  A.  V.  and  Yu,  L.,  "Sequential  Constraint  Imposition  for  Dimension-Driven  Solid 
Models",  Accepted  by  Computer-Aided  Design. 


143 


Kripac,  J.,  "A  mechanism  for  persistently  naming  topological  entities  in  history-based 
parametric  solid  models",  Proceedings  of  Solid  Modeling  '95  symposium, 
pp.  21-30,  1995. 

Kripac  J.,  "A  mechanism  for  persistently  naming  topological  entities  in  history-based 
parametric  solid  models",  Computer-Aided  Design,  Vol.  29,  No.  2,  pp. 
113-122,  1997. 

Laakko,  T.  and  Mantyla,  M.,  "Feature  modelling  by  incremental  feature  recognition". 
Computer- Aided  Design,  Vol.  25,  No.  8,  pp.  479-492,  1993. 

Lee,  K.  and  Gossard,  D.  C,  "A  hierarchical  data  structure  for  representing  assemblies: 
part  1",  Computer-Aided  Design,  Vol.  17,  No.  1,  pp.  15-19,  1985. 

Lee,  K.  and  Andrews,  G.,  "Inference  of  the  positions  of  components  in  an  assembly:  part 
2",  Computer-Aided  Design,  Vol.  17,  No.  1,  pp.  20-24,  1985. 

Lee,  K.,  "Principles  of  CAD/CAM/CAE  systems",  Addison-Wesley,  1999. 

Light,  R.  and  Gossard,  D.  C,  "Modification  of  geometric  models  using  variational 
geometry",  Computer-Aided  Design,  Vol.  14,  No.  2,  pp.  209-214,  1982. 

Light,  R.  and  Gossard,  D.  C,  "Variational  Geometry:  A  new  method  for  modifying  part 
geometry  for  finite  element  analysis".  Computer  &  Structure,  pp.  903-909, 
June,  1983. 

Lin,  V.  C,  Gossard,  D.  C.  and  R.  A.  Light,  "Variational  Geometry  in  Computer-Aided 
Design",  Computer  Graphics,  ACM,  Vol.  15,  No.  3,  pp.  171-177,  1981. 

Mantyla,  M.,  "An  introduction  to  solid  modeling",  Computer  Science  Press,  Rockville, 
Maryland,  1988. 

Martino,  T.  D.,  Falcidieno,  B.,  Giannini,  P.,  Hassinger,  S.  and  Ovtcharova,  J.,  "Feature- 
based  modelling  by  integrating  design  and  recognition  approaches", 
Computer-Aided  Design,  Vol.  26,  No.  8,  pp.  646-653,  1994. 

Mortenson,  M.  E.,  "Geometric  Modeling",  John  Wiley  &.  Sons,  1985. 

Owen,  J.  C,  "Algebraic  Solution  for  geometry  from  dimensional  constraints",  ACM 
Symposium  Foundations  of  Solid  Modeling,  ACM,  New  York,  397-407, 
1991. 


Parametric  Technology,  Pro/Engineer  Users  Guide,  Parametric  Technolgy  Corp., 
Waltham,  MA,  1989. 


144 


Peters,  T.  J.,  Demurjian,  S.  A.  and  Ting,  T.  C,  "Feature-based  modeling  by  object- 
oriented  design  with  propagation".  Proceedings  of  International 
Conference  on  Data  and  Knowledge  Systems  for  Manufacturing  and 
Engineering,  1994. 

Raghothama,  S.  and  Shapiro,  V.,  "Consistent  updates  in  dual  representation  systems". 
Proceedings  of  5th  ACM  Symposium  on  Solid  Modeling  and 
Applications,  Ann  Arbor,  Michigan,  1999. 

Raghothama,  S.  and  Shapiro,  V.,  "Boundary  representation  deformation  in  parametric 
solid  modeling",  ACM  transactions  on  Graphics,  Vol.  17,  No,  4,  pp.  259- 
286,  1998. 

Rambaugh  J.,  Blaha  M,  Premerlani  W,  Eddy  F  and  Lorensen  W.,  "Object-oriented 
modeling  and  design".  Prentice  Hall,  1991. 

Requicha,  A.  A.  G.  and  Voelcker  H.  B.,  "Solid  modeling  -  a  historical  summary  and 
contemporary  assessment",  IEEE  Computer  Graphics  and  Applications, 
Vol.  2,  No.  2,  p.  9-24,  1982. 

Requicha,  A.  A.  G.  and  Voelcker  H.  B.,  "Boolean  operations  in  solid  modeling: 
Boundary  evaluation  and  merging  algorithms".  Proceedings  of  the  IEEE, 
Vol.  73,  No.  l,pp.  30-44,  1985. 

Serrano,  D.,  "Constraint  Management  in  Conceptual  Design",  Doctor  of  Science  thesis. 
Department  of  Mechanical  Engineering,  MIT,  Cambridge,  MA,  1987. 

Serrano,  D.,  "Automatic  dimensioning  in  design  for  manufacturing",  ACM  Symposium 
Foundations  of  Solid  Modeling,  ACM,  New  York,  379-386,  1991 

Shah,  J.  and  Mantyla,  "Parametric  and  feature-based  CAD/CAM:  Concepts,  Techniques 
and  Applications",  John  Wiley  &  Sons,  1995. 

Shapiro,  V,  and  Vossler,  D.  L.,  "What  is  a  parametric  family  of  solids?".  Proceedings  of 
the  Third  ACM/IEEE  Symposium  on  Solid  Modeling  and  Applications. 
May  17-19,  1995. 

Shimizu,  S.  and  Numao,  M.,  "Constraint-based  design  for  3D  shapes".  Artificial 
Intelligence,  Vol.  91,  pp.  51-69,  1997. 

Tilove,  R.  B.,  "Set  membership  classification:  A  unified  approach  to  geometric 
intersection  problems",  IEEE  Transactions  on  Computers,  Vol.  C-29,  No. 
10,  pp.  874-883,  1980. 

Zeid,  I.,  "CAD/CAM  theory  and  practice",  McGraw-Hill,  1991. 


BIOGRAPHICAL  SKETCH 
Lichao  Yu  was  born  on  March  13,  1966,  in  Jiangxi,  P.R.China.  He  earned  his 
Bachelor  of  Science  degree  in  mechanical  engineering  from  Huazhong  University  of 
Science  and  Technology,  in  1986,  and  Master  of  Science  degree  in  mechanical 
engineering  specializing  in  finite  element  analysis  from  Huazhong  University  of  Science 
and  Technology  in  1989.  He  went  to  the  University  of  Florida  in  1996  to  study  for  his 
doctoral  degree.  He  will  earn  the  degree  of  Doctor  of  Philosophy  in  August  of  2000. 


145 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Ashok  V.  Kumar,  Chairman 
Assistant  Professor  of  Mechanical 
Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Carl  D.  Crane  m 

Professor  of  Mechanical  Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Jo 
Pro 


C.  Ziegert 
jssor  of  Mechanic 


Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Gloria  XWiens 

Associate  Professor  of  Mechanical 


Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 

Baba  C.  Vemuri 

Professor  of  Computer  and  hiformation 
Science  and  Engineering 


This  dissertation  was  submitted  to  the  Graduate  Faculty  of  the  College  of 
Engineering  and  to  the  Graduate  School  and  was  accepted  as  partial  fulfillment  of  the 
requirements  for  the  degree  of  Doctor  of  Philosophy. 


August  2000 


M.  Jack  Ohanian 

Dean,  College  of  Engineering 


Winfred  M.  Phillips 
Dean,  Graduate  School 


